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engaged in space flight and related efforts. 

Reprints of these articles may be obtained from the NASA 
Library, MSFC, Huntsville, Alabama 
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GUIDANCE AND CONTROL 
SYSTEMS AND COMPONENTS 


N65 23782 


A. INTRODUCTION TO THE KALMAN-SCHMIDT FILTER 
by John W. Kenelly and Kenneth S. Williamson 


This article provides an insight into the application 
of the Kalman-Schmidt filter to navigation for technical 
personnel with limited background in probability theory. 
The scope of the article includes the statistical deriva- 
tion of the Kalman-Schmidt filter within the context of 
spaceflight navigation. The filter may be used to locate 
a space vehicle and to determine its velocity vector from 
data containing random errors. The filter is superior to 
least squares because each observation datum does not 
have to determine the position of the vehicle. 

AuTi&R- 


is not possible. However, with reasonable instru- 
mentation, the Kalman-Schmidt filter can give results 
that will minimize errors. 

The mathematical background for the Kalman- 
Schmidt filter was performed by Kalman (Ref. 1). 
Schmidt was the first to notice that Kalman ! s results 
could be applied to space navigation (Ref. 2). Varied 
applications have been studied by the National Aero- 
nautics and Space Administration. Because the origi- 
nal paper by Kalman is highly theoretical, several 
papers have been written for the purpose of bringing 
the results as applied to navigation direct to the engi- 
neer. This article is one of these papers and includes 
necessary navigational background. The reader needs 
only a slight acquaintance with matrices and statistics. 
This article should provide insight into the data re- 
duction problems of navigation and into the Kalman- 
Schmidt filter. This insight will be of value in 
programing and studying a computer simulation of the 
navigation problem. 


INTRODUCTION 


NOMENCLATURE 


Space navigation includes the determination of t 
position and velocity of a vehicle by onboard or earth n 
based measurements of range, range rate, or various x 
combinations of star -planet, sun -planet angular posi- 1 

tion fixes. The instruments used in making the meas- 
urements are not perfect and will produce errors. 

One method of handling the errors (i.e. , data reduc- 
tion of the measurement) is the Kalman-Schmidt filter 
theory, which is a mathematical method of processing y 
data to remove errors in measurement. Since the 
filter* only receives as input inaccurate measure- 
ments, precise determination of position and velocity 


time of rfi 1 set of measurements 

state (column) vector at t^. X n 
is a perturbation vector; i. e. , 
is the difference between the actual 
trajectory and the standard or ref- 
erence trajectory at time t^. 

measurement (column) vector at 
t^. Y n is also a perturbation vec- 
tor; i. e. , Y n is the difference be- 
tween the measurables on the tra- 
jectory and the reference. 


*In this article, the term "filter" denotes a mathema- Y^ m 
tical method or computer device that receives noisy 
measurements (i.e., measurements with superimposed 
error) and puts out a more realistic approximation of 
the trajectory. 


the value of the perturbation vec- 
tor Y n of measurements that are 
indicated by the instrumentation ; 
i.e. , Y n * m is the perturbation 
of the measurement of Y n . 
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U Y n * m = Y n + Ujji i. e. , u n is a zero 

n mean random noise (column) vec- 

tor. u n is the error made by the 
measuring instruments. 

Q Q n = E(u n u n T); i. e. , Q n is the 

covariance of the noise u n . Q n is 
given information, and the noise u n 
at time t^ is given to be uncor- 
related with the noise at any time 
not tn, i. e. , any different time 

n ^n+l = ^n+l,n X n» L e * » ^n+l,n is 

* the state transition matrix, men- 

tioned previously, which translates 
Xjj into X n+ i. This matrix can be 
precomputed before a mission. 

M n Y n = M n^Si or = M n 1 Y n* M n 

is the geometry mapping matrix, 

mentioned previously and precom- 
putable, which maps the measure- 
ment perturbation vector Y n into 
the state perturbation vector x n . 

- * 

X the best estimate of X n as deter- 

mined by the filter to be derived. 
X n * is a filter approximation of 
X n and is based on the given sta- 
tistics and all measurement vec- 

5i< m 

tors Y n taken prior to time t n . 
E(X n * - X^ = 0. (E is the ex- 
pected value operator and is a sta- 
tistical average. E(X n * - X n ) = 0 
states that the average value of the 
approximation X n * of X n is X n or 
equivalently that the difference is 
zero. ) 

the projection of the perturbation 
vector X n ‘ of time into at 

timet n+ ^; i. e. , X n 4.j^ - ^ n Xn • 

Since E is a linear operator, and 
^n+1 n is a constant matrix, 
<Pn+l,n can Pass directly over the 
expected value operator E. E (X^ - 
*n> = °. ® (Xjj* ) - E(X n ) = 0, 
E(X n * ) = X n , <p n+1 E(Xu* ) = 
^n+1, n X n* and finally E(X n+1 ) = 
X^p from assumptions in the 
definition of X n * above. 


Y Y n = M n X^ ; Y n is an estimate of 

Y *m 
n 

E(Y n ) = E(M n X n * ) = M n E(Xn* ) 
E(Y n ) = M nXn 

E ( Y n ) = Y n or the average value 
of the perturbation measurement 
vector Y n . The difference between 

Y n and Y n is that Y n also contains 
instrument noise, where Y n does 
not. 

<5X error in X 

n n 

^Xn-Xn 

6Y error in Y„ 

n n 

— m * m 

6Y n =Y n- Y n = M nXn 

. , „ * m v . 

Also, E(Y n ) = Y n and since 

E(Y n )=Y n 

then E ( <5¥ n ) = E ( Y n ) - E(Y n * m ) = 
Y n - Y n =0 


n+1 

the best estimate (or approxima- 
tion) of <5^, <5X^ +1 = 3^ - X£ 

p* 

n 

given statistics; i.e., the covari- 
ance matrix of errors in the state 
vector X n 


p;=E[(<-Xn) (Xj-xj 1 ] 


DERIVATION OF KALMAN- SCHMIDT FILTER 

The Kalman- Schmidt filter as applied to a typical 
space navigation problem is derived in this section. 
In particular that part of a mission in which position 
and velocity are determined is examined. Deviations 
from a reference trajectory are assumed where the 
reference trajectory is calculated by a computer be- 
fore the mission. The advantage of using a reference 
trajectory lies chiefly in the capability of precalculat- 
ing two kinds of transition (or mapping) matrices. 
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The two matrices are the state-transition matrix 

<p . and the geometry mapping matrix M n . The 
n+ 1 ,n n 

geometry mapping matrix maps the measurements in- 
to the state (position and velocity) , and the state- 
transition matrix translates the state at time into 

the state at a later time t . . 

n+i 

The state-transition and geometry mapping ma- 
trices are approximations based on the assumption 
that the trajectory is not perturbed far from the stand- 
ard. In general, the matrices will not work for large 
deviations. 

A theorem from statistics (Ref. 3) relates some 
of the vectors by <5X^ +i = [£( 6 X ^1 <5Y n T ) ] 

tE(<5Y n 6Y n T )]“ 1 8Y n . 

This statement is the core of the derivation of the 
filter and it only remains to rearrange these terms 
into more useful relations. 

Note that the matrices <p n4 .^ n and translate 
and map any perturbation vector near the reference; 
i. e. , one definition says ^ n +i = <P n +i ^ since is 
sufficiently close to the reference. 

Figures 1 and 2 illustrate the action of a filter. 
The filter receives continuous or discrete noisy meas- 
urements of a trajectory and then puts out an estimated 
trajectory. The smooth lines represent the actual 
trajectory and are included for convenience since a 
filter would have no knowledge of the actual trajectory. 

The trajectory and reference are shown along with 
the perturbation state vector X n in Figure 3. The 
measurement perturbation vectors Y n and Y n * m are 
given in the equations. The angles shown are included 
for convenience. The star -planet etc. fixes are not 
sufficient to determine a velocity component for the 
state vector X^ in one measurement. 

Note that there is a one-to-one correspondence 
between X n and Y n and that Y n = M n X n , and it is not 
necessary to clutter the illustrations unnecessarily 
with measurement perturbation vectors’ Again, as- 
suming there is a measurement vector Y n which 
determines a state vector X n , then any illustration of 
state vectors implies a corresponding illustration of 
measurement vectors. 



FIGURE 1. FILTERING OF CONTINUOUS NOISE 




FIGURE 3. MEASUREMENTS ON TRAJECTORY 
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TABLE 1. PERTINENT DETAILS* 


T 1. Definition - A random variable is statistical 
and a constant variable is nonstatistical. 

T 2. E(A + B) = E(A) + E(B); A and B are random 

T 3. E(C) = C; C is a constant 

T 4. E ( CA) = CE ( A) ; C is constant and A is random 

T 5. E(AB) = 0, if A and B are uncorrelated 

T T T 
T 6. (AB) = B A 

T 7. AX = B => X = A -1 B 

T 8. (A + B) T = A T +B T , (DA + DB) T = (A + B) T D T 
* m 

T9. Y fl =Y n + u n 

T10. Q n =E(u n u n T ), 

Til. X . = cp . , n X n 
n+1 ^n+1 


T 12. Y n = M nXn 

T 13. E(X* - X ) = 0 
n n' 

T 14. E(S n+1 , = X o+1 

T 15. Y = M X* ; E(Y ) = M E(X* ) 
n n n n n ' n ' 


= M X ; E(Y ) = Y 
n n n' n 


T 16. 6X = X - X 
n n n 

. — — ^ m ❖ m 

T 17. 6Y = Y - Y = M X 
n n n n n 


- Y * m , E(6T )= 0 
n n' 


T 18. 6X* = X - X* 
n n n 

T 19. P* - E[(X^ -X„) (X*- V T 1 


* Numbers on the left are referred to in the following equations. 


error in 6X^ + ^ is such that variance in the estimate 
of any quantity ip = A^5X n+ ^ is minimized, or 

= <*- n T ] 

'J' jjc 

is a minimum ( ip = A <5 X n+ ^ ) . 

Table 1 show| all usable material in the deriva- 
tion. Note that P n and Q n must appear in the result 
because they are the only given statistics. 

The first step in the derivation consists in solving 
and substituting for the three factors in the equation 
from the theorem. 

Theorem, 


Part 1. 

E(6Y <5Y T ) = E[(Y - Y * m ) (Y -Y* m ) T | (T 17) 
n n ' lv n n ' ' n n ;J ' 


= E[(Y - Y - u )<Y - Y - u ) ] (T 9) 
n n n n n n 


= E[{M(X -X )-u}{(X' - XpM -u 1 }] 

1 1 n n n' n J L n ir n n ’ J 

{T 12,15,8) 

= M E[(X* - X ) (X* -X ) T 1 M T + E(u u T ) 
n n n n n J n v n n ' 

- E[M (X*- X ) u T ] - E[u (X*- X ) T ], (T 4,6) 
1 n n n n J 1 n n n J * 7 


Since E[ X 1 - X ] = 0, this reduces to 
1 n n J 


SX* + 1 = [E(« 


n+1 


6Y T )][E(6Y 6Y T )] _1 6Y 
n 1 n n n 


= M P * M T + Q . 
n n n n 


(T 10,19,5) 
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Part 2. 

E ( <5X . . 6Y T ) = E[(X . - X (Y - Y * m ) T ] 

' n+i n 7 l ' n+1 n+1 7 ' n n 7 J 

(T 16,17) 


Part 3. Parts 1 and 2 give substitution for the theorem 
and 

6X* = [<» . P* M T 1 [M P* M T + Q ] _1 6Y 

n+i LY n+i,n n n J L n n n n J n 


x , - X* = [<p P* M T 1[M P* M T + Q ] _1 
n+i n+1 lY n+l,n n n n n n n J 


w * * m v 

x (M X - Y ). 
n n n 


(T 18,17,12) 


Therefore 

X* = <» . X* - <p . P*mJ(M P^mJ+QJ " 1 

n+1 r n+l, n n r n+l,n n n n n n ^n' 

* * m 

x ( M X^‘ - Y ) 
n n n 

which is the prediction of the filter in terms of known 
matrices. 

Now some substitutions are made and is iter- 
ated into Pn+1- The covariance matrix P^+i is up- 
dated from P^ because each observation changes the 
covariance matrix. 


Rearranging X^ +1 , 

X* = [<p . . - cp P* M T ( M P* M T 

n+1 1 n+l,n n+l,n n n 1 n n n 




_ -1 * m 

+ Q ] Y 
n n 


= <p . 4 E[(X * - X ){(X* - 
^n+l,n 1 n n l ' n 

X) T M T -u T }] 
n 7 n n JJ 

Let 


(T 4,11,8,9) 

ii 

* d 
<1 

= (p ^ E[(X * - X ) (X* - 
r n+l,n 1 ' n ir n 

T T 

- x n m 1 

n 7 J n 

then 

- w , E[(X * - X ) u T l 
r n+l,n 1 n n 7 n J 

(T 5,8) 

n+1 

E(6X „ 6 Y T ) = (p . P * M T 
' n+1 n ’ ^n+l.n n n 

(T 19) 

Let 

since E[ X - X 1 = 0. 
1 n n J 


* 

‘Vi, 


A 5 " = <p . P* M T (M P* M T +Q ) 1 

n T Ti-U1 v-» ri v-k ' yr n Tl 


i* M ] X* + A*Y* m . 
n n J n n n 


then 

X* , = <p*. 4 X* + A* Y * m . 
n+1 ^n+l,n n n n 

Next P* into P^+i is derived. Observe 

i A * V * 111 « V 

X j - X - <p X + A Y - (p , X 
n+1 n+1 T n+1, n n n n n+l,n n 


+ 

<Vl,nV 

Ci.n X n> 

= 

* 

^n+l, n 

(X* - 
n 

X ) 
n 

+ <p * . 

r n+l,n 

+ 

A* Y + 

• A u 




n n 

n n 


= 

* 

^n+l.n 

(X*- 

n 

X n> 

+ (<p* 

VY n+l,i 

+ 

A* M ) 

X +, 

A* u 


n n 7 

n 

n 

n 

n+1 

- X n+1 

= ^n+i.n 

(X* - X ) 
v n n 7 




1 u 
n n 


since (<P* . - cp , 4 ) = 0 

VY n+l,n r n+l,n 7 

P* =E [{<z>* * (X* - X ) + A* u } {<p *. 4 (X* 

n+1 llv n+l,n l n n 7 n n J LY n+l,n n 


- X ) + A* u T } ] 
n n n } 1 


6 



ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


* * * * T A * ■ A T 

P n+1 ~ <Vl,n P n<Vi,n + A n Q n A n 

P n+1 ~ ^n+l.n P n ^n+l,n + ^n+l,n ^n+l.n^ 


+ A*Q A* T 
n n n 


Finally 

p* = (p* P* w + 

n+1 n+l,n n r n+l,n 

\(o . P f (p . -ip . 1 + A Q A }. 

xv n+l,n n ls n+i,n r n+l,n J n n n J 

Now, (<Pn£,n -< Pn+l,n) = °- This is shown to be true 
by the following 

<p* . P * (<p* . - <p . ) T +A*Q a* t = 0 

r n+l,n n n+l,n r n+l,n n n n 

<p* . P* (-A*M ) T + A*Q A* T = 0 
r n+i,n n v n n' n n n 

by definition of n ’ 

[~(P* . P* m t + a* Q ] a* t = 0. 

1 r n+l,n n n n n J n 


Since 

* _ A * TVJT 

^n+i,n ^n+l,n n n 

(p A P* M T - A* M P*M T -A*Q =0 
r n+l,n nn nnnn nn 

w , P* M T - A* (M p* m t + Q ) = o 
n+l,n n n nnnn n 


and since 

A* = (p . P* M T (M P* M T +Q )' 
n r n+l,n n n n n n n' 

then 


n/rT 


<p A P* M. L - (p . P* M T = 0, 

^n+l,n n n n+l,n n n 


The final result is exhibited as follows. Given 
X* , p-%and tables for <p n+1>n , M n> and Q n> then 
. an onboard computational cycle is 


A* = <p P* M T (M P* M T +Q ) _1 

n ^n+l,n n n n n n n' 


* A * iur 

ip = ip - A M 

^n+l,n n+l,n n n 

P* = cp* 4 P * cp T 4 
n+1 ^n+l.n n r n+l,n 

X* = <p* . X* + A* Y* m 
n+1 n+l,n n n n 

where the measurements of Y* m are given each cycle. 
Some insight is given to these relations if Q n or P n 
equals zero. If P n =0, 

A* = 0 
n 

^n+l,n ^n+l,n 

Cl - 0 

Y 5 * 4 ,r\ 

n+1 ^n+l,n n ’ 

These equations imply that no error in from X^ is 
expected. No statistics are involved since only the 
transition matrix appears. 


“ V 0 ' 

A n'^l,„ M » 


-1 


^n+l.n ^n+l,n ^n+l.n 


" 9^.4 . = 0 


Ci * 0 


v* iur-1 „ v* m 

X = ip . M Y = ip , , x 

n+1 n+l,n n n n+l,n n 


These equations imply that the expected value of noise 
is zero, and Y* m is Y n . Thus X^ +1 = X n+1 . 


CONCLUSION 

Since instrumentation will always contain some 
errors, some method of accounting for these errors is 
necessary. The basic cycle of the Kalman filter as 
applied to navigation is estimated trajectory, meas- 
urement, and new estimated trajectory. 
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If the four equations of the final result are re- 
called, onboard computation and updating of the sta- 
tistics are seen to be possible. 

The basic weakness of this derivation is the as- 
sumption that the actual trajectory is sufficiently close 


to the reference for the state transitions and geometry 
mapping matrices to function accurately. This re- 
striction can be relaxed by obtaining operators (pos- 
sibly a matrix of partials) to transform <P n +i n and 
M n such that the transformed values will function away 
from the reference. 
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N65-23 7 83 

B. A MATHEMATICAL METHOD TO DETERMINE THE great circles as shown in Figure i. The axes passing 
ATTITUDE OF A SPHERE by Hans F. Kennel through the great circle intersections form the orthog- 

onal sphere coordinate system. These intersections 
will henceforth be called axis pierce points. Each 
great circle is split into half circles and each circle 


331 ® 3 ' 


A method for determining ‘the attitude of a sphere 
without restricting its angular freedom is presented. Only 
the mathematical aspect of the problem is treated . Neither 
implementation nor accuracy of performance is discussed. 
The method has an application in spacecraft attitude 
simulators supported on spherical gas bearings. 

The sphere is marked with three mutually perpen- 
dicular great circles and intersected in a minor circle by 
a stationary plane. The intersections of the great circles 
with the minor circle define the attitude of the sphere 
relative to the plane. 

The orthogonal transformation matrix is developed, 
which allows the calculation of the Euler angles of any | 
_ desired sequence . 



FIGURE 1. SPHERE MARKINGS 


INTRODUCTION 

A method is needed for determining the instantan- 
eous attitude of a sphere relative to a room-fixed 
reference coordinate system. The determination must 
be made without knowledge of the history of the sphered 
angular motion. The described method accomplishes 
this and imposes no angular restrictions on the sphere. 
One application is for a spacecraft attitude simulation 
table supported by a spherical gas bearing. There will 
beother applications; therefore, theproblemis treated 
in general. 

The proposed method works as follows. The 
sphere is marked with three mutually perpendicular 


is marked or given some distinguishing characteristic. 
The six distinctly identifiable half circles prevent am- 
biguity in the attitude determination. 

The half circles are applied to the sphere in such 
a way that each axis pierce point is formed by three 
half circles; one is continuous through the pierce point 
and the other two join at the pierce point. This is 
shown in Figure 2, which gives a direct view of any of 
the six axis pierce points of the sphere. This leaves 
two choices of orientation with respect to the sphere 
coordinate system; the end points of the half circle 
that passes through one axis may touch either of the 
remaining two axes. Once this choice is made, the 
orientation of all the half circles is fixed. 


9 
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A stationary plane intersects the sphere in a minor 
circle. This plane (read plane, Fig. i) must be lo- 
cated so that the radius r of the minor circle (read 
circle) is large enough to eliminate the possibility that 
in some cases it may be intersected by only one half 
circle, thus providing insufficient information for at- 
titude determination (Fig. 3). This possibility is 
eliminated if a radius r greater than RA j~2 is selected 
(where R is the radius of the sphere). Then, four or 
six intersections will occur, depending on the position 
of the sphere. v 



Angles are measured in the stationary read plane 
from a reference on the read circle to the intersections 
between the read circle and the half circles of the 
sphere. These angles, along with the identification of 
the half circles, uniquely define the attitude of the 
sphere with respect to the reference coordinate sys- 
tem. 

The introduction of an additional coordinate sys- 
tem allows vast simplification of the problem. This 
additional coordinate system is superimposed on the 
sphere coordinate system in such a way (one of 24 
possible ways)that its attitude can be determined with- 
out the aid of the half circle identifications. That is, 
for the attitude determination of the superimposed co- 
ordinate system with respect to the reference coordi- 
nate system, the sphere only needs to be marked with 
three great circles having no distinguishing charac- 
teristics. 

The problem may now be solved in three distinct 
steps. First, the orthogonal transformation matrix 
from the reference coordinate system to the super- v 
imposed coordinate system is found. Second, the 
superimposed coordinate system is correlated with 
the sphere coordinate system, resulting in the orthog-_ 
onal transformation matrix from the reference system 
to the sphere system. Third, the appropriate direction 
cosines of the orthogonal transformation matrix are 
used to calculate the selected rotation angles. The 
selection of the rotation angle sequence is arbitrary 
with respect to the method discussed here. 

The following mathematical treatment assumes an 
ideal case, i. e. , a perfect sphere, infinitely thin lines, 
perfect intersection angle information, etc. An error 
analysis for a practical case will be left to a future 
study, as will actual mechanization of the readout, 
identification of the intersection angles and half circles, 
and high speed digital data processing to allow readout 
with minimum time lag to accommodate the dynamics 
of a vehicle control loop. This loop is usually of low 
response for control investigations on gas bearing 
supported simulation tables. 

NOMENCLATURE 

a direction cosine 

h distance along Z axis 

i subscript denoting x, y, or z 
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J 


. L 


r 

r* 


R 

sgn L 
u, v, w 


x, y, z 


X, Y, Z 

X , Y , Z 
r r 


r 


a 

u 


a 

v 




<Pl> <^ 2 » etc * 


<f> > <P 

U 1 V U 2 


^Vl’ ^V2 


0 

A0 


<p 

v 


subscript denoting X, Y, or Z 

half circle identification number 
(label) 

read circle radius 


$ 

6 

* 


rotation about the z" axis 
rotation about the y f axis 
rotation about the x axis 


>(Fig. 12) 


radius of the circle touching the 
three axis pierce points of a sphere 
octant 

sphere radius 
A V | L | 

superimposed coordinate system 

sphere coordinate system 

reference coordinate system 

coordinate system parallel to X, 
Y, Z with the Z and Z r axes coin- 
ciding and the origin at Z = -h 

angle between the -u axis and the 
XY plane 

angle between the v axis and the 
XY plane (Fig. 4) 

angle from X r , positive counter- 
clockwise 

angles of the half circle intersec- 
tions with the read circle 


DEFINITIONS 

Let x, y, z be the sphere coordinate system and 
X, Y, Z the room-fixed reference coordinate system. 
Figure 1 shows both systems coinciding. Another 
system, X r , Y r , Z r , is parallel to X, Y, Z with the 
Z r and the Z axes coinciding and the origin of X r , Y r , 
Z r the distance h below the origin of the X, Y, Z sys- 
tem. The read circle is defined by the intersection of 
the X r , Y r plane (i.e., the read plane, Fig. 1) with 
the sphere. The intersection angles are read positively 
from X r toward Y r . 


The axes of the superimposed coordinate system 
are labeled u, v, w. For zero rotation of the super- 
imposed system, u coincides with X, v with Y, and w 
with Z. The angle between u and X, as well as between 
v and Y and w and Z, will not exceed ± 7r/2. 

For identification, the sphere half axes are as- 
signed the following numbers: 

+x = +1 +y = +2 +z = +3 

-x = -i -y = -2 -z = -3. 


The half circles are also identified by numbers. 
Note that the magnitude of the half circle numbers is 
the same as the magnitude of the numbers for the axes 
perpendicular to the half circles. The following identi- 


pair of intersection angles caused 

fications were chosen (Fig. 

1): 

by half circle(s) perpendicular to 



the v axis 

Axis pierce point 



through which 

Half circle 

A 0,5 (d) + d> ) 

= vy u 2 y ur 

half circle passes 

identification 

A 0. 5 (d> - <b ) 

= U 2 Mil 

+x 

-x 

+2 

-2 

(<P v 2 + <t> vl ) 

+y 

-y 

+3 

-3 

^°- 5 <* V |- V 

+z 

-z 

+1 

-1 


pair of intersection angles caused 
by half circle (s) perpendicular to 
the u axis 
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ATTITUDE DETERMINATION OF THE 
SUPERIMPOSED COORDINATE SYSTEM 


The intersections of the great circles with the 
read circle (there may be four or six intersections, 
depending on the position of the sphere) are labeled 

0i, 0 2> etc. , in the order in which they appear when 
the angle 0 is increased from its reference X^.. 

Any two intersections formed by the same great 
circle are considered a pair. Each pair will allow the 
calculation of the three direction cosines associated 
with the axis perpendicular to the great circle causing 
the intersections. 

The first intersection pair defines the.v axis and 
is called 0 V1 and 0 V2 . The second pair (0 ul and 0 U2 ) 
defines the u axis. If there are six intersections, two 
are disregarded as explained in the next section. 

The following quantities are calculated for the v 

axis: 

A 0 A 0.5(d) - 0 ) 0 A 0.5(0 +0 ) 

= V2 VI V= V2 VI 

cos A0 

V 

2 + cos 2 A0 

V 

, (h/r) . ( 1 ) 

V + n/ (h/r) 2 + cos 2 


sin a = 


\l (h/r) 


z 



A<Ml<*y.-<*>v,) 


FIGURE 4. 


♦4(F) 


cos A<£ v 
'+ cos 2 A<£ y 

(fr) 

+4(^) 2 +cos 2 A^ y ' 


ANGLE CALCULATIONS 


cos a = 


(h/r) 


U + nT (h/r) 2 + cos 2 A0 


( 2 ) 


The orthogonal transformation matrix may be written 
now. 


This is illustrated in Figure 4. 
Similarly for the u axis: 


u = 


sin a = 


^U2 ^ul' 
cos A 0 


u = 



V 


1 

p 

e 

X 

a v 
uY 

a uZ 


X 

5 (<f> + <6 ) 

r u2 r Ul 

V 

1 


a 

vX 

a __ 
vY 

a vZ 


Y 


w 


a 

a 

a „ 


Z 


J 


wX 

wY 

wZ 


. 


U + \/ (h/r) 2 + cos 2 A 0 


with the values shown below. 


(3) 


a , = -cos a cos 0 
uX u ^u 


a = -cos a sin 0 
uY u u 


a „ = -sin a 
uZ u 


a v = cos a 
vX v 


a = cos of 
vY v 


a ^ = sin a 
vZ v 


cos 0 

v 

sin 0 

v 


a _ a _ r a a a 

wX uY vZ uZ vY 

— a ~ a „a — — a a 
wY uX vZ uZ vX 

a ^ = a __ - a a 
wZ uX vY uY vX 
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The difference in sign between the direction 
cosines for u and v is explained by the fact that the 
angle 0 V is the angle of the projection of the +v axis, 
on the read plane, whereas <p u is the angle of the -u 
axis projection. 

As mentioned earlier, the orthogonal transforma- 
tion matrix from the reference coordinate system to 
the superimposed coordinate system is found without 
recourse to any great circle split or identification 
other than proper pairing. 

DETERMINATION OF THE COORDINATE SYSTEM 
CORRELATION 

The superimposed coordinate system may have 
any one of 24 different superpositions with respect to 
the sphere coordinate system. The proper superposi- 
tion is identified by combining two separate pieces of 
information. First, the absolute value of the half 
circle identification number indicates which axis is 
correlated. Second, the sign sequence of the half 
circle intersections indicates the sequence of the signs 
for the axes coinciding with u, v, w, respectively. 

There are eight possible sign sequences for the 
four intersection case (Fig. 5). This can be seen in 



FIGURE 5. FOUR INTERSECTION CASE 




+y +y 



FIGURE 6. POSSIBLE FOUR INTERSECTION CASES 


Figure 6 by starting the sequence clockwise in each 
quadrant of the six circles shown. Each sign sequence 
appears three times with different number magnitudes. , 
One of the three possibilities for the sequence +, 

+ , - is used as an example. The eight sign sequences 
are shown in Table 1. For each of these sequences, 
there is a sign sequence for the axes coinciding with 
u, v, w. Labels L^, Ly 2 , etc. , are the half circle 
identifications. 

For the case of six half circle intersections, there 
are several subcases: zero, one, two, or three axis 
pierce points may be visible within the read circle 
(Figs. 7,8,9, and 10). These subcases will be called 
6 0 , 6i» 6 2 , and 6 3 , respectively. Cases 6 0 and 6 3 are 
mutually exclusive, depending on whether the read 
circle radius is smaller or larger than r* = R sf 2/3, 
where r T is the radius of a circle touching three axis 
pierce points belonging to the same sphere octant. 
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FIGURE 7. SIX INTERSECTION CASE (NO AXIS PIERCE POINT IN FIGURE 9. SIX INTERSECTION CASE (TWO AXIS PIERCE POINTS IN 

READ CIRCLE) READ CIRCLE) 



FIGURE 8. SIX INTERSECTION CASE (ONE AXIS PIERCE POINT IN FIGURE 10. SIX INTERSECTION CASE (THREE AXIS PIERCE 

READ CIRCLE) POINTS IN READ CIRCLE) 


14 











ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


TABLE 1 . FOUR INTERSECTION CASE 


TABLE 2 . SIX INTERSECTION SUBCASE 6 0 


sgn L y 

sgnL u 

sgn L 

w 

V 

u 

w 


+ 1 

+ 1 

+1 

- 

+ 

- 

+i 

+ 1 

-1 

- 

- 

- 

+1 

-1 

+1 

+ 

+ 

- 

-1 

+ 1 

+1 

- 

+ 

4 - 

-1 

-i 

-1 

+ 

- 

+ 

-1 

"I 

+1 

- 

- 

+ 

-1 

+ 1 

-1 

+ 

- 

- 

+i 

-1 

-1 

+ 

+ 

+ 



To determine the subcase by inspection of the half 
circle intersections, it should be remembered that 
they are labeled <f > l9 etc, as they appear . A 
pair (belonging to the same great circle) can be un- 
separated (</> 1 0 2 > 0203» ••• > 0 G 0i), or separated by one 

(0103> 0204> • • • 9 0602) or 1^° (0104> 0205> • ■ • » 0603) 
intersections belonging to other pairs. These proper- 
ties are used to determine the proper subcase. 

Subcase 6 0 has three unseparated pairs (Fig. 7 ). 

Subcase G^as one unseparated pair and two single- 
separation pairs (Fig. 8). 


Subcase 6 0 must be treated separately; Table 2 applies 
with Lyj = Ly 2 = Ly, etc. This is also illustrated in 
Figure 11 . 

Using the correlation between the superimposed 
coordinate system and the sphere coordinate system, 
the reference-system/ sphere-system orthogonal 
transformation matrix may be established. This is 
done by writing a correlation matrix and multiplying 
it with the reference-system/ super imposed- system 
matrix. Consider, for example, the following sequence 
of identification numbers: 

Lj = +3; L 2 = +2; L 3 = -3; L 4 = +2; 


Subcase 6 2 has two single -separation pairs and 

one double -separation pair (Fig. 9). or after proper pairing. 


Subcase 63 has three double- separation pairs (Fig. 

10 ). 

For the matrix calculation, subcases 6*, 6 2 , and 
63 are reduced to the four intersection case by disre- 
garding two intersections. In 64 the unseparated pair 
is disregarded, in 6 2 the pair with the smallest A0 is 
uis regarded, and in 63 the third pair is disregarded. 


L - +3; 

Vi 


L 

V2 


- 3 ; L = +2; L = + 2 . 
ul u2 


The second sign sequence in Table 1 applies, with the 
following result: 


v = -z u = +y w = -x 


and accordingly 
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+ Z +Z 



BOTTOM VIEW 

FIGURE 11. POSSIBLE CONFIGURATIONS FOR CASE 6 Q 


X 


1 

o 

o 

! 

M- , 


u 

y 


+10 0 


V 

z 


0-1 0 


w 


X 


0 

0 

-1 

y 


+ 1 

0 

0 

z 


0 

-1 

0 


X 


a wX 

-a __ 
wY 

y 

= 

a uX 

a V 

uY 

z 


~ a vX 

" a V 

vY 


a uX 

a T7 

uY 


a uZ 

a vX 

a V 

vY 


a vZ 

a \vX 

a 

wY 

E wZ 

a wZ 


X 


a ^ 
uZ 


Y 


-a ry 

vZ 


Z 



DETERMINATION OF EULER ANGLES 


In general, the orthogonal transformation matrix 
has the form: 


X 


a xX 

a v 

xY 

a xZ 

y 


a V 

yx 

a v 
yY 

a yz 

z 


_ a zx 

a zY 

a zZ_ 


X 

Y 

Z 


(4) 


Any order of rotation may be treated. The sphere 
axes are designated x n , y", z" (.identical with X, Y, Z) 
before rotation. The first rotation through the angle 
<i>, for example, is about the z" axis, resulting in x T , 
y f , z", The second rotation through the angle 0 is 
about the axisy f , resulting in x, y ! , z ! . The third 
rotation through the angle * is about the x axis, re- 
sulting in the final attitude x, y, z. This is illustrated 
in Figure 12 and expressed in matrix notation in equa- 
tion 5. . Multiplication of the individual rotation ma- 
trices gives the orthogonal transformation matrix 
expressed in equation 6. 


X 


1 

0 

0 


cosO 

0 

-sin0 


cos<£ 

sin$ 

0 


X 

y 

- 

0 

cos^ 

simJ' 


0 

1 

0 


-sin<£ 

cos$ 

0 


Y 

z 


0 

-sin ^ 

cos^ 


si nO 

0 

cos 0 


L 0 

0 

1 


Z 
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X 


y 

z 


COS # COS0 

sin# cos0 

-sin0 

cos# sin0 sin# 

sin# sin0 sin# 

cos0 sin# 

-sin# cos# 

+cos# cos# 

cos# sin0 cos# 

sin# sin0 cos# 

COS0 cos# 

+sin# sin# 

-cos# sin# 


( 6 ) 


z" 



X 


ORDER OF ROTATION 

1. <t> CCW ABOUT z"- AXIS 

2. 9 CCW ABOUT y'-AXIS 

3. ¥ CCW ABOUT x-AXIS 

FIGURE 12. ROTATION ANGLES 

Comparison of the direction cosines between 
equations 4 and 6 leads to 

4> = arc tan (a„/a„) 

A X XA 

0 = arc sin (-a xZ ) 

* = arc tan (a yZ /a zZ ). 


To make the solution single valued, the following 
limits must be imposed: 

-7T < +7T 

-7r/2 s 9 =£ +7i/2 

-7T < #< +7 T. 

These limits force cos 0 ^ 0 allowing the determi- 
nation of the. proper quadrant for # and # by 

(sgn a xY )/(sgn a^) 

and 

(sgna yZ )/(sgn a^). 

One sample was calculated and is presented in 
the appendix. 


CONCLUSIONS 

It appears that the mathematical solution presen- 
ted can be implemented with reasonable cost while 
rendering adequate accuracy for the problems to be 
investigated on gas bearing supported simulation 
tables. An exact error analysis, however, to deter- 
mine the accuracy of the attitude readout is left to a 
future study. 

The orthogonal transformation matrices and the 
Euler angles can be calculated with standard digital 
computer equipment. 

A method for the readout of the intersection angles 
must be devised and the necessary equipment developed 
and mechanized. 
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APPENDIX 

DERIVATIONS OF CORRELATION FORMULAS 

A. FOUR INTERSECTION CASE 

The correlation information contained in Table 
1 for the four intersection case may be condensed into 
a mathematical formulation. This is demonstrated 
for the axis coinciding with v. Two modifiers are in- 
troduced: 0.5(1 + sgn Ly^ and 0. 5 ( 1 - sgn L^) . 
Each can only assume the value zero or +1. One is 
present when the other vanishes and vice versa, thus 
affording a selection. For sgn h Vi = + 1, we can see 
from Table 1 that sgn v = sgn L U2 ; and for sgn h vl = - 1, 
sgn v = sgn Lmj. The combination of the above state- 
ments with the modifiers and the absolute value of the 
label results in 

v = -Lvi [0.5 (1 + sgn hyi) sgn L^ 

+ 0. 5 (1 - sgn Lyj) sgn L Ul ] . 

Similarly the correlation for u is 
u = L U1 [ 0. 5 ( 1 + sgn L^) sgn L^ 

+ 0. 5 (1 - sgn \ x ) sgn L^] . 

Axis w is correlated by the formula: 

W = -0. 25 1 12 - jLyj | - | Lu 2 | - | L V1 1 - l^l ] L 

with Xsgn L = sgn L U1 + sgn L^ + sgn L vl + sgn Ly 2 f 

which can be explained as follows. The sum of the 
absolute values of all half circles is 12. For the four 
intersection case under discussion here, there are 
never more than four intersections showing. There- 
fore, subtracting the magnitude of the visible labels 
from 12 leaves twice the magnitude of the not-visible 
labels, which provides the magnitude of the axis with 
which w correlates. Summing the label signs deter- 
mines whether w coincides with a plus or minus axis 
(compare Fig. 6 ). 

B. SIX INTERSECTION CASE, SUBCASE 6 0 

When six intersections between the half circles 
and the read circle occur and no axis pierce point is 
visible within the read circle (Fig. 7), the following 
formulas apply. These formulas were derived from 


Table 2, which summarizes the eight possible sign se- 
quences. When sgn is positive, multiplication by 
sgn I^y gives the proper sign for u. When sgn L u is 
negative, it must be multiplied by sgn Ly and the re- 
sulting sign inverted. This is achieved by the formu- 
la: 

U = ‘ L u [ °‘ 5 ( 1 " Sgn V Sgn L v 
-0.5 (1 + sgn LJ sgn LJ . 

Again two selective modifiers are used to cause the 
unwanted part of the equation to drop out. The equa- 
tions for v and w are similar: 

v = L [0.5(1- sgn L ) sgn L 
v 1 v w 

- 0.5 (i + sgn L y ) sgn LJ 

w = L w [°* 5 (i - sgn L w ) sgn L u 

- 0.5 (1 + sgn L w ) sgn LJ . 

SAMPLE CALCULATION 

Assume four intersections <j> i through 0 4 identified 
by the angles from the reference X r and by the labels 
of the intersecting half circles. 

<Pl = 45. 4° 02=61.0° 03=165.5° 0 4 = 301. 1° 

Lj = -1 L>2 - +3 L 3 = +1 L 4 = +3 

To find the direction cosines of the axes perpen- 
dicular to the intersecting great circles, the angles 
0 V , A 0 V , 0 U , and A0 U must be identified. 

0 v = 0. 5 (165. 5 + 45. 4) = 105. 45° 

A 0 v = 0. 5 (165. 5 - 45. 4) = 60. 05° 

0 u = 0. 5 (301. 1 + 61. 0) = 181. 05° 

A 0 u = 0. 5 (301. 1 - 61. 0) - 120. 05° . 

Then 

sin 0 = + 0 . 96386 
v 

cos 0 = -0. 26640 
v 
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sin Ad> = +0. 86646 
v 

cos Ad> = +0. 49924 
v 


The correlation between the superimposed and the 
sphere coordinate system is found from Table 1 , fourth 
sequence, to be: 


sin d> = -0. 01832 
r u 

cos <b = -0. 99983 
^u 

sin Ad) = +0. 86559 
r u 

cos Ad) = -0. 50076 . 
u 

Substituting into equations 1 and 2 results in (with h/r = 1) 

sin a = +0. 44667 since =-0.44776 

v u 

cos a = +0. 89470 cos a = +0. 89416 


v 


u 


u = -z v = +x w = -y. 

With this correlation, the above matrix can be trans- 
formed into the final orthogonal transformation matrix: 


X 


-0.23835 +0.86237 +0.44667 


X 

y 

= 

+0.37882 -0.50605 -0.77487 


Y 

z 


-0.89401 -0.01638 -0.44776 


Z 


This results in 


and allows the orthogonal transformation matrix of 
equation 3 to be written. Because the intersection 
angle values were derived graphically, the matrix is 
not precisely orthogonal. 


u 


+0. 89401 +0. 01638 +0. 44776 


"x" 

V 

= 

-0. 23835 +0. 86237 +0. 44667 


Y 

w 


-0. 37882 +0. 50605 +0. 77487 


Z 


tan $ = +0. 86237/-0. 23835 = (+1/-1) • 3. 6181 
sin 0 = -0. 44667 

tan * = -0. 77487/-0. 44776 = (-1/-1) • 1. 7305. 

The sign relationship of the direction cosines indicates 
that the angle 4> is in the second quadrant and that the 
angle ^ is in the third quadrant. The final result is 
$ = +105° 37’ , 6 = -26°32\ and * = -120° 1\ 
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C. DEVELOPMENT OF SQUEEZE FILM BEARINGS by 
Peter H. Broussard, Jr., John Burch, and Emmett L. Martz 


N65-23784 


% 


flirt 


In the past, the term "gas bearing" has referred to 
bearings In which the load-carrying member is supported 
on a gas film resulting either from external pressurization 
or from the high-speed relative lateral motion of the 
bearing surfaces. In the bearing described in this arti- 
cle, the necessary gas pressure is obtained by vibrating 
one bearing surface in a normal direction at a high fre- 
quency, thus alternately squeezing and releasing the gas 
film between the surfaces. This bearing does not require 
either external gas supply or high-speed rotation. 


The first part of the article describes the theoretical 
development of the squeeze bearing using Reynolds’ equa- 
tion. The second part describes the design of a working 
model in which a piezoelectric sleeve and end plates 
provide motion normal to the supported member at all 
points. 




INTRODUCTION 


The most commonly used gas-lubricated bearings 
have been either hydrostatic or hydrodynamic. The 
theoretical solution of either follows from Reynolds’ 
equation. For the hydrostatic type, no relative motion 
of the bearing surfaces is necessary, as a pressure 
gradient is maintained by external pressurization. For 
the hydrodynamic type, a pressure gradient is main- 
tained by lateral, relative motion of the bearing sur- 
faces. Also contained in Reynolds’ equation is a term 
pertaining to the rate of change of the product of gap 
height and density with respect to time; this is the 
type of motion treated in this article. It is shown that 
a net load carrying capacity is generated by cyclic 
squeeze motion of two parallel surfaces. 


squeeze film bearings in the gyro field occurred to 
Dr. B. Sternlicht of Mechanical Technology, Inc. 
Astrionics Laboratory commenced work in squeeze 
films in November 1963, after having become aware 
of the work reported in Reference 1. 


NOMENCLATURE 


R i 

r 2 

r 

e 

^o 

h 

Ah 

P 

P 

M 

t 

co 

T 

R 

c 

Pa 


inner radii of annuli 

outer radii of annuli 

radial coordinate 

angular coordinate 

mean gap height 

instantaneous gap height 

peak amplitude of excursion 

pressure 

density 

viscosity 

time 

frequency 

absolute temperature 
gas constant 
ambient pressure 


ANALYSIS 

Consider two parallel flatannuli as shown in Fig- 
ure 1 > the top annulus to be fixed and the bottom an- 
nulus to be capable of executing vertical cyclic mo- 
tion. Let the inner and outer radii of the annuli be R t 
and R 2 , respectively. Denote the instantaneous gap 
height between the annuli by hand let the bottom annu- 
lus execute a sinusoidal motion (Fig. 2) such that the 
instantaneous gap height is given by 


According to Reference 1, the earliest mention of h = h - Ah cos cot 
gas squeeze films seems to have been made in 1954. 

At least as early as 1960, the idea of utilizing gas or, if made dimensionless by division of h Q , 
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FIGURE 1. PARALLEL FLAT ANNULI 


P = £ — , R = ~ C - — , H = , t = cut, cu = frequency 

p rv>— n 

a z 1 o 

where = ambient pressure 

and by assuming that the process is isothermal: 

p 

p = 9 f* c = gas constant, T = absolute temperature 

c 

Substituting into equation 1 and noting that 



FIGURE 2. BOTTOM ANNULUS IN SINUSOIDAL MOTION 


H = — =1-6 cos cut 
h 

O 


a _ i a , a a 

9r R 2 - Rj 9R &nd 9t _ W 9 t 


_L 



Ah 

equation 1 becomes 

~T 

1 <L 

R 8R 

r rh 3 p n 

|_ 9R 




9(PH) 


9t 


( 2 ) 


where cr = squeeze number = 


12ptw(R 2 - R t ) 

p h 2 
a o 


where e = ~r — 

h 

o 


h^ = mean gap height. 


The dynamic behavior of the gas is described by 
Reynolds' equation in cylindrical coordinates: 


i_ d_ h^ ap _ a(ph) 

r 8r r 12 p P 3r at 


where r 


radial coordinate 


h = gap height 


p = pressure 


( 1 ) 


For a specified motion of the bottom plate, H = H(r) , 
the pressure P can be foundas a function of R and t in 
theory. Practically, however, equation 2 cannot be 
solved in closed form and recourse must be made to 
solution either by numerical methods or by perturba- 
tion techniques, either one of which is relatively leng- 
thy. 


Fortunately, an approximate method exists (Refs. 
1 and 2) that has given excellent agreement with ex- 
periment where the squeeze number a approaches in- 
finity (100 or more) . As can be seen from the defini- 
tion of a, this can happen in several different ways, 
but the cases in which either co °° or h — 0 are ap- 
propriate here. 

Let both sides of equation 2 be integrated with re- 
spect to t over one cycle, 27 r: 


p = density 
/ x = viscosity 



3R 


rh 3 p 


3P 

m 


dr 


r 2w a/PHl 

a j ^_^- dT 

o 


3 T 


t = time. 


0 , 


This equation can be made dimensionless by letting 


since PH is periodic of period 27r. Thus 


(3) 
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or 


r 2tt 9_ 

J 8R 

o 


- ap 

L RH P 9rJ 


dr = 0. 


Interchanging the order of integration and differentia- 
tion, 


Is / 2 '[™’ p If] * - 0 


and integrating 


j 2 ’ [ Ri,,p i] * ■ c 


or, since R is not a function of r, 


f [ H3p £] 


a T . - . 


(4) 


Since H is not a function of R, the integrand can be 
written as 


5 Is < Hlp2) • 


equation 4 can be written as 

r [« Hip! ] dT ■ f ’ 


and, interchanging the order of integration and differ- 
entiation, 


_8_ 

8R 


C M 


dr = 


R 


Integrating, 


C M 


dT = C< In R + C<> 


(5) 


R 

Applying the boundary conditions P - 1 at R = 1 — 

R2-R1 


and P = 1 at R - 


R, 


R 2 -Rl 


gives two equations: 


J 27r H 3 dT = Cj In 


R 2 - Rj 


/ 2 * hj dr - c, 11 


+ Co 


+ Co . 


These two equations imply that C t = 0. Thus 


C 2 = J 2n H 3 dT = J Z7r ( 1 - e cos t ) 3 dr 


27 r 


o o 

or C 2 = 2 tt ( 1 + 3/2 e 2 ) . 
Therefore, from equation 5, 

J 2ir H 3 P 2 dr = 2n ( 1 + 3/2 e 2 ) , 


(6) 


Equation 6 is used now for the* case in which the 
squeeze number cr is very large, i.e. , a-* From 
equation 2, if a— * <», then 

MPH) = 

8 r 


or, integrating, 
PH =f(r). 


(7) 


where f(r) is to be determined. The fact that PH is a 
function of r is explained as follows: H can be a func- 
tion only of t, since the plates are parallel. P is a 
function of H, t, and r. P is not a function of angular 
position on the annulus 0 , since there is no circum- 
ferential variation in pressure. Therefore, P can 
depend only on H, t, and r. This means that the pro- 
duct PH can be a function of, at most, two variables, 
r and t . If 

a (PH) 
dr 

then it follows that PH is a function of, at most, one 
variable, r. Substituting equation 7 into equation 6 
gives 


f 2ir H 3 P 2 dr= J H( PH) 


27T 


= f * f (r) 2 Hd r = f(r) 2 J Zir (1 - e cos r ) dr 
o o 

= 27r f(r) 2 . 


dr 


2?r 
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That is, 

f 2n H 3 P 2 dr = 2 tt f(r) 2 . (8) 

o 

Equations 6 and 8 together imply that 


p a (1 + 3/2 c 2 ) ' /z 

p - 

K I - € COS a>t 
€ = 0.5 


27r f (r) 2 = 27T ( 1 + 3/2 e 2 ) 
or 

i 

f(r) = (1 + 3/2 e 2 ) 2 " . 

Using equation 9 in equation 7 gives 

1 

p = UlL = (1 + 3/2 € 2 )2 

H i - e cos cut 


(9) 


( 10 ) 


or in dimensional form, 
JL 

p (1 + 3/2 e 2 ) 2 
a 

^ 1-6 COS CUt 


(ID 


Equation li shows that p is a function of time but not 
of R. '' The plot of p versus time is shown in Figure 
3 for e = 0. 5. 

As can be seen from the plot, the average pressure in. 
the gap is higher than ambient pressure. This means 
that there is a net force. The force F caused by the 
pressure in the gap is given by 


F = f Ri 27rrpdr = 27rp ( 1 + 3 / 2 g2 > 
J a i - c cos cot 


jR 2 


rdr 


or 


F = 


7T (R 2 2 - Rj 2 ) 


(R/-R! 2 ) p a (i + 3/^€ 2 ) 2 


1 - € COS cot 


( 12 ) 


* It has been pointed out that this pressure distribution 
fails to meet the boundary conditions at the edges, 


P = P ( 


am 



FIGURE 3. VARIATION OF PRESSURE WITH TIME, £=0.5 


The average force over one cycle is given by 

^ = 2tt f 2IF Fdr = n ^ " Rl ^ P a < 1 +3 A 
o 


1 r 2ir 

2 7T J 

O 


d t 

1 - € COS T 


or 

F = 7T (Rj 2 - R t 2 ) P a 


1 + 3/2 € 2 \ h 
l-c* ) 


(13) 


since 



dr 

1 - € cos t 


27 r 

1 

( 1 - € 2 ) 2 


Equation 13 shows that the average force F is pro- 
portional to the ambient pressure, p . Hence, for a 
fixed e, the average force can be increased by raising 
the arnbi ent press ur e in wh ich the squeeze bear ing op - 
erates. Figure 4 shows a load displacement curve which 
is applicable to either a disc or an annulus of the 
same area. 
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F 

^ Ro Pa 


FIGURE 4. AVERAGE LOAD CARRYING CAPACITY VS 
ECCENTRICITY RATIO 


JOURNAL SQUEEZE FILM BEARING 



The expression for the force developed by the 
squeeze film can be used to determine the load carry- 
ing capacity of the journal bearing when certain as- 
sumptions are made. 

In Figure 5, the inner cylinder is subjected to a 
constant radial force W. If a squeeze film force of 
sufficient magnitude is being generated, the inner cyl- 
inder will be displaced an average distance d from its 
concentric position relative to the sleeve and will float 
in this position. The system is in equilibrium; that is, 
the squeeze film force is exactly equal to the applied 
force W. 

Assume that instead of being circular, the inner 
cylinder and journal are right cylindrical polygons 
having N sides each. This assumption will not intro- 
duce any appreciable error provided N is sufficiently 
large. Also, assume that the sleeve is subjected to a 
sinusoidal vibration and that the vibrating frequency 
is sufficiently large to give a squeeze number of 100 
or more. There are two possible means of vibrating 
the sleeve relative to the inner cylinder; each will be 
treated separately. 

i. The case of the sleeve vibrating in the direc- 
tion of the applied force W will* be considered first. 
The sleeve does not change size; that is, it does not 
expand and contract. The squeeze film force is gen- 
erated over N separate pads. Now, the component of 


FIGURE 5. JOURNAL BEARING DISPLACEMENT GEOMETRY 


the squeeze film force in the direction of the applied 
force W that is generated over the i^ 1 pad will be de- 
termined. 

th 

The gap height at the i pad is given by 


h = h. (1 
J 

where h 


h. 

] 


- € COS ) 

th 

= the gap height at the i pad 
= the radial clearance = R 2 - Rj 
= the eccentricity of the inner cylinder = 


6 , 

l 


d 

R 2 - R, 

th 

the angular position of the i gap. 


Since the gap height atany pad is considered con- 
stant across the pad, the expression developed for the 
thrust bearing to the journal bearing can be applied, 
remembering that only the component of the squeeze 
film force in the direction of the applied force W is of 
interest. 


th 

Therefore, the force generated over the i pad in 
the direction of the applied force W is given by 
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“♦f 


W. = P A< 
1 am 


Ah cos 0. 


h.( 1 - e cos 0 . ) 

. y }_ 


i - 


’ Ah cos 0. 


cos 0 . 

i 


h.( 1 - € cos 0.) 
J 1 


where W. = the squeeze film force generated over 
the i pad in the direction of the ap- 
plied force W 


P 

am 

A 


27tR 1 L 


ambient pressure 

the area of the i^ pad = 

L = the length of the inner cylinder 


N 


where 


Ah = the amplitude of vibration. 


Therefore, the total force generated over the en- 
tire sleeve in the direction of the applied force W is 
given in equation 14. 

2. The second means of vibration is given by the 
sleeve expanding and contracting sinusoidally with re- 
spect to time. This mode of vibration would be given 
by using a piezoelectric sleeve. In this mode of vibra- 
tion, for all practical purposes, the sleeve is moving 
normal to the inner cylinder at all points around the 
inner cylinder. For this mode of vibration, the total 
load carrying capacity of the journal bearing is given 
in equation 15. 

Computer runs were made to determine the load 
carrying capacity of the journal bearing using both 
modes of vibration. Separate runs were made for both 
the AB-5 and AB-3 bearings to determine the effect of 
varying bearing clearance sand amplitude of vibration. 
The results are plotted in Figures 6, 7, 8, and 9, 


W - 


2ttR x LP 
N 


am 


i=N 

i 

i=l 


2h. 2 (i - c cos 0.) 2 + 3Ah 2 cos 2 0.) 
J l }' 

2^ 2 ( i - e cos 0.) 2 - Ah 2 cos 2 0^j ^ 


cos 0. 

l 


(14) 



RADIAL LOAD (NEWTONS) 


FIGURE 6 . RADIAL ECCENTRICITY VS LOAD, AB-3 JOURNAL SQUEEZE FILM BEARING (PIEZOELECTRIC SLEEVE) 
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RADIAL CLEARANCE 

1.9 x ICf 3 cm 
1.5 x I0 -3 cm 
1.3 x ICT 3 cm 
1.0 x I0 -3 cm 


AMBIENT PRESSURE = 1.01 x I0 5 N m ~ 2 ABSOLUTE 


Ah = 5.6 x I0' 5 cm 


FIGURE 7. 


RADIAL LOAD (NEWTONS) 

RADIAL ECCENTRICITY VS LOAD, AB-3 JOURNAL SQUEEZE FILM BEARING (UNIDIRECTIONAL SLEEVE MOVEMENT 
AS WOULD BE OBTAINED ON A SHAKE TABLE) 


RADIAL CLEARANCE 
1.9 x I0 " 3 cm 
1.5 x I0 -3 cm 
■ 1.3 x I0 -3 cm 
- 1.0 x I0 -3 cm 

AMBIENT PRESSURE = 1.01 x I0 5 N m " 2 ABSOLUTE 
^ \ 7 Ah = 5.6 x I0 -5 cm 


RADIAL LOAD (NEWTONS) 

FIGURE 8. RADIAL ECCENTRICITY VS LOAD, AB-5 JOURNAL SQUEEZE FILM BEARING (UNIDIRECTIONAL SLEEVE MOVEMENT 
AS WOULD BE OBTAINED ON A SHAKE TABLE) 
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FIGURE 9. RADIAL ECCENTRICITY VS. LOAD, AB-5 JOURNAL SQUEEZE FILM BEARING (PIEZOELECTRIC SLEEVE) 


• A standard AB-3 assembly was vibrated on the 
shake table in an effort to compare theory with experi- 
ment for this mode of vibration (no piezoelectric 
sleeves were available at the time) . Several factors 
such as resonance effects, cocking of the inner cylin- 
der, inadequate means of loading, etc,, made it im- 
possible to accurately check theory against experiment 
for this mode of vibration. It was found that a suffi- 
cient squeeze film force was generated to float the in- 
ner cylinder for an applied radial force equivalent to 
1 g acceleration. 


EXPERIMENTAL 


A. Squeeze Film Generator. After checking the 
theory of the squeeze film bearing and proving that it 
works by tests on the shake table, the next step has 
been the development of self-energized experimental 
models. Three types of power were considered: 

1. Voice coil 

2. Magnetostrictive 

3. Piezoelectric 


Requirements of minimum size, weight, power, 
magnetic field, and self-heating have made the piezo- 
electric ceramic transducer our choice for this appli- 
cation. 

The ceramic is composed of a multitude of cry- 
stals, originally in random orientation. Its properties 
are the sum of the properties of all these crystallites. 
The ceramic is then poled by a strong dc field, where- 
by a majority of dipoles of the component crystals 
are aligned. After poling, if any electric field opposite 
in sense to the aligned dipoles is applied at higher and 
higher voltage, the dipoles of some crystals reverse 
their direction. To do this, the ions that form the di- 
pole suddenly shift position slightly and a net expansion 
or contraction of the ceramic takes place. A rapidly 
changing electric field causes an in-phase changing of 
physical dimensions of the ceramic. Conversely, a 
rapid squeezing of the ceramic produces an in-phase 
changing in the polarity of the ceramic. This motor- 
generator property of the poled ceramic is known as 
piezoelectricity. 

B. There are several problems associated with 
the use of piezoelectric ceramic. Some of these are: 
i. Elliptical Radial Airgapu Tests conduct- 
ed at MSFC with capacitive and photonic probes have 
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shown that when a piezoelectric ceramic sleeve is en- 
ergized by an alternating current across curvilinear 
electrodes the radial excursion is 2. 5 times as much 
at the center as near the ends. This asymmetric mo- 
tion generates an elliptical radial airgap. 

2. Resonant Frequency. By exciting the ce- 
ramic at its natural mechanical resonant frequency, 
the excursion can be amplified as much as 500 times. 
Resonance of a ceramic sleeve of equal length and 
diameter is about 1.6 times higher infrequency in the 
axial than in the radial direction. Although the reso- 
nant frequencies can be changed by changing the di- 
mensions, simultaneous axial and radial resonance 
cannot be obtained. 

3. Power and Heating. Power i§ limited by 
the physical strength of the ceramic and by the limit 
to which the material may be electrically driven be- 
fore de-poling begins. In a squeeze film bearing, the 
power may also be limited by internal heating, which 
is a partial function of frequency. 

4. Aging. The parameters of ferroelectric 
ceramics change with time after major disturbances 
such as temporary heating to high temperatures for 
electroding or the poling process. A reverse trend 
sets in after poling and is usually described as "aging 
per decade." The aging process is a logarithmic func- 
tionof time but can be speeded up by exposure to radi- 
ation. 

5. Electrodes. Electrodes should be good 
conductors, should not peel, and should be resistant 
to wear and corrosion. Electrolysis nickel plating, 
50 to 75 p thick, meets these requirements if the ce- 
ramic surface is carefully prepared in advance. 

6. Attachment of Endplates. Endplates can 
be mounted either rigidly or flexibly to the ceramic 
sleeve. Rigidly mounted endplates have the same ex- 
cursion and frequency as the sleeve-plus-endplate 
mass combination. This resonant frequency is slightly 
less than for the sleeve alone. Flexibly mounted end- 
platescan be "tuned" to have an excursion much larger 
than the sleeve and to resonate at a lower frequency. 

C. Conventional TestModels. There are several 
ways to design the ceramic transducers to float the 
conventional right inner cylinder. Three of these 
methods are described. 



1. Three-Transducer Model (Fig. 10). 
This model has not been built. Its advantages are that 
maximum excursion can be maintained in the radial 
direction and a smaller excursion, proportional to 
sleeve length, can be maintained in the axial direction 
with all three ceramic sleeves operating at their indi- 
vidual resonant frequencies, A large amount of heat 
is generated but heat transfer is good. Disadvantages 
are the large size, alignment problems, and power re- 
quired for three piezoelectric ceramic sleeves. 

2. Two-Transducer Model (Fig. 11). This 
model was built and tested at MSFC. Its advantages 
are that maximum excursion can be maintained in both 
axial and radial directions. The axial airgap is gen- 
erated by the endplates attached to the outer ceramic 



FIGURE 11. TWO TRANSDUCER SQUEEZE FILM BEARING 
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sleeve; the radial airgap is generated by the inner 
sleeve. Both sleeves are operated at their individual 
resonant frequencies. The two sleeves can be sepa- 
rated by O rings (as shown), by foamed metal, or even 
by gas if the two ceramic surfaces are lapped to close 
dimensional tolerances. The ability of a ceramic 
sleeve to support itself and an inner cylinder was dem- 
onstrated by Mechanical Technology, Inc. Disadvan- 
tages of this model are alignment problems, power 
required for two piezoelectric ceramic sleeves, and 
the relatively poor heat transfer. 


3. One-Transducer Model (Fig. 12). This 
model was built and tested at MSFC, Its advantages 
are that maximum excursion in both axial and radial 
directions is obtained with only one ceramic sleeve. 
This is done by exciting the piezoelectric ceramic 
sleeve at resonant frequency in the radial direction 
and amplifying the smaller axial excursion by means 
of "tuned" endplates, held by flexible mounts. The 



ENDPLATE 

FLEXIBLE MOUNT 


INNER CYLINDER 


CERAMIC SLEEVE 



FLEXIBLE MOUNT 
ENDPLATE 


FIGURE 12. ONE TRANSDUCER SQUEEZE FILM BEARING 

power required and heat generated are less for this 
model than for the other two. The disadvantage of this 
model is the added complexity of making the flexible 
mounts and maintaining them at the correct tension. 

4. Other Configurations. The models shown 
in Figures 13 and 14 were built and tested at MSFC. 
Their purpose was to demonstrate the feasibility of 
spherical and conical shapes for squeeze film bearings. 
Components of both bearings were held together by 
epoxy. 


a. The cup and ball used in our first 
spherical squeeze film bearing came from existing gas 
bearing parts. The ball weighs 1. 81 kg and is 7. 62 cm 



SPHERE 


CUP 


CERAMIC 

SLEEVE 


BASE 


FIGURE 13. SPHERICAL SQUEEZE FILM BEARING 

in diameter. A flat groove was cut on the under 
side of the cup so that it could be cemented to a ce- 
ramic sleeve. The sleeve measures 7. 62 cm in length 
and inside diameter and 0.51 cm in thickness. The 
ball -floats freely at the resonant frequency of 17, 000 
Hz, 80 V, and 1. 2 W power. 

b. The second model of the spherical 
squeeze film bearing is considerably smaller than the 
first and uses a standard steel ball bearing which weighs 



CERAMIC 

SLEEVE 


FIGURE 14. SPHERICAL SQUEEZE FILM BEARING WITH 
AMPLIFYING DISC 
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0.23 kg and is 3. 81 cm in diameter. On this model 
the thickness of the disc supporting the cup was calcu- 
lated to make the disc resonate at the same frequency 
as the ceramic sleeve and to amplify its motion. The 
sleeve measures 5. 08 cm in length and outside diam- 
eter and 0. 64 cm in thickness. The ball floats freely 
at the resonant frequencies of 9,200 and 13,600 Hz, 
60 V, and 0. 54 W power. 


CONCLUSIONS 


Knowledge gained from both experimental and the- 
oretical work suggests that certain limitations of ce- 
ramics be observed. Since the ceramic is a heat 
source, it should be used in a squeeze film bearing in 
such a manner that its own heat is dissipated and that 
it does not interfere with heat dissipation of other 
sources. Ceramics can undergo small but important 
dimensional changes from both self-heating and aging. 
These dimensional changes can result in bearing sei- 
zure when the ceramic is used as a bearing surface 


if proper care is not taken in matching the thermal 
expansion characteristics of this ceramic bearing sur- 
face and the other bearing surface. Finally, the prob- 
lem of providing the ceramic with an adequately hard 
metal surface presents itself if the ceramic is to be 
used as a bearing surface. 

Itappears that most, if not all, of these problems 
are avoided if the ceramic is used only as a source of 
excitation for a conventional metallic bearing surface. 
This arrangement allows use of mechanical amplifica- 
tion of the ceramics relatively small (0.000127 cm 
is typical) excursion by proper suspension. 

The existence of a squeeze film force, under the 
conditions described earlier, has been verified many 
times, both theoretically and experimentally. The 
feasibility of constructing a squeeze film bearing for 
gyro applications able to comply with the imposed con- 
straints of weight, power required, heat generated, 
static and dynamic behavior, etc. , remains to be de- 
termined. As the present time, however, there ap- 
pear to be no insurmountable problems in the con- 
struction of an adequate squeeze film bearing suitable 
for gyro applications. 
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D. RECENT DEVELOPMENTS IN ADAPTIVE TRACKING 
NOTCH FILTER TECHNIQUES by Michael T. Borelli, 
Stanley N. Carroll, and Hans H. Hosenthien 


N65 23785 s 

• • 

Some recent advancements in the development and' 
analysis of the adaptive tracking notch (ATN) filter are 
presented. The ATN fi lter was developed to help stabi- 
lize the structural bending signals encountered in the 
control loop of large aerospace vehicle boosters. First, 
an ATN filter that utilizes a simple differential ampli- 
fier is described, and its equivalence to the ATN filter 
that utilizes a chopper-stabilized operational amplifier 
is presented. Second, the presence of a quadrature com- 
ponent in the output signal of the ATN filter is shown 
to exist, and its effect on the notch characteristics is 
presented. A means of minimizing this effect is given 
along with a theoretical expression that approximately 
describes the magnitude of the in-phase and quadrature 
f components when the input and commutating frequencies 
are the same. Finally, an analytical expression for the 
transfer function of the coupled ATN filter is presented, 
b and its frequency response characteristics are discussed 
and compared to experimental data. Also, the relation- 
ship showing the dependence of the output harmonic 
frequencies of the ATN filter on the number of commu- 
tated capacitors employed and the respective phasing of 
the commutating functions is summarized. 




This article describes some of the work completed 
in the continued development of the ATN filter and in- 
clude s ( i) a mechanization simplification by the sub- 
stitution of differe ntial am plifiers for chopper stabi- 
lized amplifiers, (2) investigations of compensation 
for the filter’s quadrature component of the structural 
bending mode signals, and (3) analytical studies that 
include (a) the derivation of a suitable transfer 
function for the RC commutated network utilized in the 
ATN filter, (b) the determination of the sideband sig- 
nals generated by the RC commutated network, and (c) 
the calculation of the optimum value of feedback adjust- 
ment for the RC commutated network. 


NOMENCLATURE 


*2k 


D 


phasing parameter (equal to 1 or 2) 
Fourier coefficients of output signal 
attenuation of filter at frequency 2kc o 


b^ Fourier coefficients of output signal 


c 


special function 
e voltage function, volts 

G transfer function 

h(t) , H(s) linear network, (time, frequency domain) 

j 


\! -l 


INTRODUCTION 


i The adaptive tracking notch (ATN) filter was de- 
veloped to help stabilize the structural bending signals 
encountered in the control loop of large aerospace ve- 
hicle boosters.^ In two previous articles (Refs. 1 and 
2) , the ATN filter was described in detail and results 
from its use in analog computer simulations of the Sat- 
urn V flight were presented. Although an ATN filter 
breadboard autopilot has operated very satisfactorily 
in all laboratory tests and simulations conducted thus 
far, the present design is by no means considered op- 
timum. 


k integer 

K potentiometer setting 

K’ feedback gain 

m integer 

n integer 

N number* of commutated capacitors (also 

equal to number of paths) 

P (t) commutating function 
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P ,P, 

m k 


q 

s 

t 

T 

o 

x,X 


y,Y 


a 

o 


P 

y 

6 

^>2k 

T 

CO 


Fourier coefficients of commutating 
functions 

integer 

complex variable, radians/second 
time, seconds 

period of commutating functions 

input signal (time, frequency domain) 

output signal (time, frequency domain) 

reciprocal of equivalent time constant, 
1/seconds 

normalized frequency 

dimensionless parameter 

phase angle, radians 

time shift, seconds 

phase shift of filter at frequency 2kco^ 

time constant, seconds 

angular frequency, radians/second 
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FIGURE 1. RC COMMUTATED NETWORK FOR CHOPPER- 
STABILIZED OPERATIONAL AMPLIFIER 
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co angular frequency of commutating func- 
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FIGURE 2 . RC COMMUTATED NETWORK FOR 
DIFFERENTIAL AMPLIFIER 

differential amplifiers are shown in Figures 1 and 
2, respectively. The derivation of the equation repre- 
senting the commutated network output voltage e 2 of 
Figure 1 is obtained in the following way. 


DIFFERENTIAL AMPLIFIER ATN FILTER 


In the ATN filters developed previously, chopper- 
stabilized operational amplifiers were used in the RC 
commutated networks. To simplify the mechanization 
of the ATN filter, investigations were conducted to de- 
termine if differential amplifiers could be used in lieu 
of the chopper-stabilized operational amplifiers. Both 
analytical and experimental investigations have shown 
that equivalent performance can be* obtained. 

Schematic diagrams for the RC commutated net- 
works as used with chopper-stabilized operational and 


ii + i 2 + i 3 =0 

where u = and i 3 = e 2 
K l 

therefore 



The voltage across the commutated capacitors 



( 1 ) 


(2) 


is 


e 2 = E 

n=l 


f / P n *' 


dt 


(3) 
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where P n is the square-wave commutating functions 
for each capacitor. 

Now substitute equation 2 into equation 3 to obtain 
the desired relation 


e 2 


4 P 
n=i 


r _ 

e l 

f 1 + J_\ 

f p 

; n 

^ +G2 

[ R 2 R 3 J 


— 



dt . (4) 


The equation for the output voltage of the RC commu- 
tated network of Figure 2 is obtained in much the same 
manner. 


* 1=12 + h 


(5) 


where ij = 

therefore 

e l 

■* 2 = ¥ 4 


e t - e£ 

andi 3 = 7T ; 


Ri 


1 [ 1 
r 4 r 5 


(6) 


The voltage across the commutated capacitors is 

(7) 


>'* - t C* / P n dt - 


n=l 


Now substituting equation 6 into equation 7, 


•* “ t ^ 

n=l 


®i 


- e 2 I — + ~ 

2 \R 4 R, 


dt. 

(8) 


By the proper selection of the resistors, equations 4 
and 8 can be seen to have the same form. For ex- 
ample, in equation 4 let R 2 = Rj and in equation 8 let 
R 4 = Rj and R 5 = R 3 ; then equations 4 and 8 have the 
same form except for some important differences in 
signs. The differences in signs are resolved when the 
two circuits are utilized in a notch filter circuit as 
shown in Figures 3 and 4. 

From Figure 3 theoutputof the ATN filter utiliz- 
ing a chopper-stabilized operational amplifier com- 
mutated network is 



FIGURE 3. ATN FILTER WITH CHOPPER— STABILIZED 
OPERATIONAL AMPLIFIER 

Substituting equation 9 into equation 4 to eliminate the 
voltage e 2 and letting R 2 = R t , the output voltage of the 
ATN filter is obtained. 


e *H e ‘ + £ / P n 

n=i 6 


R 3 ^ 

e i + (1 + zr) e 3 
Rl ^ 


dt). (10) 


Theoutputof the ATN filter utilizing the differen- 
tial amplifier commutated network of Figure 4 is 


e 3 “ e 2 ” e i 


(U) 


Substituting equation 1 1 into equation 8 to eliminate the 
voltage e£ and letting R 4 = R t and R 5 = R 3 , the output 
voltage of the ATN filter is obtained. 


e «H e ‘ + i W f P n 


n=l 


R 


1 


e 4 + (1 + — )e 3 


dt; • (12) 


c - , 



K- r 


C x X p > 

i- --c 


c j'A B. 


e 3 = “( ©i + e 2 ) . 


(9) 


FIGURE 4. ATN FILTER USING A DIFFERENTIAL AMPLIFIER 


33 


ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


Thus the equivalence of the two ATN filter circuits is 
established. Inaddition to the simplification achieved 
by replacing a chopper stabilized amplifier with a dif- 
ferential amplifier, one less amplifier is required for 
the ATN filter of Figure 4 than for the ATN filter of 
Figure 3. 

To verify these analytical results experimentally, 
two equivalent circuits were built and tested. Although 
the values of the circuit parameters (R t , R 2 , C, etc.) 
for the two circuits are not equal, the values of R 3 C 
and R/Rt are equal. Normalized amplitude frequency 
response plots of the two circuits tested are shown in 
Figure 5. Good correspondence (approximately 5 per- 
centerroror less) between the two circuits is obtained 
below and around the notch. By trimming the values 
of R 3 C and R 3 /R t in one of the circuits, a closer cor- 
respondence between the dynamic performance of the 
two circuits can be obtained. 



FIGURE 5. ATN FILTER AMPLITUDE, FREQUENCY RESPONSES 
FOR FOUR CAPACITOR NETWORKS 

The equivalence between the differential and the 
chopper-stabilized amplifier circuits have thus been 
shown, but some specific limitations are imposed on 
networks employing transistorized differential am- 
plifiers. These limitations are attributed to the inher- 
ent characteristics of the transistorized differential 
amplifier. For example, the temperature sensitivity 
of the amplifier is proportional to the value of the in- 
put and feedback resistors. This means that smaller 
resistances shouldbe used and therefore larger capac- 
itances to obtain the same dynamic performance. 

QUADRATURE COMPONENT OF THE 
CONSTRUCTED BENDING MODE SIGNAL 

The ATN filter attenuates a bending mode signal 
by generating a reversed phase bending signal which, 


when added to the control sensor output signal, will 
suppress the amplitude of the bending mode component. 
Complete suppression cannot be realized because the 
generated signal does not have an exact 180-degree 
phase relationship with the input signal. Based on the 
transfer function of the commutated network defined 
and discussed in the next section, the frequency re- 
sponse locus in the complex plane is approximately a 
circle as shown in Figures 6A and 6B. The input to 
the commutated network is the vector x(w ) , and the 
output component at the same frequency as the input 
is Ky(a; n ) , where n = 0 and 1 (co 1 < a> 0 ) . . . and K is 

the feedback potentiometer setting or the gain adjust- 
ment. In Figure 6A, co 0 is the commutating frequency. 
The angle 6 is greatly exaggerated to show clearly the 
presence of a quadrature component Ky (^o)q °f the 
bending mode signal at coq. 

The notch filter output R(oo n ) is obtained by the 
vector addition of x(co n ) andKy(c*; n ). Figure 6B shows 

that the notch will always occur below the commutating 
frequency o) 0 , or at the frequency Wj as shown. Also, 
the smaller the angle 6, the smaller the quadrature* 



FIGURE 6A. VECTOR DIAGRAM AT THE 

COMMUTATING FREQUENCY o>. 



FIGURE 6B. VECTOR DIAGRAM AT THE 

FREQUENCY OF THE NOTCH (0 1 
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component will be and the deeper the notch. If 6 = 0, 
the notch depth would be infinite, where the notch depth 
is defined as the reciprocal of the normalized gain at 
the notch frequency. The factor K, the feedback po- 
tentiometer setting, is adjusted to minimize the re- 
sulting vector R(oj 0 ) . The minimum value of vector 

R(o; 0 ) occurs when R(co 0 ) is perpendicular to the 
Ky(co q) vector and will yield its minimum at the com- 
mutating frequency u; 0 . 

One technique successfully employed in reducing 
the effects of the quadrature is a phase lead network 
in the path containing the commutated network as shown 
in Figure 7 A (or equivalently, to include a phase lag 
circuit in the unity gain path, Figure 7R) . The phase 
shift rotates the vector Ky(a; 0 ) toward the negative 
real axis, thereby reducing the phase angle 6. The 
angle 6 is small enough at the commutating frequency 
cu 0 that the phase shift network does not appreciably 
alter the length of the vector, but does provide the 
proper phase. Once the amount of phase shift has been 



FIGURE 7A. UNCOUPLED ATN FILTER WITH 

LEAD QUADRATURE COMPENSATION 



FIGURE 78. UNCOUPLED ATN FILTER WITH 

LAG QUADRATURE COMPENSATION 


determined, it is always possible to change the com- 
mutated network gain slightly, if necessary, to yield 
the maximum notch depth. When the angle 6 is small 
( 5 degrees) , the compensated phase shift will vary 
only slightly throughout the first bending mode fre- 
quency range. This characteristic provides good im- 
provement to the notch depth throughout the desired 
tracking range. 


Two RC commutated network configurations, 
coupled and uncoupled (Ref. 1), have been studied and 
their equivalent circuits are shown in Figures 8A and 
8B, respectively. For the uncoupled case , an analyt- 
ical expression for the’quadrature shows it to be en- 
tirely dependent on the RC time constant and the com- 
mutating frequency. This can be established by first 
deriving the fundamental output of a four capacitor un- 
coupled' case. For a sine wave input of frequency u> 0 



FIGURE 8A. EQUIVALENT DIAGRAM FOR 

COUPLED RC COMMUTATED NETWORK 



FIGURE 8B. EQUIVALENT DIAGRAM FOR 

UNCOUPLED RC COMMUTATED NETWORK 
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and maximum amplitude x 0 , the fundamental output can 
be shown to be 


y(t) = ~ 8k o 

X 0 - 7T 2 

8k 0 

sin co 0 t + — p 

7T 


°° 2( 4k 2 + 1) A 2k cos <£ 2k 

z - 

k=l 


oo 8k A 2k sin <p 2k 

z 

k=i 


(4k 2 - 1) 


(4k 2 - l) 2 

cos Cj 0 t ( 13) 


where 

A _ 1 

A 2k — • - — 

\/l + 4k 2 ( RC) 2 co 0 z 

0 2 k = arc ten 2kRCo) 0 

k 0 = dc gain of smoothing filter (first order) 
o; 0 = commutating frequency. 


A 2k and 0 2k are the attenuation and phase of the inte- 
grator of Figure 8B with the negative resistive feed- 
back that acts as a first order smoothing filter meas- 
ured at frequency co = 2kcu 0 (only even harmonics of a> 0 
are present at the smoothing filter because of demod- 
ulation by the first multiplier) . The magnitude of the 
quadrature component is the coefficient of the cosine 
term and the magnitude of the in-phase component is 
the coefficient of the sine term. The signs associated 
with equation 13 verify that the vector representing 
y(t) will be as indicated in Figure 6. The angle 6 can 
be computed by taking the inverse tangent of the ratio 
of the magnitudes of the quadrature to in-phase com- 
ponent. The two infinite series can be expressed in 
closed form (Ref. 3) after substituting for A 2k and 
0 2 k* The final result is 


6 = tan 


-l 


JL 


t Trty 2 + 1) -4y coth 2yJ 


> (14) 


7T (y 2 + 1) + 2y(y 2 - 1) coth ~~ 
^ 2y 




where y = RCu> 0 = ro) 0 . 


If the commutating frequency is much greater than the 
break frequency of the first order smoothing filter, 
i. e.,TCO 0 »l , a good approximation to equation 14 is 


6 = tan 


0. 188y 
1 + 0.405^ 


(15) 


For y = 1, 6-17 degrees; for y = 2tt, 6 = 4. 2 degrees. 
This shows that if the break frequency is decreased 


co o 

from co 0 to — the angle 6 is reduced from 17 degrees 
27 r 

to 4 degrees. The above results indicate how the notch 
frequency depends on the RC time constant through the 
magnitude of the phase angle 6. 


Although equation 14 was derived for four capaci- 
tors, it can be shown (Ref. 4) that the equation holds 
true regardless of how many resistor-capacitor pairs 
are employed, provided the RC values are the same. 
In contrast, for the coupled configuration, the angle 6 
becomes smaller as the number of capacitors is in- 
creased. As stated previously, 6 must be reasonably 
small to obtain a good notch depth. The reason a 
correction circuit is required in the uncoupled case 
to reduce the quadrature component as opposed to using 
a large time constant is that the larger time constant 
yields a notch too narrow to satisfactorily suppress 
the bending mode signal. The effects of the time con- 
* stanton the notch filter frequency response are shown 
in Figures 9A and 9B for the coupled case. For a 



FREQUENCY (Hz) 


FIGURE 9A. ATN FREQUENCY RESPONSE, AMPLITUDE 
—THEORETICAL 
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FIGURE 9B. ATN FREQUENCY RESPONSE, PHASE 

-THEORETICAL 
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large time constant, it was found that the response of 
the uncoupled case is approximately the same as that 
of the coupled case. 

The quadrature signal in the coupled configuration 
is much smaller than in the uncoupled configuration. 
This is attributed to the cross-coupling terms appear- 
ing in the equations as a result of the feedback loop 
with gain K T . The block diagrarn of the coupled case 
can be represented as the product of the uncoupled 
case and a network containing only the cross-coupling 
terms. Based on this equivalent block diagram, it is 
concluded that the cross- coupling terms tend to rein- 
force the in-phase component and attenuate the quad- 
rature component, Although a general expression for 
6 has not been derived for the coupled case ( similar 
to equation 14) , it is possible to calculate the value of 
6 from the transfer function (equation 16) for given 
valuesof RC, K T , and co 0 , In the present design of the 
coupled ATN filter, the angle 6 has been calculated to 
be less than one degree. 


can be written as the sum of the transfer function and 
the spectrum of the sideband frequencies. 


The transfer function for the circuit shown in Fig- 
ure 10 is given by equation 16 (Ref. 4). 


Y„(s) 

xTST 


= G(s) 



+ «o D 


tan (3 
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(16) 
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D = 


K ! o 


s + K^q 


7TS ( S + K^fl) 
8K’ olh co n 


- coth 


7r ( s + K’a 0 ) 
8u; n 


coth 


7TS 
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ANALYTICAL STUDIES OF 
RC COMMUTATED NETWORKS 

This section reviews three important character- 
istics of RC commutated networks as applicable to the 
ATN filter. First, the transfer function and the side- 
band spectrum for the commutated network in the ATN 
filter are discussed. Second, it is shown that the side- 
band frequencies depend only on the number of com- 
mutated capacitors employed and the phasing between 
the commutating driving functions. Third, a deriva- 
tionfor the optimum value of the feedback gain is given. 


Transfer Function' For The Coupled Case 

If the RC commutated network in the ATN filter is 
subjected to a single -frequency input, the total output 
signal is comprised of a primary frequency component 
plus certain sideband, or complementary, frequency 
components. The primary frequency is the same as 
the input frequency; the sideband frequencies occur at 
multiples of eight times the commutating frequency 
plus and minus the primary frequency for four com- 
mutated capacitors. An analytical expression has been 
derived (Ref. 4) for the output frequency spectrum of 
the commutated network. The term ’ transfer function” 
is used here to denote the ratio of the output compon- 
ent at the primary frequency to that of the input driv- 
ing signal. Using this definition, the output spectrum 


t “ RjC 
4 


K = potentiometer setting 
R i 

K' = feedback gain, — K 
R f 

cj 0 = commutating frequency 


Y 0 = primary frequency component. 


c c c c 



FIGURE 10. FOUR CAPACITOR COMMUTATED NETWORK 
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The equation relating the sideband amplitudes to 
the input is 

Y^(jco + jq8w 0 ) 

X(jco) “ 



jcj + K T a 0 


where q is all integers except zero, G(ju; + jq 8co 0 ) is 
the value of the transfer function (equation 16) evalu- 
ated at frequency (w + 8q c o 0 ) , and Yq( jcu 4- jq8u) 0 ) is 
the sideband component at frequency (w + 8qcu 0 ). 

The term D in equation 16 is present because the 
sideband frequencies in theoutputare fed back through 
K\ which results in an output component at the pri- 
mary frequency. For frequencies between zero and 
2. 5 cu 0 , D can be neglected without introducing more 
than approximately five percent error (Ref. 5). 

A comparison between the theoretical and experi- 
mental frequency responses for the notch filter is 
shownin Figure ii. There are two important param- 
eters in the notch filter: the RC time constant and 

the feedback gain K\ The time constant is used to 
control the width of the notch, as shown in Figures 9A 
and 9B. The time constant also slightly affects the 
frequency at which the minimum gain value (notch) 
will occur as previously discussed. Variations in the 
feedback gain K’ about its nominal value have negligi- 
ble effect on the notch width, but mainly determine 
the minimum gain value at the notch frequency. 



0.1 02 0 4 06 08 10 2 0 4 0 6 0 80 100 


FREQUENCY (Hz) 

FIGURE 11. ATN FREQUENCY RESPONSE, AMPLITUDE 
-THEORETICAL AND EXPERIMENTAL 


Determination of Sideband Frequencies 


As indicated in equation 17 the sideband frequen- 
cies depend on the term q8o> 0 where the number eight 
results from two factors: the number of capacitors 

utilized and the respective phasing between the com- 
mutating driving functions. The dependence of the 
sideband frequencies on these two factors can be seen 
by investigating the fundamental equation describing 
the N capacitor configuration shown in Figure 12, For 
N = 4, the circuits in Figures 10 and 12 are equivalent. 
This equation (Refs. 4 and 6) is 


-Y(s) = Z E 




P P H( s-jc^ok) [X( s-jcook-jcoom) 
m k 

+ K T Y( s-jcuQk-jcoom)^ ( 18) 


where 



Fourier coefficients of the square wave 
commutating functions 


H(s) 


-1 

TS 


T 0 = period of commutating functions 
N = number of commutated capacitors 



aN 

a = phasing parameter (equal to 1 or 2) . 



FIGURE 12. EQUIVALENT CIRCUIT FOR 
N COMMUTATED CAPACITORS 
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The exponential factor occurs in equation 18 be- 
cause the nth commutating function is shifted by X(nrl) 
seconds relative to P t (t), i.e. P (t) = P^t-A. (n-1)] . 
The significance of the frequency shift in the argu- 
ment of the input term, X, in equation 16 is that the 
sideband frequencies will be at o) + u> 0 (k + m) . It will 
be shown that equation 18 is identically zero unless 
(k + m) is a multiple of aN, thereby establishing the 
dependence of the sideband frequencies on the phasing 
parameter and the number of commutated capacitors. 


The summation over n is seen to be a geometric 
series. Since 


corX = -~r , then 
v aN 

2 7T 

N -j (n_1) (k + m) 


S e 

n=l 


= 1 + e 


2tt .. , 

-i xs (k + ra) 


using a complete circle (2tt radians) or a half circle 
( 7 r radians) for phasing the commutating functions. 
The admissible values of q depend upon whether N is 
even or odd and the two possible values of "a. M Be- 
cause the commutating functions are square waves, 
k + mis always an even number; thus the number q(aN) 
must also be even. If N is even, then either value of 
the phasing parameter "a" results in all integers for 
q being admissible. Likewise, all values of q are ad- 
missible if "a" is even and N is odd. However, when 
both N and "a" are odd, then q must assume only even 
values. The final result is that for odd N the sideband 
frequencies will be at co + q(2N) , and for even N the 
sideband frequencies will be at co + q(aN) , q = 0, ± i, 
±2, ... . For the ATN filter, a = 2 was chosen be- 
cause it was desired to eliminate as many low sideband 
frequencies as possible. 


2 „ Optimum Value of Feedback Gain K f 

+ e J ^ (k + m) 2 + _ +e -j (k+m) (N-1). (19a) 


For k + m = q(aN) , where q assumes certain in- 
tegral values, each of the N terms in the series is 
unity. If k + mX q(aN) , the series can be expressed 
in closed form as 


N .27 T . . W1 x 

l e - J ^ (n - 1)(k + m) 

n=l 


— j 2 7T (k+m) 


1 - e 


— j 2 7r ( k + m) 
aN 


(19b) 

^ k + m , , , , 

Because — was assumed not to be an integer in 

aN 

equation 19b, the denominator is non-zero. In con- 
trast, the numerator is identically zero since + n - 
is always an integer. 


Therefore 


N . 27r 

Yj e' ] aN ( + m) (n- ) _ N for k + m = q(aN) (19c) 

n=l 

= 0 otherwise. 

With the constraint k + m = q(aN) , and replacing 
s by jcu in equation 18, it is seen that the sideband fre- 
quencies will be at c o + q(aN) u? 0 . The two values, 1 
and 2, of the phasing parameter fT a M corresponds to 


It was stated in a previous section that the phase 
angle 6 (Fig. 6) is approximately zero for the coupled 
configuration of the commutated network. If this angle 
is assumed to be zero, itis apparent thatexact cancel- 
lation between the bending mode component contained 
in the sensor signal and the reversed phase signal de- 
pendsonlyon the feedback gain K\ Using steady state 
analysis, an expression for the optimum K* can be de- 
rived and is seen to depend only on the order of har- 
monics in the commutated network output. 

With reference to Figure 13, the average value of 
the input to each of the four integrators must be zero 



FIGURE 13. EQUIVALENT CIRCUIT FOR FOUR 
COMMUTATED CAPACITORS 
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in a steady state condition. The path containing the 
first integrator will be considered since (t) is in 
phase with the input signal. The Fourier series for 
Pi (t) is 


p, <«> - - E 

n=0 


sin ( 2n + 1)gl? 0 t 
2n+i 


The input and output signals are given as 
x(t) = x 0 sin co 0 t 


y(t) = Z (a^ cos kco 0 t + sin kco 0 t) . 
k=0 


( 20 ) 


( 21 ) 

( 22 ) 


The amplitude of the harmonic components given by 
equation 25 corresponds to the Fourier coefficients 
"b" in equations 22 and 24. Equation 26 implies that 
bg , b 5 , b tl , b 13 , ... are. zero since the third, fifth, ... 
harmonics are not present. Also, the amplitude of 
the seventh, ninth, . .. harmonics are proportional to 
i/7, 1/9, ..., respectively. Because the output is 
180 degrees out of phase with the input, each of the 
Fourier coefficients, b^ + in equation 24 is nega- 
tive as shown in Reference 4. Using equation 26 in- 
equation 24 gives, for the four capacitor case, 


K’ 


1 



1 

+ 15 2 



= 0.95. 

(27) 


The average value of the signal into the first integrator 
will be the dc component that results from multiplying 
Pj(t) by x(t) + K r y(t). Since theproductof a sine and 
a cosine term contributes no dc component, the net dc 
component results from multiplying the corresponding 
sine terms of like frequencies. It is seen that the total 
dc is 


2x 0 

dc = + 

7 r 


2K' 

7 r 


OO 


E 

n=0 


b 2n + 1 
2n + 1 * 


(23) 


In general, if ru) ( y»l, a good approximation for the 
feedback gain is 


K’ 


1 


Z 



(28) 


where the index n corresponds to the order of harmon- 
ics which are present in the output. For the two capac- 
itor case, all odd harmonics are present and equation 
28 yields 


or since equation 23 must be zero, 


K' 


~ x o 


oO 


E 

n=0 


b 2n + 1 
2n + 1 


(24) 


Using equations 16 and 17, for s = jco 0 , and assuming 
x 0 = -b 1 to have exact cancellation, it can be shown that 
the amplitudes of the harmonic components are 


Y 

q 


cu 0 (l-8q) 


x 0 


O?0 ( 1 - K ') 

K r a Q + jco 0 ( l-8q) 


(25) 


If the time-constant is sufficiently large that a oK f 
can be neglected relative to cu 0 ( l-8q) , and by knowing 
that K' is close to unity, then equation 25 can be ap- 
proximated by 

Y Oo<l-8q)] t 

x 0 a > q = ±1 » ±2 * ••• < 26) 


K' = -dr \ J . (29) 

Z (2n+l) 2 
n=0 

Equations 27 and 29 are in excellent agreement with 
experimental results. 


CONCLUSIONS 


The results obtained from these investigations 
have (1) permitted a simplification of the mechaniza- 
tion of the ATN filter for some applications, (2) im- 
proved the performance of the ATN filter, and (3) 
achieved an improved understanding of the operation 
of RC commutated networks. The derivation of the 
transfer function for the RC commutated networks is 
a definite extension of the state of the art for this type 
of circuit. The use of this new analytical tool will be 
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reflected in a considerable redaction in the required 
experimental efforts in further development work. 

Future development efforts will include replacing 


the present mechanical relays for the commutated ele- 
ments with electronic switches. Investigations will 
also be made to improve the identification function or 
the tracker system of the ATN filter. 
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E. REAL TIME ANALOG COMPUTER SIMULATION OF 
THE ELASTIC SPACE VEHICLE INVOLVING REAM 
CELL TRANSMISSION MATRICES by George L. von 
Pragenau 

^ N65-23 786. j, 

A flying space vehicle behaves like a freely moving 
beam responding with translation, rotation, and oscilla- 
tions to external forces. The latter especially compli- 
cates the analysis and design of a control system used to 
stabilize the attitude of a space vehicle. Certain meas- 
ures yield a stable operating control system; however, 
its design requires agood knowledge of the bending modes 
of the space vehicle. 


The electronic model described is one method to 
obtain such knowledge or to simulate the space vehicle's 
response in real time on an analog computer, together 
with actual control system hardware. The method allows 
observation of traveling waves; e.g., a disturbance 
travels in approximately one half of a second from the 
tail to the nose of the Saturn V space vehicle. Unlike 
other methods, this process considers the effect of dis- 
tributed damping of the vehicle's structure and therefore 
represents more realistically the dynamic effects of 
elastic vehicles. An electronic model of the Saturn V 
vehicle was successfully operated together with a con- 
trol loop to stabilize the vehicle's attitude. The lateral 
deflection of the beam was displayed by the lateral ac- 
celeration at nine stations. Three bending modes or the 
standing waves resulting from sinusoidal excitation, a 
traveling wave resulting from unit step excitation, and 
a reaction of the Saturn V vehicle in the controlled 
operation mode when excited by a unit step wind force 
are shown. The frequency response of the model is 
plotted for sinusoidal excitation at the tail of the vehi- 
cle for the frequency range of 0.5 to 10 Hz. A block 
diagram of the vehicle model and its matrix formulation 
are included. K 


y 


INTRODUCTION 

The dynamic behavior of a flexing beam can be 
described by partial differential equations. Partial 
differential equations approximated by finite differ- 
ence equations can be solved by analog computers as 


shown by several authors. The method treated in this 
article differs from these solutions by applying a ma- 
trix formulation not directly related to the partial dif- 
ferential equation. The concept was developed by the 
Astrionics Laboratory to simulate flexible space ve- 
hicles while maintaining good physical insight. 

The simulation utilizes standard analog computer 
equipment to form an electronic model of a space ve- 
hicle and allows real time simulation including the 
hardware of the flight control system. This method 
simulates transmission characteristics, distributed 
structural damping, distributed wind forces, and mul- 
tiple sensor locations. Transfer functions and bend- 
ing modes are readily obtainable and parameter studies 
are possible without difficulty. The flexibility of the 
method will save time and cost in developing control 
systems and in studying control concepts which are 
rather complex. 

The control system contains attitude sensors, 
signal conditioners and filters, amplifiers, and a 
thrust vector actuator system. Attitude sensors are^ 
usually placed at accessible locations in a vehicle and 
detect local attitudes which contain the attitude of the 
whole vehicle plus the undesirable effect of the vehi- 
cled flexibility. The latter needs special attention 
and if not properly processed could result in wild os- 
cillations finally leading to failure. The vehicle and 
the control system compose a feedback loop; if the 
loop gain is larger than one, instability is possible. 
Signal conditioners are utilized to shift the phase angle 
of the sensor signals to damp oscillations artificially, 
and filters are used to attenuate the gains at the res- 
onance frequencies of the space vehicle. 

Commonly the methods of Rayleigh-Ritz or Stodola 
(Refs, i, 2, and 3) are employed to determine bend- 
ing mode shapes and bending frequencies, thus de- 
scribing the dynamic characteristics of a flexible 
space vehicle. Structural gains at bending resonances 
are obtained to design the control system by assuming 
a certain overall damping. The flexible space vehicle 
is simply replaced by damped mass-spring resonators 
for each natural frequency. 

Another way to find dynamic vehicle characteris- 
tics is ground vibration tests (Ref. 4) on a full size 
space vehicle. By suspending it as softly as possible 
to simulate free flight conditions, valuable data can 
be gained, e.g., damping coefficients and transfer 
functions at different sensor locations. The test ob- 
ject is vibrated sinusoidally at different frequencies 
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by a shaker device. Since this method is a test on the 
actual vehicle, it has a certain value, especially in 
producing damping data. Checking on local resonance 
effects is another advantage of ground vibration tests; 
for very complex systems such as clustered space ve- 
hicles, it is probably the only way to obtain a true 
picture of dynamic effects. However, such tests are 
costly and difficult to schedule. 

Data obtained by these conventional methods are 
utilized in an analog computer setup. The computer 
reproduces resonance effects to checH out control 
hardware under simulated dynamic conditions. This 
is a less direct approach than using an electronic 
model of the space vehicle. The term. model applies 
here more to analog copying than to a less accurately 
general representation of dynamic qualities. The 
model directly simulates the space vehicle based only 
on mass, stiffness, and damping information; the elec- 
trical signals of the model can be fed directly to the 
control hardware in checkout. With respect to the 
number of amplifiers of the conventional and the pro- 
posed analog computer arrangement for a similar task 
"(six resonances and deflections at six different sta- 
tions of the vehicle) , no difference of computer vol- 
ume is required; the volume for each case is approx- 
imately 32 amplifiers and 12 integrators. However, 
for special applications where fewer resonances and 
deflections are required, the conventional analog com- 
puter setup employs a smaller computer volume (for 
4 resonances and 1 deflection, only 13 amplifiers and 
8 integrators) , while the electronic model must still 
maintain a larger volume to avoid inaccuracy. 

A possible classification of dynamic data produc- 
ing methods appears as follows. The ground vibration 
test represents the empirical method, the electronic 
model represents the semi-empirical, and the bending 
mode technique represents the purely analytical meth- 
od. The electronic model has a dual function; it pro- 
duces dynamic data and can be used to simulate the 
overall dynamic behavior of the control system. The 
method does not require elaborate means to produce 
data needed for a conventional analog computer sim- 
ulation. 

Stiffness, damping, and mass distributions are 
essentially all that is needed to generate certain ma- 
trix expressions used for the electronic model pro- 
graming. From all the data, damping is the least well 
defined information; equivalent viscous damping is 
employed to approximate dynamic test results. How- 
ever, it is anticipated that damping-oriented tests 


will clarify the damping characteristics much better 
in the near future. By dividing the beam -like object 
into cells that are represented by matrices, trans- 
mission characteristics are preserved. The ideali- 
zation of the classical methods (Refs. 2 and 3) that 
consider only "conservative" systems in resonant 
modes was avoided. A good physical insight was ob- 
tained without too much abstraction, thus aiding the 
usage of the electronic model. Therefore, the direct 
approach of. modeling electronically results in good 
adaptability. Mass, elasticity, and damping data can 
be changed readily by potentiometers, and the results 
are obtained at analog computer computing rates. 
Scheduling is not difficult and expenses seem to be 
much lower than for any of the other methods. 

The electronic model was tried at first as a rather 
crude representation of the Saturn V space vehicle by 
lumping the mass characteristic of each stage about 
its particular center of gravity. The space vehicle 
was divided into four mass fields and three intercon- 
necting flexibility fields. Each mass field contained 
information about mass, moment of inertia, and cen- 
ter of gravity. Each flexibility field described the 
flexibility contributions of elongation, bending, and 
shear stiffness. This lumped mechanical model was 
then described by matrix equations (Ref. 1) which 
were further used for programing an analog computer. 
Because the variables were available at the bounds of' 
each field, boundary conditions were easily satisfied. 

Other electronic model techniques of elastic beams 
and partial differential equations are found in the lit- 
erature (Refs. 5 through 9). These techniques are 
based on the approximation of partial differential equa- 
tions by finite differences which make boundary con- 
ditions difficult to apply. 

The work of Charles Elms (Ref. 10), who built 
the signal conditioner to show the mode shapes, is 
gratefully acknowledged, as is the encouragement of 
my colleagues. 


NOMENC LATURE 


A cross section of beam structure 

a 0 dc gain of (j> filter of control loop (Fig. 7) 
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dc gain of (p filter of control loop ( Fig. 7) 

D. flexibility matrix of the beam structure be- 

tween stations i and i + i (equation 12) 

DT 1 inversed matrix of D. 

i l 

E modulus of elasticity 

F. generalized force on flexibility segment at 

station i, Laplace transformed 

F. same force not transformed (equations 2 

1 and 3) 

G modules of rigidity (equation 12) 

G. mass matrix at station i (equation 4) 

G. -1 inversed matrix G^ 

g center of gravity eccentricity of mass m 

from a station (Fig. 1) 

g 2 dc gain of y filter of control loop (Fig. 7) 

I moment of inertia of lumped mass 

J cross-sectional moment of inertia of beam 

structure 

k rigidity efficiency of cylindrical shell (equa- 

tion 12) 

L field length of beam cell (equation 5) 

L i lever arm matrix from station i to i + 1 

(equation 5) 

L.^ 1 inversed lever arm matrix (equation 5) 

L i transposed lever arm matrix 6. 

L^ 1 inverse of transposed lever arm matrix 

M moment, component of generalized force 77. 

(equations 2 and 3 and Fig. 1) 

m mass of lumped mass at a station (equation (p 

4 ) 


p Laplace's differential operator 

s variable station on vehicle axis (equation 12) 

T longitudinal force, component of general- 

ized force (equations 2 and 3) 

T c controlled thrust force ( Fig. 7) 

U. Laplace transformed deformation of flexing 

beam cell at stations i (equations 10 and 11) 

V lateral force or shear force, component of 

generalized force (equations 2 and 3) 

generalized wind force vector attacking 
mass G- at station i 

• • 

X i space-fixed acceleration vector at station i 

(equation 1) 

x space-fixed position component (equation 1) 

y space-fixed position component (equation 1) 

Z. body-oriented position vector at station i, 

Laplace transformed 

Z. body-oriented acceleration vector at station 

i (equation 8) 

gimbal angle of controlled rocket motors 
(Fig. 7) 

body-oriented acceleration signal for accel- 
erometer control loop (Fig. 7) 

space-fixed rotation angle, component of 
position vector (equations 1 and 7) 

rotation angle, a component of body-oriented 
position vector Z^ Laplace transformed 
(Figs. 3, 4, and 5) 

lateral deflection, a component of body- 
oriented position (Figs. 3,4, and 5) 

angular control signal of space vehicle's 
angular position (Fig. 7) 
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FORMULATION 

Generally, rotations of a space vehicle are rela- 
tively slow; therefore, it is sufficient to consider only 
planar motions for the design of the attitude stabiliza- 
tion. Consequently, only these motions will be treated 
here, while the circumstance is appreciated that the 
resulting equations are much simpler than those of a 
three dimensional formulation. 

First, equations will be derived relating the 
space-fixed planar accelerations (two translatory, one 
angular) of a rigid mass (Fig. 1) to external forces 
(two normal forces, one moment). As indicated, ac- 
celeration is represented by a vector consisting of 
three components; one in the x direction, one in the y 
direction, and one being the angular acceleration 
(equation i) . The force is also a vector composed of 
three components; one in the longitudinal direction, 
one normal to the longitudinal axis, and one being the 
moment (equation 2). As defined, the terms accel- 
eration and force are applied generally to simplify the 
terminology in connection with the vector equations. 
The first equation relates the vectors of the forces 
and accelerations appearing on a rigid body ( Fig. 1) . 
The equation can be considered a dynamic equilibrium 
condition between the external forces Fj and F 2 and 
the D'Alembert's forces of the mass (equation 6). 

A detailed derivation of the equations is beyond 
the scope of this article; therefore, only the essential 
equations and steps are shown. Equations of a flexi- 
ble but massless beam will follow the equations of the 
massy but rigid beam. 



FIGURE 1. FORCES AND DISPLACEMENTS ON THE ELONGATED 
RIGID BODY 


The acceleration vector is defined by 



The force vectors are defined by 



The effect of the D'Alembert forces is the pro- 
duct of the acceleration of equation 1 projected into 
body-oriented coordinates and a mass matrix G 2 . 



0 

m 

mg 



( 4 ) 


The components oi G 2 contain the mass m of the 
body, the static moment mg about point 2 (Fig. 1), 
and the moment of inertia I + mg 2 about point 2. 

The beam functions also jis a lever arm of length 
L, thus relating the forces Fj and F 2 . This is ex- 
pressed by the lever arm matrix 



For very slow motion, the equilibrium condition F 2 = 
Lf 1 Fj is valid; for faster motions this relation must 
be modified by subtracting D’Alembert’s forces 

M 

( g 2 0 2 x 2 ) . 

• m 

F 2 - Li 1 Fi - G 2 n 2 X 2 . (6) 
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The matrix £2 2 projects the space-fixed accelera- 
tions X 2 into body-oriented coordinates. 



sin 6 
cos <5 
0 



(7) 


Defining a body-oriented acceleration Z 2 as sensed 
possible by accelerometers 


Z 2 


f2 2 X 2 




(8) 


and applying Laplace transform to a system initially 
at rest, equation 6 can be reduced to 


F 2 = Lf 1 F, - G 2 p 2 Z 2 . (9) 

Next, the deformation of a massless but flexible 
beam shall be expressed. For small deformations, 
nonlinear terms such as Coriolis force and centrifu- 
gal force are negligible, and the difference between 
an acceleration at one beam end and the acceleration 
of the opposite beam end multiplied by the lever arm 
matrix L 2 can be used to formulate a body-fixed ac- 
celeration p 2 U 3 which is equivalent to deformation. 

p 2 U 3 = p 2 Z 3 - L 2 p 2 Z 2 . (10) 

The deformation U 3 is the deflection of the right 
beam end measured against the left beam end and thus 
can be related to the deformation of a cantilever beam. 
For certain frequency and amplitude ranges, the 
structural damping can be approximated by viscous 
damping. Applying this method, the deflection of a 
beam with damping can be derived from the deflection 
of a spring-like beam by dividing the deflection of the 
latter by 1 + pr where r is the quotient of a dashpot 
constant and a spring constant. 

P 2 u 3 = l 2 d 2 f 2 . (11) 

The matrix D 2 represents the flexibility of a can- 
tilever beam clamped at the right end and deflected at 
the left by a force F 2 . 
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2 

w 3 ,f 3 
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W 4 ,F 4 
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FIGURE 2. LUMPED MODEL OF A SATURN V SPACE VEHICLE 


The flexibility matrix D 2 is defined by 


L 


l±- 

f ea 




0 


° 

1 

s 2 

' 

-s 

kGA + 

EJ 

EJ 

s 


J_ 

' EJ 


EJ 


/ 2 


( 12 ) 


This concludes the formulation relating the forces 
and accelerations of the mass cell (stations 1 to 2.5) 
and the flexibility cell ( stations 2 to 3) . This is re- 
peated for all the mass cells (stations 2.5 to 3.5, 3. 5 
to 4. 5, and 4. 5 to 6) and all the flexibility cells (sta- 
tions 3 to 4 and 4 to 5) ( Fig. 2) . Wind forces will be 
considered by essentially adding a force W 2 to the 
right half of equation 9. 


PROGRAMING 


As mentioned previously, the space vehicle was 
divided into four lumped mass fields (G 2 , G 3 , G 4 ,G 5 ) , 
the S-IC stage, S-II stage, S-IVB stage, and the Apollo 
spacecraft. Three flexibility fields (D 2 , D 3 , D 4 ) were 
interconnecting the four masses as shown in Figure 2. 
As displayed in Figure 2 the variables are available 
at six different stations along the longitudinal axis. 
External forces are the tail force F t and the distrib- 
uted wind forces W 2 through W 6 . Internal forces are: 
the forces F 2 through F 5 . For simplicity it is assumed 
that the forces F t and W 6 will not appreciably deform 
the structure between stations 1 and 2 and between 5 
and 6. Now the equations of the lumped model can be 
arranged. 
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p 2 Z t = Lf 1 p 2 Z 2 Lever arm effect (13) 

- p z Z 2 = G 2 -1 ( Lf 1 Fj + W 2 - F 2 ) Mass effect ( 14) 

Flexibility effect 

F 2 = ( 1/p 2 ) ( 1 + PT 2 ) Df 1 (p 2 Z 2 - Lf 1 p 2 Z 3 ) ( 15) 

p 2 Zj = G. - 1 ( Lf ^ Fj i +Wj-F.) Mass effect (16) 

Flexibility effect 

F. = ( l/p 2 )(l .+ pTj) D“‘ (p 2 Z. - Lf 1 p 2 Z. +1 ) ( 17) 

p 2 Z 5 = Gs " 1 ( L 4" 1 F 4 + W 5 - F 5 ) Mass effect ( 18) 

F 5 = L 5 W 6 Lever arm effect ( 19) 

p 2 Z 6 = L 5 p 2 Z 5 Lever arm effect ( 20) 

Equations 13 through 20 describe completely the 
lumped model of Figure 2 by expressing all accelera- 
tions and internal forces dependent on external forces, 
* thus satisfying the boundary conditions of a free body. 

Since all the matrices G, D, and L (equations 4, 
5 , and 12 ) can be partitioned into similar diagonal 
matrices, the longitudinal components may be sep- 
arated from the remaining components of the vectors 


-w 2 



FIGURE 3. ANALOG COMPUTER DIAGRAM FOR STATIONS 1 
THROUGH 3, EXCLUDING POTENTIOMETERS 


w 3 



FIGURE 4. ANALOG COMPUTER DIAGRAM FOR STATIONS 2 
THROUGH 4, EXCLUDING POTENTIOMETERS 


-w 4 w 5 



FIGURE 5. ANALOG COMPUTER DIAGRAM FOR STATIONS 3 
THROUGH 6, EXCLUDING POTENTIOMETERS 


F, W, and Z. The remaining components are of spe- 
cial interest and only these components will be simu- 
lated here because they contain the information (lat- 
eral and angular accelerations) utilized by the control 
system. 

Figures 3, 4, and 5 depict the arrangement of 
operational amplifiers and integrators programed to 
model the dynamic behavior of the Saturn V space ve- 
hicle. The integrators with the triple lines deviate 
from the usual integrator setup to include the time 
constant t in their transfer function, (1 + pr)/p. 
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SIMULATION 

The model program was realized on one analog 
computer AD-2-64-PBC. The variables were scaled 
properly to keep the amplifier gains between one and 
ten 0 Figure 6 shows how the cells were connected and 
where the signal generator was attached to simulate a 
shaker on the tail of the Saturn V vehicle. Figure 7 
illustrates the arrangement of the electronic model 
together with the control system. Also indicated are 
interpolators to generate the lateral acceleration sig- 
nals between the station pairs 2-3, 3-4, and 4-5, thus 
increasing the number of stations to 9. Nine acceler- 
ation signals and one force signal were displayed on 
an oscilloscope (Ref. 10) and recorded on Brush re- 
corders. Photographs were made from the displayed 
mode shapes and traveling waves shown in Figures 8 
through 11. The transfer function of the model (Figs. 
12 and 13) was recorded by a frequency response an- 
alyzer (Industrial Measurement Corp.). 

The photographs are self-explanatory. The mode 
shapes can be seen in Figure 8; however, for modes 
higher than the third, more points have to be used. 



STATION 



FIGURE 7. SIMULATION OF SATURN V MODEL WITH ATTITUDE 

STABILIZATION LOOP CLOSED EXCITED BY UNIT 

STEP WIND 



FIGURE 6. SIMULATION OF SATURN V MODEL WITH SIGNAL 
GENERATOR SIMULATING A SHAKER 


The "0" mode is the acceleration distribution for a 
steady force at the tail (it has one node, which is the 
center of rotation) . All photographs were taken for a 
vehicle at 40 seconds flight time and the bending mode 
frequencies were 1. 08, 2. 06, and 3. 05 Hz, thus being 
4.5, 8, and 4 percent above the values of officially 
available digital data (Ref. il). 


Of special interest is Figure 9> which shows a 
disturbance traveling through the 105 m length of 
the Saturn V vehicle. The velocity is apparently 
210 ms” 1 and agrees well with the firstbending mode 
resonance considerations. About half the time is re- 
quired for a disturbance starting at the center of pres- 
sure (Fig. 10) . 

The photographs shown in Figures 10 and 11 
are the traveling waves excited by a unit step wind 
force at the center of pressure while the attitude of 
the whole vehicle model is stabilized by a control loop 
(Fig. 7) . Here it is obvious that the disturbance rings 
out faster, indicating the damping of the first bending 
mode by the control system. The photograph at 10 
seconds after the excitation at the center of pressure 
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2ND MODE ACCELERATION EXCITED BY SINUSOIDAL FORCE OF 
2 x 10 5 N pp AT THE TAIL . 




3RD MODE ACCELERATION EXCITED BY SINUSOIDAL FORCE OF 
4 x 10 5 Npp AT THE TAIL. 



"0” MODE ACCELERATION EXCITED BY UNIT STEP FORCE OF 1ST MODE ACCELERATION EXCITED BY SINUSOIDAL FORCE OF 

2.2 x 10 6 N AT THE TAIL PUSHING TO RIGHT AND TAKEN 70 s 2 x 10 5 N pp AT THE TAIL OF SA-V MODEL. 

AFTER STEP STARTED. 

FIGURE 8. MODE SHAPES OF ACCELERATIONS PERPENDICULAR TO AXIS OF SATURN V MODEL. UNITS: VERTICAL, 15 m/div; 
HORIZONTAL, 5 m/s 2 /div. 
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clearly shows how the vehicle accelerates to the right, 
pushed by the wind and the attitude stabilizing compo- 
nent of the thrust. At 10 seconds the rotation of the 
vehicle has completely stopped and only a lateral 
drifting prevails. 

The recordings in Figures 12 and 13 show the 
transfer function of the model for the signal ratio of 
nose acceleration to the tail excitation force. Note 
the four resonances with increasing damping towards 
the higher frequencies. The critical damping factors 
were approximately 1, 2, 3, and 6 percent. 


CONCLUSIONS 

An electronic model of a flexible space vehicle 
was successfully used for simulation of traveling wave 
effects, bending modes, rigid body motions, and 
transfer functions. It promises to be a good tool to 
gain better insight in the dynamic behavior of flexible 
space vehicles. It is a very efficient method to simu- 
late with real hardware, especially if it is necessary 
to vary structural vehicle parameters. Available an- 
alog computers can be used, thus demonstrating the 
great versatility of the model. 
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F. CONSTRUCTION OF ELASTIC BODY MODELS FROM NOMENCLATURE 

SATURN DYNAMIC TEST RESULTS by John H. George 

and Mauriee E. Singley & 

x angular deflection at station x, deg 


F 


, N65 U7 8Z 

* -> 

~ ~ m 

This article presents a mathematical 
obtaining mathematical models for dynamic systems. 
This procedure was used to obtain mathematical models 
of the elastic body dynamics for Saturn cluster mode 
vehicles. 


? Asr- 

procedure 


for 



M 

u 


N 

u 


Theoretical procedures were not available for ade- 
quately describing the elastic body dynamics. For this 
reason, it was necessary to utilize this # procedure for 
obtaining mathematical models directly from test data. 
Examples were taken from flight vehicle dynamic testing . 




Cl) 

U 


force input, N 

angular deflection at station x caused 
by the vfi 1 bending component, deg 

gain term obtained by partial fraction 
expansion of the vfi 1 bending component, 
deg/ s 2 -N 

lead term obtained by the partial fraction 
expansion of the u^ bending component, 
deg/ s-N 

damping ratio associated with the u** 1 
bending component 

circular frequency of vP 1 bending com- 
ponent, 1/s 


<p. 

xa 


amplitude of angular velocity at station 
x, deg/s 


INTRODUCTION 

During the early days of Saturn research and de- 
velopment, it was necessary to supplement the 
theoretical elastic body data by a dynamic test program. 
The dynamic test results were to determine local ef- 
fects, non-planar coupling, and structural damping 
that were not given or explained by theoretical con- 
siderations. A method was needed that would convert 
dynamic test data into a mathematical model of the 
flexible vehicle dynamics, including the complicated 
multiple tank behavior of the Saturn vehicle. 


<p. phase of angular velocity at station x 

ip referenced to the oscillator, deg 

F^ amplitude of load cell force, N 

F phase of load cell referenced to the oscil- 

^ lator, deg 

MATHEMATICAL MODEL CONCEPT 


The procedure described in this article was de- 
vised to convert empirical frequency response data 
obtained from the Saturn dynamic test program into 
a mathematical model describing elastic body behavior. 
The technique was derived in general and is valid for 
any station of the vehicle, with the accuracy depending 
only upon the input data. The resulting mathematical 
model includes all local effects, coupling, damping, 
and the effects of multiple tank behavior. The mathe- 
matical models obtained by this procedure were used 
in the design, synthesis, analysis, and final checkout 
of control systems of the Saturn SA-5, SA-6, and SA-7 
flight vehicles. 


In general, any transfer function may be written 
as a ratio of polynomials, 


<P.( s) 

F(s) 


a + a<s + . . . 
o i 


b + b s + . . . 
o i 


+ a s 
n 


+ b, 


.2m 


2m" 


n < 2m 


(i) 


where the exponents of the polynomials, n and m, are 
determined by the order and the coupling of the physi- 
cal system. Expanding equation 1 by a standard par- 
tial fraction expansion and assuming only complex 
zeros in the denominator yield the normal mode form 
of the transfer function as 
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c p As ) m M + N s 

Y V ' _ y U U 

F(s) " s 2 + 2£ w s + w 2 ' 
' ’ u=i u u u 


( 2 ) 


The second order ordinary differential equation of the 
elastic body behavior caused by the first m bending 
modes may be obtained from equation 2 by taking the 
inverse Laplace transformation, 

w + 2 1 c 0(p. + (p . = M F + N F 
Y iu u uMu u Y iu u u 

(3) 

u = 1, 2, . . . , m. 

The coefficients of equations 2 and 3 are constructed 
from the dynamic test frequency response data using 
the following method. 

The input force, measured by a load ceil, is re- 
corded as an amplitude F a and a phase F p . The output 
(j?i of a rate gyro at station i has an amplitude <^ a and 

<p ip . The amplitude function is calculated from the 
dynamic test data as 



(4) 


where co = 2?rf and f is the frequency of the forcing 
function. The phase relationship is obtained from the 
expression 



w - (F + 90°). 
ip P 


(5) 



FIGURE 1. DYNAMIC TEST SETUP 


This numerical transfer function is used as an input 
to Levy's curve fit routine (Refs. 1 and 2) which has as 
an output a ratio of polynomials, or a transfer function, 
in the form of equation 1. This is expanded by a par- 
tial fraction expansion method to produce equation 2. 

SATURN DYNAMIC TEST 

A typical dynamic test vehicle with test instru- 
mentation is shown in Figure 1. The vehicle j.s sus- 
pended in the dynamic test tower by a network of cables 
attached to the fins at one end and to the test tower 
through a spring package at the other. The number of 
springs in each package may be varied for different 
suspension effects for each simulated time of flight. 
Currently two different spring configurations are being 
used for each simulated time of flight. 


The input force is applied at the engine swivel 
point (Fig. 1) and is sinusoidal with a frequency that 
is controlled by a crystal oscillator (Ref. 3). 

The flow of data from the dynamic test to a form 
that is usable by the control system design and synthe- 
sis group is shown in Figure 2. The output of the load 
cell, the oscillator, and each of the sensors on the 
dynamic test vehicle is recorded on digital tape and 
strip charts. A primary reduction is performed on 
the digital tape and the printout is compared with the 
strip charts for determining unusual wave shapes, 
proper instrument polarity, and correct calibration. 
The results of this comparison are reported to the 
operations group, who make the appropriate changes. 

A secondary reduction is performed on the digital 


56 



ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 



FIGURE 2. DYNAMIC TEST DATA FLOW DIAGRAM 


tape yielding; (a) frequency response data of the form 
of equations 4 and 5, and (b) computer tape of the 
frequency response of certain selected instruments. 
The frequency response data from (a) are documented. 
The tape from (b) is used as an input to Levy's curve 
fit routine and expanded by the partial fraction method, 
thereby yielding a mathematical model matching the 
amplitude and phase of the elastic vehicle response in 
the sense of least squares (Fig. 3). 

A linear extrapolation is performed by considering 
the bending model coefficients versus the number of 
springs in the two different suspension systems for 
each simulated time of flight. The bending model co- 
efficients found by extrapolation characterize a sus- 
pensionless vehicle which approximates a free-free 
system. 



FIGURE 3. MATHEMATICAL MODEL 
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APPLICATION TO SATURN DYNAMIC TEST 

Tables 1,2, and 3 and Figures 4 through 9 typify 
results obtained during dynamic testing of the Saturn 
dynamic test vehicle. The coefficients (Tables 1 , 2, and 
3) were computed using the procedure outlined in the 
preceding section. They describe the first four modes, 
as measured from the exciting force at the tail to a 


rate gyro mounted in the instrument compartment. 
Table 1 shows the computed coefficients for the test 
using 8 supporting cables and 12 springs. Table 2 
shows the computed coefficients for the test case using 
8 supporting cables and 8 springs. Table 3 shows 
the coefficients obtained by linear extrapolation to the , 
free-free case. Figures 4 through 7 are comparisons 
of the frequency response functions computed from the 


TABLE 1. BENDING MODEL COEFFICIENTS USED IN FIGURES 4 AND 5. 


f 

i 

M 

N 

1.73 

0. 0104 

8. 69 

-0. 0459 

3.27 

0. 0080 

1. 75 

0. 0142 

4. 64 

0.0065 

-7. 37 

0.0063 

5. 31 

0. 0045 

2. 39 

0. 0396 


.90x10* 5 



FREQUENCY (Hz) 



FREQUENCY (Hz) 


FIGURE 4. AMPLITUDE FUNCTION COMPARISON, HARD SUSPENSION FIGURE 5. PHASE FUNCTION COMPARISON. HARD SUSPENSION 
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models of Tables 1 and 2 with respect to the original 
measured data. Figures 8 and 9 are frequency re- 
sponse functions of the model given in Table 3. 

Extrapolation from a vehicle with soft suspension 
to a suspensionless vehicle resulted in an average 


modal frequency shift of less than 1 percent. The re- 
sulting change of the gain associated with each modal 
frequency averaged less than 30 percent. The example 
cited from the Saturn SA-6 dynamic test had an average 
frequency shift of. 74 percent and gain change of 15. 1 
percent. 


TABLE 2. BENDING MODEL COEFFICIENTS USED IN FIGURES 6 AND 7. 


f 


M 

N 

1.74 

0.0101 

9. 13 

-0. 0351 

3.28 

0. 0070 

1. 68 

0.0073 

4.65 

0. 0063 

-7. 81 

0. 0386 

5.30 

0. 0045 

2.21 

0. 0199 


.90x10- 5 



FREQUENCY (Hz) 

FIGURE 6. AMPLITUDE FUNCTION COMPARISON, SOFT SUSPENSION 



FREQUENCY (Hz) 


FIGURE 7. PHASE FUNCTION COMPARISON, SOFT SUSPENSION 
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G. SIGNAL CONDITIONING CIRCUIT FOR AN OSCILLO- 
SCOPE DISPLAY OF BENDING by Charles P. Elms 

_ N65-23788 

t 

The problem of displaying the bending, 
space vehicle or a structural beam is treated by employ- 
ing a switching circuit that samples the outputs of ac- 
celerometers fixed to the frame of a vehicle and displays 
them on the screen of an oscilloscope. The sampled data 
are displayed as a vertical column of points represent- 
ing the sensors. Vertical spacing is proportionate to the 
spacing of the sensors on the vehicle. The bending mo- 
tion is represented by horizontal motion of the points 
proportionate to the bending motion sensed. 



action of a 


thevehicle. Forthiscase, ten such dots are generated; 
nine are accelerometer outputs and one displays the 
driving force. This is achieved by time sharing the 
vertical and horizontal inputs of the oscilloscope (Ref* 
1) . The time sharing is done at such a rate that all 
the dots seem to appear on the screen at the same time , 
and time shift between points is negligible compared to 
the shortest expected period of bending oscillation to 
be displayed. 

NOMENCLATURE 


A amplitude of sine wave input 

B(t) blanking signal 


A formula for calculating the minimum required 
sampling rate for an input of chosen maximum frequency 
is derived, based on the amount of maximum tolerable 
horizontal display error. 


, The circuit can be used with any single beam os- 
cilloscope having external intensity modulation capa- 
bility. It is of simple design and is useful in applications 
other than bending where more than one variable is to 
be displayed on one instrument. 


The circuit has been built and used to display the 
bending motions of the Saturn V vehicle simulated by an 
electronic elastic beam model . K v & 

W 



d 


E n (tl 

f 

c 

f. 

1 


H 

n 


H 

n 


n 


N 

P 


maximum horizontal deflection error be- 
tween points n = U and n = N^ where t = 0 

n th switching pulse 

clock rate 

frequency of sine wave input 


horizontal deflection at location n for sine 
wave input 

H when t = 0 
n 

(1, ... 10) index of ring counter stages 

number of points displayed (also number 
of ring counter stages) 


INTRODUCTION 


t time 


V 

In the simulation of the bending or actual dynamic 
testing of thevehicle, it is advantageous if one can ac- At 
tually see the movements taking place. 

• 6 

This article describes a technique for the nearly 
simultaneous display of accelerometer outputs as a 
function of their vertical location on the vehicle. 


vertical dot position bias 

time shift between points n & and n &+1 

acceleration 


SYSTEM DESCRIPTION 


By generatinga dot corresponding to each location 
of an accelerometer and having its movement horizon- 
tally displayedon the screen proportional to the accel- 
erometer outputs, one can actually see the bending of 


The system contains a clock generator (Ref. 2) , 
a ten-stage ring counter (Ref. 3)_with buffers, 20 gate 
circuits, and a blanking pulse circuit. 
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FIGURE 1. VEHICLE AT REST 


FIGURE 3. SECOND BENDING MODE 



FIGURE 2. FIRST BENDING MODE 



FIGURE 4. ENVELOPE OF FIRST BENDING MODE 


Figure i depicts the display as seen on the oscil- 
loscope when the vehicle is at rest. Figure 2 shows 
the vehicle in the first bending mode, and Figure 3 
shows the second bending mode. Figure 4 is a time 
exposure of the first bending mode, thus showing its 
envelope. 

A block diagram of the circuit is presented in 
Figure 5. The clock is stepping a ten^stage ring 
counter. Each stage of the ring counter in succession 
turns on a set of gates. The horizontal gates (the 
a-gates) accept the acceleration signals and produce 
time sharing of the horizontal deflection of the oscil- 
loscope. The vertical gates (the b-gates) set up the 


vertical bias levels corresponding to the location of 
the accelerometers along the vehicle axis and effect 
time sharing of the vertical deflection of the oscillo- 
scope. Potentiometers for zero setting, amplitude 
scaling, and vertical bias level are integral parts of 
the gate circuits. 

The sequential switching from gate to gate must 
be done fast enough that there is virtually no time de- 
lay between samples. When this is done, an undesir- 
able spike occurs. The spike is eliminated by a blank- 
ing pulse generated by differentiating the clock gene- 
rator output which modulates the cathode of t Tektronix. 
RM 561 A Oscilloscope (Ref. 4). 
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When 27r — n is small, 
c 


H = A 2tt — n 
n max f 

c 


(5) 


To find the error caused by time shift between • 
points n & and n^ , 


(H - H„ ) = A 2tt — 
n b n a max f 


(n b- n a } • (6) 


If n^ = Np and n^ = 1 , then equation 6 yields the 
horizontal display error d . 


d = (H - H t ) = A 2 tt — (N - i) (7) 

v N max f p 


FIGURE 5.* SYSTEM BLOCK DIAGRAM 


Selection of the sampling rate (clock rate) is based 
on two considerations. First, the clock rate must be 
so fast that the eye cannot detect the sampling. Second, 
the rate must be fast enough to make the error between 
samples, caused by time shift, negligible. It is shown 
that the first requirement is fulfilled by the necessarily 
high rate specified by the second consideration. 


From equation 7, a clock rate f may be computed 
to satisfy a maximum input frequency F max to mini- 
mize the maximum horizontal display error resulting 
from time shift. 


For example, let A =4 , N 

max cm p 


10. f. 


i max 

10 Hz, and d = .2 cm. Then from equation 7, f = 
22,608 Hz. ° 


Consider the horizontal deflection for sine inputs 
of frequency f. max at locations n as 


max 


c 


H = A sin 27r f. (t + At n) 


imax 


( 1 ) 


The worst case condition will be when t = 0, the steep- 
est slope of the sine wave. 


H = A sin 
n max 


2ttF At 
imax 


but 

A t = 1/f 


n 


( 2 ) 


(3) 


Then equation 2 becomes 


f. 

i 


H = A sin 
n max 


2tt y n 
c 


(4) 


CONCLUSIONS 


The circuit has been installed and is functioning 
very satisfactorily. Some results of its use can be 
found in Reference 5. 

The main advantages are its compatibility with any 
standard laboratory oscilloscope having external inten- 
sity modulation capability , the simplicity of the circuits 
involved, and its usefulness in dynamically displaying 
the bending of a space vehicle. 

The discontinuous line (or array of dots) is the 
main disadvantage. It wouldhe more desirable to have 
a continuous line to display the bending. Consideration 
of this is presently being given and development work 
is being done to implement it. 


64 






ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


REFERENCES 

1. Elms, Charles P. ; Signal Conditioning Circuit for Displaying Pitch, Roll, 
and Yaw on an Oscilloscope. Astr ionics Research and Development Report 
No. 2, NASA TM X-53044, May 1, 1964. 

2. Hurley, Richard B. ; Asymmetrical Free-Running Multivibrator , Junction 
Transistor Electronics , p. 421, John Wiley and Sons , Inc., New York, N.Y. 
1960. 

3. Marcus, Ira R.; Design Procedure for a Transistor Ring Counter. TR-1154, 
United States Army Material Command, Harry Diamond Laboratories, Wash- 
ington 25, D. C., October 15, 1963. 

4. Instruction Manual; Tektronic, Inc., P. O. Box 500, Beaverton, Oregon. 

5. von Pragenau, George L. ; Real Time Analog Computer Simulation of the 
Elastic Space Vehicle Involving Beam Cell Transmission Matrices. Pub- 
lished in this report. 


65 


ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


H. SIMULATION OF FLEXIBLE BODY DYNAMICS BE- 
HAVIOB by Glen D. Bitter 

$e>s^ N 65 23 789 

A method for computing flexible body natural fre- 
quencies, mode shapes and slopes, and generalized mass 
during the real time simulation of large launch vehicle 
flight dynamics is discussed. Derivation of the equations 
and evaluation of the proposed method through compari- 
son of results with existing theoretical data computed by 
conventional techniques are described. 

Although the method was derived for the real time 
simulation application, it could be used in the routine 
computation of elastic body data. The method is useful 
because it avoids the inherent cumulative errors and the 
difficulties of high ordered determinant manipulation of 
more conventional approaches. I \ 


INTRODUCTION 


In the treatment of large flexible launch vehicle 
control problems, a recurring question is asked by the 
control system engineer: Could an analog computer be 
adapted to the solution of the flexible vehicle dynamics 
in such a way that it would yield the significant modes 
quickly enough for a real time simulation? The meth- 
od discussed here answers this question affirmatively. 
The approach used extends the result of Mori (Ref. 1) 
by providing a method, based on well known principles 
of ordinary differential equations, for quickly deter- 
mining the proper end conditions. This method could 
be extended to a ’’branched’' beam structure such as 
the Saturn I and I-B multiple beam structure. This 
method could also be used, for example, in the prob- 
lem discussed by Fettis (Ref. 2) of a two beam struc- 
ture. 

In obtaining equations of lateral vehicle motion, 
it has been established (Refs. 3 and 4) that structures 
the size of a space vehicle must be assumed nonuni- 
form, and the effects of shear and rotary inertia must 


be included. These effects are included in a proce- 
dure that should be easily adapted to an analog simu- 
lation of the flexible vehicle dynamics. The procedure 
is evaluated by comparison with data generated from 
an accepted technique in use at NASA’s Marshall Space 
Flight Center. 

Although the method was derived for the real time 
simulation application, it could be used in the routine 
computation of elastic body data. The method is use- 
ful because it avoids the inherent cumulative errors 
and difficulties of high ordered determinant manipu- 
lation of more conventional approaches. 


NOMENCLATURE 


A 

cross-sectional area 

m 

Ci» C2 

constants 


El 

flexural rigidity 

N m 2 

M 

generalized mass 

kg 

I 

cross-sectional moment of inertia 

m 

kAG 

shear rigidity 

N 

L 

length of launch vehicle 

m 

M 

bending moment 

N m 

Q 

shear force 

N 

t 

time 

s 

u 

potential energy 

N m 

X 

measure of length along vehicle 
axis (station) 

m 

y 

lateral deflection 

m 

7 

mass per unit volume 

kg/m s 

* 

slope of deflection curve caused by 
bending, neglecting shear 

rad 

CO 

angular frequency 

rad/s 

p 

slope of deflection curve caused 
by shear 

rad 
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DESCRIPTION OF FLEXIBLE VEHICLE EQUATIONS 


To obtain the equations of motion of a flexible 
space vehicle, consider a finite element of the vehicle 
( Fig. 1). The following assumptions are made*. 

1. The vehicle has an axis of symmetry. 

2. The vehicle is not deflecting beyond its elastic 
limit. 

3 . All vibrations are in the plane under consider- 
ation. 

4. At natural modes, the motion is simple and 
harmonic. 

It is assumed that a space launch vehicle can be con- 
sidered as an infinite set of uniform elements. The 
equations of motion become partial differential equa- 
tions by a limit process and are then valid for a non- 
uniform vehicle. 

♦ 

The equations of lateral motion for a uniform beam 
section as shown in Figure 1 can now be derived, as- 
„ suming variable moments and shear forces acting on 
both ends of this section. These equations in partial 
differential equation form are given (Ref. 5). 


3*(x,t) 

M (x, t) 

(i) 

dx 

EI(x) 

9Y(x, t) 
3x 

+ ^( X t) 
kAG(x) ( ’ ' 

(2) 

3Q(x.t) 

3x 

yam d > 

(3) 

3M(x, t) 
dx 

dx + Q(x) dx - Iy(x) 

dx. (4) 


The displacement Y, the slope neglecting shear effects 
moment M, and shear Q must be a function of time 
and distance along the vehicle axis- to satisfy the set of 
partial differential equations, 1, 2, 3, and 4. In gen- 
eral, any vibration of a system can be resolved into 
the natural modes of vibration. When a system is vi- 
brating at one of its natural modes, the system expe- 
riences simple harmonic motion; therefore, particular 
solutions of the equations can be found (Ref. 5). As- 
suming that the vehicle under consideration is vibrat- 
ing at a natural mode, solutions of the partial differ- 
ential equations are given as 


ELEMENT 

DISPLACEMENTS 



Y ( x, t) = Y(x) sin (cot + <p) (5) 

*(x,t) = *(x) sin (cot + (p) (6) 

M(x,t)= M(x) sin (cot + cp) (7) 

Q(x, t) = Q(x) sin (cot + cp) • (8) 


Differentiating equations 5 and 6 twice with respect to 
time, 


, ( * = -Y(x) co 2 sin (cot + cp) 

and 


d 2 ^(x,t) 

dt 2 


-'I'(x) co 2 sin (cot + cp ) . 


(9) 


( 10 ) 


Substituting equations 5 through 10 into equations 1 
through 4 and factoring out the sin (cot + cp) terms., 
the ordinary differential equations are obtained. 


d^ (x) 
dx 

dy(x) 

dx 

dQ(x) 

dx 

dM(x) 

dx 


M(x) 

EI(x) 


Q(x) 


kAG(x) 

= -yA(x) co 2 Y(x) 


+ ^(x) 


= Q(x) + Iy(x) co 2 *(x). 


( 11 ) 

( 12 ) 

(13) 

(14) 
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These equations will now yield a solution for any 
angular frequency cu. However, only the solutions that 
satisfy certain boundary conditions are of interest. 


normal mode frequencies, and this is treated by the 
usual method of linear equations. The equations have 
a solution if and only if 


The co values for which the particular boundary 
conditions are satisfied are the natural frequencies of 
the flexible vehicle described by equations il through 
14. At these frequencies the natural mode shapes can 
be obtained from the solution of these equations. This 
process will now be discussed for a free-free flight 
vehicle. 

For the free-free condition, the moment and 
shear at both ends are zero. A frequency go will now 
be determined that will satisfy the boundary conditions. 
For a given co, two solutions will be determined with 
the following initial conditions: 

Y(0) ¥( 0) M( 0) Q ( 0) 

Solution denoted 

by subscript 11 0 0 0 

Solution denoted 

by subscript 2 0 10 0 


Mj(L, co) M 2 (L,co) 


F(c o) 


Qi ( L, a;) 


Q? ( L, co) 


= 0 , ( 21 ) 


or, expanding the determinant, if 


F (co) = Mj(L, co) Q 2 (L,co) - Qj(L,co) M 2 (L,co) = 0. 

( 22 ) 


This defines a continuous function of to. To obtain the 
normal mode frequencies, it is necessary to determine 
the zeros of this equation. Several methods could be 
used: a three point extrapolation method or, as used 
in this article, Newton's method with the derivative 
F T (c o) approximated by the first difference quotient. 
Equation 22 possesses all of the infinite frequencies 
associated with this nonuniform beam problem. After 
a particular frequency co t is determined that satisfies T 
equation 22, it is still necessary to determine the 
mode shapes. At co = coj, equation 19 yields 


The two sets of solutions are independent solutions _ M^L,^) 

of thegiven differential equations, and thegeneral so- C 1 M 2 (L,co t ) 

lution can be written: 

and by choosing 


(23) 


Y(x) = C^Y^co) + C 2 Y 2 (x,co) (15) 

*(x) = C^^x, co) + C 2 * 2 (x,co) (16) 

M(x) = C^^co) + C 2 M 2 (x,co) (17) 

Q(x) = CjQ^x, co) + C 2 Q 2 (x,co). (18) 


This solution has enforced the free-free boundary 
condition at the end x = 0, by the selection of the ini- 
tial conditions. The boundary conditions at this end 
are 


C 1= 1, 

the mode shapes, slopes, moment, and shear are eas- 
ily determined from equations 15 through 18 at the 
frequency co - coj. 


GENERALIZED MASS 


M(0,t) = M(0) sin (cot + cp) - 0. 

Q(0,t) = Q(0) sin (cut + <p) = 0. 

To satisfy the remaining boundary conditions at 
to other end of the vehicle (x = L) , the following 
equations, from equations 17 and 18, must be satis- 
fied. 

CjM^L, co) + C 2 M 2 (L, co) = 0 (19) 

CjQj (L, co) + C 2 Q 2 (L, co) = 0. (20) 

The solution to these equations now yields the 


While the natural frequencies and mode shapes 
are important to control system design, the general- 
ized mass, another variable (Ref. 3) needed in con- 
trol analyses, can be found from the modal solution of 
the equations of motion. 

The differential expression for potential energy 
of a finite section of the vehicle is 

dU (x) = — M(x) • dx + |Q(x)/3(x) dx. 


68 



ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


Integrating this equation yields 


' U=| J M(x) • dx + | J Q(x) -/3(x) dx. 

o o 

3^ 

Substituting —7— from equation ii, and since 
dx 

P ( x ) = (from Fig. 1 and equation 2), 

kACj( x) 

the potential energy equation becomes 


U 



/ M 2 (x) 
\ EI(x) 


\ 

kAG(x) J 


dx 


The maximum kinetic energy of a conservative 
system at a particular normal mode frequency w is 
defined by 

T = |moj 2 

where M‘is the generalized mass of the normal mode 
-frequency o>. 

Under the assumption that at a normal mode the 
\maximum potential energy of a system is equal to the 
maximum kinetic energy, 


REFERENCE Mt9 
CALCULATED 


GENERALIZED MASS 
CALCULATED 
REFERENCE 


1970 


1900 + 


1800 4- 


1700 + , 


/ + 1600 


/ 

/* 

/ 4- 1500 


\ l 


1300 


\ 


30,790 kg 
33,364 kg 


1200 FREQUENCY 

CALCULATED 

REFERENCE 


1100 

1066 




9.99 Hz 
9.98 HZ 


— Mco 2 = - 
2 2 


J 

o 


EI(x) 


1U1V J ^ J 


The generalized mass M is therefore 


m = ^ r + 

CO 2 J lEIfxl 
o 


Q 2 (X) 
EI(x) ' kAG(x) 


dx. 


(24) 


M(x) and Q(x) are obtained from the preceding sec- 
tions. Equation 24 can then be solved to yield the gen- 
eralized mass at any particular natural mode. 


FIGURE 2. S-IV WITH JUPITER NOSE, FIRST MODE 

of 9. 99 Hz, compared to 9. 98 Hz obtained from bend- 
ing data.* The mode shape obtained was also close 
to the mode shape obtained from the bending data ( Fig. 
2) . The second mode frequency was calculated to be 
21.56 Hz, compared to 24.2 Hz from the reference 
data. The second mode shape is also close, as can 
be seen from Figure 3. The first mode generalized 
mass for the S-IV stage with a Jupiter nose cone pay- 
load calculated by this method is 30,790kg, compared 
to 33,364 kg of the reference data. 


CONCLUSIONS 


The results of this method compare closely with 
those of other theoretical methods. The solution of 
equations by this method, considering the S-IV stage 
with a Jupiter nose cone, yields a first mode frequency 


It is concluded that the method described will be 
an effective normal mode analysis technique for non- 
uniform space launch vehicles because of its accuracy 
and the ease of application to relatively complicated 
problems. One advantage of this method is that the 
higher modes are not dependent on the solution of the 
lower modes; therefore error is not accumulated from 


* Unpublished Aeroballistic theoretical SA-5 bending 
data. 
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REFERENCE •••• 
CALCULATED 


2070 + 


2000 + 


1900 + 


1800 


4 - 1700 


+ 1600 


4 1500 


1400 4 - 


1300 -K 


GENERALIZED MASS 
CALCULATED 5,331 kg 

REFERENCE 6,444 kg 



FREQUENCY 

CALCULATED 

REFERENCE 


« 

/ 

-- 1100 
1066 


21.36 Hz 
24.2 Hz 


FIGURE 3. S-IV WITH JUPITER NOSE, SECOND MODE 


mode to mode as is the case with matrix methods. A 
disadvantage in using this procedure on a digital com- 
puter is the computation time necessary to numerically 
integrate the system of differential equations. 


While the method is completely amenable to digi- 
tal computation, use of an analog computer to inte- 
grate the system of ordinary differential equations ii 
through 14 should result in a substantial saving in 
overall computation time. Use of analog computation, 
with machine time t representing the variable x, would 
also permit properties of the vehicle to be stored on 
either function generators or a cam device, thereby 
permitting a complete characterization of continuous 
vehicle properties. An analog block diagram for so- 
lution of the system of ordinary differential equations 
11 through 14 is shown in Figure 4. Figure 5 is a flow 
diagram summarizing the overall technique. 



FIGURE 4. ANALOG SCHEMATIC DIAGRAM 



FIGURE 5. SOLUTION FLOW CHART 

One possible area of further investigation is the 
manner in which the zeros of F (co) are obtained. New- 
ton's method, replacing the derivative F'(cu) by the 
difference quotient, gave entirely satisfactory results 
in the cases considered. However, it is quite possible 
that other procedures would be easier to implement 
and faster on an analog computer simulation. . This 
possibility is also currently under investigation and 
will be discussed in a future article. 
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INSTRUMENTATION SYSTEMS 



A. A PROPELLANT SLOSH MEASURING SYSTEM FOR 
SPACE VEHICLES by B. G. Bynum and John F. Hamlet 




N65-23790 


A knowledge of the slosh characteristics of liquid 
propellants in space vehicles is necessary to determine 
control, structural, and pressurization parameters. This 
article describes the theory and configuration of a ca- 
pacitance system that may be used to measure the first 
mode of propel I ant slosh . Two pairs of parallel flat-plate 
capacitance probes extending the length of the propel- 
lant tank are connected into two bridge circuits. Differ- 
ence in liquid height on diametrically opposed probes 
causes a difference in capacitance between the two 
probes. This difference unbalances the bridge circuit 
and causes a voltage output from the electronics package 
proportional to the height of slosh. The electronic theory 
of the system ispresented with an explanation of the cir- 
cuit design and operation. The design of the probes is 
described, and supporting theory is presented. Circuit 
diagrams and illustrations of the electronics network, 
probe construction, and vehicle installation are includ- 
ed. Tests proved that the system measured slosh accu- 
rately within 1.5 cm. y\ 



INTRODUCTION 

A knowledge of the slosh characteristics of liquid 
propellants in space vehicles is necessary to deter- 
mine control, structural, and pressurization param- 
eters. This report describes a plan whereby capac- 
itance sensors may be used to measure the first 
mode of propellant slosh. 

Slosh is normally measured by differential pres- 
sure transducers or computed from liquid level meas- 
urement data. The technique presented here utilizes 
four liquid level measuring sensors mounted on the 
inside of a circular tank at 90-degree intervals and 


extending the length of the tank. Differential liquid 
level between pairs of opposing sensors is computed 
and transmitted to earth. 

The slosh sensors are flat plate capacitors in 
which the propellant and ullage gas are the dielec- 
trics. The capacitance of a sensor is a function of 
propellant height and dielectric constant of the media 
in which the sensor is submerged. Slosh is measured 
by measuring the difference in capacitance of oppos- 
ing probes. Figure 1 shows the location of the slosh 
sensors in a propellant tank. The advantages of this 
scheme are (1) simplification of data reduction on 
the ground, (2) conservation of telemetry bandwidth, 
and (3) reduced circuitry. Two prototype sensors, 
each 10. 16 meters inlength, and a breadboard circuit 
were built and tested. 

TOP VIEW 



FIGURE 1. TANK INSTALLATION OF RP-I 
SLOSH MEASURING PROBE 
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NOMENCLATURE 
C sensor capacitance 

C g stray capacitance of sensor 

C a active capacitance of sensor 

h height of liquid 

H length of sensor 

e dielectric constant of liquid 

L 

e G dielectric constant of gas 

angular velocity of excitation vector 
E Q output after amplification 

E. bridge output voltage 

A amplifier gain 

R i 

-E 2 E 1 bridge excitation voltage 

k constant 

Cf full capacitance of a sensor 

C e empty capacitance of a sensor 

a slope of liquid surface in h-x plane 

x z h tank coordinates 

b slope of liquid surface in h-z plane 

ELECTRICAL CHARACTERISTICS 

The capacitance of a sensor is given by equation 
1; the sensor is assumed to be linear to the required 
accuracy oyer its entire length. 

c = C s + Ca + e G 

= C s +C a [e G + (e L -e G ) 


Figure 2 shows a bridge circuit which may be 
used to translate the capacitance of two sensors into 
differential liquid height. Cj and C 2 are the capaci- 
tances of two diametrically opposed sensors. C 3 , C 4 , 
and C 5 are cable capacitances and sensor capacitances 
to ground. E t and E 2 are bridge voltages for capaci- 
tances Cj and C 2 , respectively, of the two sensors. 
Ri is a variable resistor used to balance the bridge by 
adjusting the magnitude of E^ 

If Ri is chosen small enough that no significant 
phase difference exists between Ei and E 2 , the follow- 
ing equation defines Ei, the voltage resulting from 
unbalance of the bridge. 

^ _ Ei jttfCj - E^jcoC? 

E i “ 1 * 

j w (Ci + c 2 + c 3 ) + — 

If the conductance “ is much greater than the 

susceptance cv{C i + C 2 + C 3 ), the above equation may 
be simplified to 


E 0 , the output of the circuit representing a dif- 
ference in capacitance between the two probes, is ob- 
tained by amplification of E^ The gain of the ampli- 
fier is represented by A. 

Eg — AE| = AR 2 (EijtoCj “ E 2 jcoC 2 ) 

= K (EjCj - E 2 C 2 ) (2) 

where K = Aja>R 2 . 



FIGURE 2. SCHEMATIC OF SLOSH SENSOR ELECTRONICS 


bridge zero adjust potentiometer 


Ei - R 2 (EjjuJCi E 2 jcoC 2 ). 


H-h \ 
H ) 

— 1 


( 1 ) 
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Substituting equation 1, which represents the ca- 
pacitance of a probe, for C\ and C 2 in equation 2 gives 
the following. 

E 0 = K (c Sj + [e G + (e L - e Q ) J 

" E 2 + C \ [ e G + (e L " e G ) if ]) 

h h 

= K [E 1 C a ^(e L - e G ) - E 2 C a ^ -^~(e L - e G ) 

+ Ei ( C s + C a e G ) E 2 ( C s + C a ® G )^J • , 

For hj = h 2 = 0, 

E 0 = K ^Ej ( Cg + C a 6 q) - E 2 ( Cg + C a e G )^ . 

If Ej is adjusted by Rj so that E 0 = 0 when hq = 
h 2 = 0, 

Ei (C s + C a e G ) = E 2 (C s + C a e G ) (4) 

l 1 2 2 

and equation 3 becomes 
K 

E 0 = ^ E l^a^ h^ - E 2 Ca b 2 ) ( e L ~ e G^ ’ ( 5) 

The stray capacitance of a sensor is proportional 
to the sensor active capacitance; the proportionality 
constant is identical for all sensors of identical con- 
struction. This relationship is expressed by the fol- 
lowing equations. 


C S1 = Bc al 

(6) 

C s ~ BC a « 
2 2 

(7) 

Substituting equations 6 and 7 
the following. 

into equation 4 gives 

Ej BC a ^ + ^"a 2 e G ^ a 2 


E 2 BC a ^ + C a e G C aj 


E i Ca = E 2 C a . 
1 2 

(8) 


Substituting equation 8 into equation 5 gives a final 


expression of circuit output. The output E 0 may 
be filtered and biased to obtain the proper telemetry 
voltage range. 

E 0 = “ Ca^Ej (h* - h 2 ) (e L -e G ). (9) 

It can be seen from equation 9 that E 0 will be pro- 
portional to h t - h 2 if C a (e^ - e G ) is known. 

The effect of the variables e^ and e G may be 
eliminated by subtracting the empty capacitance C e 
(C e - C s + C a e G ) f rom the full capacitance C£ 

< c f = C s + C a e I> Thus 


c f - c e = C a (e L - e G ). 

which is the unknown in equation 9. The only variable 
in equation 9 is, therefore, (hj -h 2 ). 

The full capacitance of a sensor may be measured 
by using the signal conditioning circuit. If a known 
fixed capacitor is switched into the circuit in place of 
one. of the sensors, the output then indicates the dif- 
ference between the capacitance of the remaining 
sensor and that of the known capacitor. The empty 
capacitance of the sensors is known from previous 
measurements of the sensor in air. The difference 
between the dielectric constants of air and ullage gas 
is not significant for non-cryogenic propellants. For 
cryogenic propellants, a correction factor should be 
added to account for this difference. 


The slosh measuring system described in this 
article is a point measuring device; therefore, it can 
only be used when the shape of the liquid surface is 
known. The orientation of a planar liquid surface, as 
when slosh frequency is low relative to the first res- 
onant frequency, is determined as follows (Fig. 3). 

h - h 0 = ax + bz; (10) 

h 0 is the average propellant height, a is the slope of 
the liquid surface in the z plane, and b is the slope of 
the liquid surface in the x plane. The deviation of any 
point on the surface from the average propellant 
height, h - h 0 , is to be determined. 
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FIGURE 3. THREE-DIMENSIONAL REPRESENTATION OF 
PLANAR LIQUID SURFACE 

The axes x and z in Figure 3 are drawn along the 
^ bases of the two planes of liquid slosh measurement. 
Liquid height data are thus provided by sensors at 
known distances r along x and z from the origin (cen- 
ter line of the tank) ; the output of one pair of sensors 
in one plane is h 1 - h 2 ; the output of the other pair in 
the other plane is h 3 - h 4 . With these data, a and bmay 
be determined by the following equations. 

For the two sensors located at z = 0, x = +r and 
-r, 



For the two sensors located at x = 0, z = +r and 
-r , 



If the motion of the liquid surface is a form of 
periodic motion such as simple harmonic motion, the 
outputs of the two pairs of probes readily provide in- 
formation on slosh frequency, peak amplitude, and 
orientation. 


SENSOR CHARACTERISTICS 

The conventional capacitance sensor consists of 
rigid concentric tubes, open top and bottom, supported 
and insulated from each other by solid insulators. 
This is not the optimum design if fast response is re- 
quired as in slosh measurements. 

A parallel plate configuration was selected for the 
slosh sensor for a practical design to obtain fast re- 
sponse, minimum weight and size, and ease of as- 
sembly. 

A cross-section of the sensor is shown in Figure 
4, and a section of the sensor is shown in Figure 5. 
The two outer plates are connected as one electrode, 
and the center plate forms the other electrode. The 
shields confine the field mainly to the area between 
the capacitance plates and form a rigid structure for 
support of the probe. The parameters of the parallel 
plate construction cannot be accurately defined ana- 
lytically as can a concentric tube construction; but 
from the general equation for a parallel plate capaci- 
tor, the major problem can be seen to be in maintain- 
ing the plate area and spacing constant as a function 
of sensor length. The unique design of the sensor 
minimizes all extreme tolerance requirements with- 
out sacrificing performance. The use of a center 



FIGURE 4. CROSS SECTION OF PROBE 
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plate doubles the capacitance but requires very little 
additional weight. The spacing between the outside 
plates is large ( approximately 2. 5 cm) and is con- 
trolled by molded spacers. The center plate is nar- * 
rower than the outer plates; this feature minimizes 
the effect of lateral movement of the center plate in 
relation to the outer plates. The center plate is slot- 
ted to eliminate the effects of spacers on linearity. 
Relatively large displacements of the center plate to- 
ward the outside plates, which could occur under vi- 
bration environments, can be tolerated with very 
small effects on capacitance. This can be seen from 
the curve of capacitance versus displacement of cen- 
ter plate (Fig. 6). A 5 percent displacement causes 
only a 0.25 percent change in capacitance. The active 
capacitance of the probe is approximately 0. 4pF/cm. 
The capacitances to ground of the probe are relatively 
large, but their effects are eliminated by connecting 
the sensor as a three-terminal capacitor with the large 
shunt capacitors placed across very small imped- 
ances. 

Spacers are located at intervals of 2. 5 cm along 
die probe to maintain adequate spacing and structural 
integrity. These spacers form conductive paths be- 
tween the two elements of the capacitor that drastically 



FIGURE 6. SENSOR CAPACITANCE VS DISPLACEMENT OF 
CENTER PLATE 
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* lower the Q of the sensor. This problem is elim- 
inated by making conductive guards in the center of 
the spacers. The guards are grounded to the shields 
•by the bolts in the spacers. This technique eliminates 
the effect of conductive currents by shunting them to 
the bridge guard point in parallel with very low im- 
pedances. 

The sensor is designed in segments for easy 


installation. The segment interconnections are de- 
signed so no discontinuities are at segment joints. 

Two 10. 2 m sensors and a breadboard circuit 
were built and tested. The tests show- that an ac- 
curacy of approximately ± 1.5 cm is achieved. 

The technique is to be used on the S-IC vehicle 
to measure fuel slosh. 
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B. AN EMERGENCY HELIUM RUBBLING CONTROL 
SYSTEM FOR THE PREVENTION OF GEYSERING IN 
THE LOX SUCTION LINES OF S-IC by James A. Power 
and A. Turel* 

-—/Hssrrt 

An emergency helium bubbling system was devised 
for geyser suppression in the LOX suction lines of S-IC. 
Flight model pressure and temperature transducers were 
compared to more accurate transducers to determine the 
best sensor configuration. Transducers were electroni- 
cally simulated to calibrate the circuit design. A de- 
velopment system was tested for feasibility on a full- 
scale LOX flow line test facility. 

AlcthoAj 

INTRODUCTION 

Large space vehicles with long LOX suction lines 
to the engines have a common danger problem in that 
vaporization in the bottom of the lines may produce a 
geyser during the oxidizer loading and prelaunch pe- 
riod. A geyser occurs when the heat input to the bot- 
tom of the line causes vapor to form in excess of the 
normal bubble release. The increased rate of vapori- 
zation causes expulsion of the liquid and sudden grav- 
itational refilling of the line from the tank above, re- 
sulting in a pressure surge large enough to damage 
feed lines, valves, and line supports. 

A cross-feed circulation technique, thermal self- 
pumping, is used to prevent geyseringin the LOX suc- 
tion lines of S-IC. The circulation is initiated by bub- 
bling helium into one line only. The colder LOX flows 
down one line and through interconnect lines to a line 
with warmer LOX, displacing the warmer LOX up- 
ward to be cooled by the large volume of colder LOX 
in the tank. 

If in emergency cases thermal self-pumping is in- 
adequate to prevent geysering conditions, a method is 
needed to suppress geysering. One technique would 
be to continuously bubble helium into the lines. Enough 
helium would have to be injected to cope with the worst 
possible conditions, which for the entire prelaunch 
period would require a considerable amount of helium. 
Helium is in short supply and is considered a critical 
material. 


Astrionics Laboratory was requested to investigate 
the feasibility of using a system of flight instrumen- 
tation transducers and electronic circuitry to provide 
the necessary helium bubbling control. * 

The following goals were established for the sys- 
tem: 

a. LOX suction line ambient pressure should 
be monitored between 22.8 and 51.7 N/cm 2 , with max- 
imum accuracy in the range of 34. 5 to 5i. 7^ N/cm 2 . 

b. LOX suction line internal temperature 
should be monitored between -184. 4° and -156. 7° c. 

c. Helium bubbling should be initiated if the 
internal temperature of the suction line rises to with- 
in 0.56° C of the LOX vaporization temperature. 
(Since the initiation of this program, tests have indi- 
cated that the actual differential temperature used will 
be greater than 2. 78° C. ) 

The mechanical layout of the system is shown in 
Figure 1. 



HELIUM FLOW 


If a suitable control system could be devised, the 
helium would only be used in the case of an emergency. 


FIGURE 1. MECHANICAL LAYOUT OF HELIUM 
BUBBLING CONTROL SYSTEM 


Brown Engineering Company 


78 


# ¥ 

ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


PRESSURE 
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t 


1 


i 


REGULATED 
POWER 
SUPPLY 
PROVIDES 
DC VOLTAGES 
FOR SYSTEM 


COMPARATOR 
AMPLIFIER 
OUTPUT 
CAN ACTUATE 
A RELAY 


I f 
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^ 0 TO 5 VOLTS 
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.660V TO 1.720V 


SIMULATED 
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21 TO 69 N/cm 2 

DIODE 

FUNCTION 

GENERATOR 

-|74°C TO -159° C 

.75V TO 2.5V 

DERIVES 

VAPORIZATION 

TEMPERATURE 

.660V TO I.720Y 
ERRORiO.r C 



FIGURE 2. BLOCK DIAGRAM OF HELIUM BUBBLING CONTROL SYSTEM 


The following method is proposed for sensing the 
danger of geyser ing (Fig. 2): 

a. Pressure in the LOX suction line would be 
measured. 

b. The LOX vaporization temperature would be 
derived from the output of the pressure transducer. 

c. The temperature of LOX in the suction line 
would be measured. 

d. The internal temperature and the vaporiza- 
tion temperature would be compared. If the difference 
between the two is small, indicating that vaporization 
is likely, an alarm signal would be initiated. 

SYSTEM DESIGN DETAILS 

A study was conducted to determine whether meas- 
urements of pressure and temperature inside the LOX 


suction lines could be obtained with sufficient accu- 
racy, using flight transducers, for use in the helium 
bubbling control system. A circuit was designed for 
testing with flight transducers; then performance char- 
acteristics of the circuit and flight transducers were 
determined prior to operational testing. 

Pressure Transducer Evaluation. The pressure 
transducer specified for flight measurements is a 
bellows/potentiometer sensor with a range of 137. 9 
N/cm 2 (Fig. 3). The performance characteristics of 
this model transducer in several ranges are given for 
comparison in Table 1. Vibration and acceleration will 
contribute no significant errors during the prelaunch 
period when helium bubbling control is needed. 

A 2. 76 N/cm 2 error in the 137. 9 N/cm 2 transdu- 
cer is equivalent tol.ll°Cin the derived vaporization 
temperature. A diode function generator (DFG) is to 
be used with the pressure transducer to derive the 
vaporization temperature of LOX at a given pressure. 
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FIGURE 3. COMPONENTS USED IN FEASIBILITY TEST SYSTEM 


TABLE 1. PERFORMANCE CHARACTERISTICS OF 
THE PRESSURE TRANSDUCER DURING 
SATURN V PRELAUNCH PERIOD 


Pressure 

Static 

Static Error 

Error to 

Range 

Error Band 

Band Plus 

be Ex- 

( N,/cm 2 ) 

<%) 

Temperature 
Error Band 
(%) 

pec ted 
CN/cm 2 ) 

1.71 

±1.2 

± 2. 2 

± 1. 14 

68.95 

±1.0 

± 2. 0 

± 1. 38 

137.9 

±1.0 

± 2. 0 

± 2.76 


This DFG will contribute an error of ± 0. li°C over any 
50° C temperature range. Total error of the derived 
LOX vaporization temperature measurement is, there- 
fore, ±1.22°C. Use of a lower pressure range trans- 
ducer (Table 2) or use of a strain gauge transducer 
would improve accuracy. 

Temperature Transducer Evaluation. The tem- 
perature transducer used for Right measurements is 

a platinum resistance sensor (Fig. 3). The value of 
resistance of the unit is specified to be 1400 ohms 
±0.5 percent at 0°C; at any temperature, the repeat- 
ability of the unit is specified to be within ±0. 5 ohm. 


A dc amplifier assembly is specified for use with this 
temperature sensor. Amplifier performance specifi- 
cations for a period of 8 hours are as follows*. 

Linearity = ± 1% 

Gain Stability = ± 1% 

Output Di'ift = ± 0. 5% 

RSS Error = ±1.5%. 

Present Astrionics procedures allow for the tem- 
perature measurement to be obtained in several ways. 
One method of making the platinum resistance sensor 
measurements is to use uncalibrated gauges and an 
associated dc amplifier and to depend upon the above 
specifications to limit the error to the following ex- 
pression: 

Total RSS Temperature Error = 

( error ^ /error \ 2 /errorv Avire N “ 

of R 0 I ± (of read-j + j of dc ] +( leads! 

at 0° c) \ing J \amp J \ervox) 


(0. 5) 2 ± (0. 12) 2 ± ( 1. 5) 2 ± (wire leads error) 2 

~ 2 . 0 %. 

The wire leads connecting the platinum resistance 
element to the bridge produce an error that is difficult 
to evaluate. A second method of making the tempera- 
ture measurement is to use a three-wire bridge that 
puts a signal lead in each side of the bridge and re- 
duces the error effect of the leads. For this discus- 
sion, the non-caiibrated method using a three-wire 
bridge will be assumed to produce an accuracy of 1. 6 
percent. 

An improvement in the accuracy of the measure- 
ment should be obtained by calibrating the total meas- 
urement, both temperature gauge and dc amplifier, 
and by using a three-wire bridge. The total error for 
the calibrated temperature measurements would be as 
follows: 

Total Temperature Error =\/(0. 12) 2 + ( 1. 50) 2 

= 1.5%. 
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TABLE 2. PRESSURE ACCURACY ATTAINABLE 


Range and Type 
Of Transducer 

Total Error 
(±%) (± N/cm 2 ) 

Resulting Error 
in Vaporization 
Temperature (±°C) 

Total Error 
with DFG 
(±°C) 

137. 9 N/cm 2 

potentiometer 

type 

2.0 2.76 

i.i 

1.2 

69 N/cm 2 ' 

potentiometer 

type 

2.0 1.38 

0,6 

0.7 

69 N/cm 2 
strain gauge 

1.0 0.69 

0. 3 

0.4 

TABLE 

3. TEMPERATURE ACCURACY ATTAINABLE 

Methods of 
Measurement 

Error 

(±%) 

Total Error for 
19. 4° C Range 
(±°C) 

Total Error for 
27. 8°C Range 
(±°C) 

Non-calib rated 
two-wire bridge 

2. 0 

0.39 

0. 56 

Non-calibrated 
three-wire bridge 

1. 6 

0.31 

0. 44 

Calibrated 
two-wire bridge 

1. 8 

0.35 

0. 50 

Calibrated 
three-wire bridge 

1. 5 

0.29 

0. 42 

Calibrated 
three-wire with 
special amplifier 

0. 5 

0. 10 

0. 14 


Tf a separate temperature gauge and a higher qual- 
ity dc amplifier were used, an accuracy of ±0. 5 per- 
cent could be obtained. 

Effect of Transducer Accuracies on the Control 
System. The resulting accuracies of the pressure and 


temperature measurements are tabulated below tc 
study their effect on the sensitivity of the helium bub- 
bling control system. 

Note in Table 3 that, by decreasing the temper- 
ature measurement range, an accuracy improvement 
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in terms of degrees is achieved. The resultant error 
in vaporization temperature in the pressure channel, 
Ep, is the sum of the pressure transducer error and 

the DFG error. The error in the LOX temperature 
measurement, E , is the same as the error of the tem- 
perature transducer. Total error for the control sys- 
tem would be the sum of the errors of the pressure and 
temperature channels. 

Control Error = E + E,. 

P t 

To prove the feasibility of the system, transduc- 
ers were selected which would represent those speci- 
fied in the S-IC-i measuring program. A 137. 9N/cm 2 
pressure transducer and a platinum resistance ther- 
mometer using the standard flight model dc amplifier 
were selected. The control error would be pressure 
error plus temperature error, or 1.66°C, 

It is desired that the suction line temperature be 
no closer than acertain differential temperature (AT) 
from the vaporization point before helium bubbling is 
initiated. The AT originally selected was 0. 56° C. 
The control set point is the sum of AT and the control 
error, or 2. 22° C. Actual initiation of helium bubbling 
could occur at any point in the alarm range, which is 
within plus or minus the control error of the control 
set point (2.22° ± 1,66°C). The straight addition of 
pressure and temperature errors adds a measure of 
reliability to the accuracy calculations. 

The transducer evaluation indicated that, by using 
standard instruments while reducing the range of the 
measurement, an improvement in the control accuracy 
can be made. By using a 69 N/cm 2 pressure trans- 
ducer and reducing the range of the temperature meas- 
urement to 19. 4° C, the control errorwouldbe 0.96°C. 
The control set point would be 1.52°C, and the alarm 
rangewould be 1. 52° ± 0. 96°C. If expenditures were 
increased for the transducers and signal conditioning 
equipment, considerable improvement in the accura- 
cies could be obtained. Using a 69 N/cm 2 strain 
gauge pressure transducer, pressure error would be 
0. 39° C; and using a special high accuracy amplifier 
for the temperature transducer, temperature error 
would be 0. 11 °C. The control error , obtained by add- 
ing pressure error and temperature error, would be 
0. 5°C. The control set point would be 1.06°C, and 
the alarm range would be 1. 06° ± 0. 5° C. This alarm 
range represents the best accuracy that can be reas- 
onably expected using the proposed design for the he- 
lium bubbling control system. 


Recent tests have indicated that helium bubbling 
will be initiated at a differential temperature greater 
than 2. 78° C; therefore, flight measurement transduc- 
ers will probably be sufficiently accurate to obtain 
satisfactory results. 

Description of the Helium Bubbling Control Cir- 
cuitry. The helium bubbling control circuitry ( Fig. 3) 
includes the following components (Fig. 2) : 

Regulated power supply (two) 

Voltage follower (three) 

Diode function generator 
Comparator- amplifier 
Coefficient and control setting network 
Voltage selection and calibration network 
Miscellaneous circuitry 


Two identical transistorized regulated power sup- 
plies provide positive and negative 15 volts required 
by voltage followers and the comparator-amplifier; 
excitation voltage for the pressure transducer; and 
current for the diode function generator, internal cal- 
ibration network, and coefficient and control setting 
circuit. 

The three voltage followers used in the circuit are 
designated Al, A2, and A3 in the schematic (Fig. 4) . 
All three are solid state operational amplifiers and are 
used in the isolation amplifier ( gain =1) configuration 
to provide very high input and very low output imped- 
ances, thus eliminating the effect of loading on the op- 
eration. Adjustable resistors R40, R42, and R 55 are 
external offset zero trim potentiometers. 

The nonlinear relationship between the pressure 
and the evaporation temperature is derived by means 
of a diode function generator formed by resistors R46 
through R56 and the two diodes CR10 and CR11. A 
10,000 ohm thermistor compensates the ambient tem- 
perature variations. 

Resistors R15, R28, and zener diode CR9 with 
amplifier A4form the comparator-amplifier. A4, like 
Al, A2, and A3, is an operational amplifier with a gain 
of approximately 10,000. The theory of operation of 
the comparator-amplifier is as follows. Assume an 
open loop gain of 10,000 and -t 12 volts saturation for 
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the amplifier; if the diode were not across the ampli- 
fier, any voltage difference of 12/10,000 = 0,0012 
volts or greater would give an output of positive or neg- 
ative 12 volts depending upon the polarity of the input 
signal. The diode limits the negative output to approx- 
imately -0, 5 volts or less. If the diode has a zener 
voltage of 10 volts, output in the positive direction c an- 
notexceed 10 volts. During the actual operation of the 
system, the signal (vaporization temperature) applied 
to the amplifier through Rl5 is normally higher than 
the signal (temperature) applied through R28 by at 
least 1.2 millivolts. Therefore, the output of the 
amplifier is approximately -0. 5 volts. As the emer- 
gency condition is approached, the magnitude of the 
signal coming through R28 increases; this increase, 
if further amplified by the comparator, causes the 
output to swing from -0. 5 to + 10 volts almost in- 
stantly. 

The coefficient and control setting network is 
formed by an attenuator consisting of resistors R24, 
R25, R26, and R27, and a biasing arrangement con- 
sisting of resistors R29 through R34. R28 and R35 act 
as adding resistors by isolating the two networks. By 
means of S5, four discrete temperature differential 
settings can be obtained; by means of R36, any tem- 
perature differential setting between 0° and 5, 6° C can 
be obtained. 

All important voltages in the operating system can 
be selected and read with a voltmeter at S3. The cal- 
ibration network formed by resistors R16 through R21 
is used to verify the proper operation of the circuit 
and to adjust the amplifier. The two switches desig- 
nated as HI CAL and LOW CAL are used to calibrate 
the temperature transducer/amplifier combination. 

Voltage dividers R44, R45 and R22, R23 provide 
an attenuation of 1000 for recording purposes. Re- 
sistors R37 and R38 also attenuate the output signal 
by a factor of approximately 200 for recording pur- 
poses . 

System Performance Characteristics. The 
137.9 N/cm 2 potentiometric flight~model pressure 
transducer was calibrated over its operational range. 
The resulting calibration curve is presented in Fig- 
ure 5. 

The platinum resistance flight-model temperature 
transducer was also calibrated over its operational 
range, -185° to -155°C, and found to be nonlinear. 



FIGURE 5. FLIGHT PRESSURE TRANSDUCER CALIBRATION 

The resulting calibration curve is presented in Fig- 
ure 6. 

The system's circuitry performance character- 
istics were obtained during operation with the trans-^ 
ducers electronically simulated. The pressure trans- 
ducer was simulated by a 2000 ohm potentiometer in 
series with a 3000 ohm resistor; 5.076 volts was. 
applied across this series combination at the source 
end of the cable. The potentiometer output was ad- 
justed manually to give different voltage values corre- 
sponding to pressure values from 20.7 to 55.1 
N/cm 2 . The temperature transducer was simulated 
by a decade resistance box. At each pressure setting, 
the value of resistance was varied until an output sig- 
nal (approximately half scale deflection) was observed 
on the front panel meter. When the signal appeared, 
the voltage at pin 7 of S4 was read and recorded. The 



FIGURE 6. FLIGHT TEMPERATURE TRANSDUCER CALIBRATION 


84 



ASTRIONICS RESEARCH AND DEVELOPMENT REPORT NO. 3 


corresponding temperature was obtained from the 
calibration curve of the temperature transducer . This 
procedure was repeated at differential temperature 
settings of 0° , 3°, 4.2°, and 5.8°C at each pressure 
value from 20.7 to 55.1 N/cm 2 . System calibra- 
tion data are plotted in Figures 7A and 7B. Notice 
that all four differential temperature setting curves 


do not maintain a constant difference of temperature 
in relation to the vaporization curve; these differences 
get smaller as pressure increases. This condition is 
caused by the nonlinearity of the temperature meas- 
urement. When the diode function generator was de- 
signed, the temperature measurement was assumed 
to be linear. 



FIGURE 7 A. SIMULATED SYSTEM CALIBRATION 



FIGURE 7B. EXPANDED PORTION FROM FIGURE 7A 
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OPERATIONAL TESTING OF THE 
DEVELOPMENT SYSTEM 

Operational tests were performed to demonstrate 
the feasibility of the entire helium bubbling control 
system. The flight measurement pressure and tem- 
perature transducers were mounted in a LOX suction 
line test facility that simulated the vehicle LOX suction 
lines and tank. Three successive control setpoints, 
5. 8° , 4.2°, and 3° C below LOX vaporization temper- 
ature, were established by a switch in the system cir- 
cuit. At each set point, the system correctly emitted 
an alarm signal within the prescribed alarm range as 
LOX inside the test suction lines warmed toward va- 
porization temperature. In an operational system, 
the alarm would have released gaseous helium into the 
suction lines to reduce the danger of LOX vaporization 
and subsequent geysering. These results are pre- 
sented in Table 4. 


CONCLUSIONS 

The results of the theoretical study and the feas- 
ibility test indicate that a system can be developed 
which would reliably initiate helium bubbling into the 
S-IC LOX suction lines to prevent geysering. 


In a test using the least sensitive of the systems 
considered, i. e. , a 137.9 N/cm 2 pressure transducer 
and a 28° C range for the temperature measurement, 
a maximum error of 0. 35° C resulted. This result in- 
dicates that the calculated values could be depended 
upon over the necessary period of time and over a 
large variation of ambient temperature. 

A development system produced for a feasibility 
test does not necessarily, however, incorporate all 
the features that an operational system should have. 
Features not included in the feasibility model would be 
required in the operational system. The units should 
be self-checking in case of malfunction; an alarm 
should be initiated, the malfunctioning unit automati- 
cally turned off, and the system switched to manual 
operation. A pressure transducer and a temperature 
transducer should probably be installed in each suction 
line. A "signal processor" may therefore be desirable 
to select the transducer combination approaching or 
nearest the alarm condition and to provide an analog 
display of the two signals. All transducers should, 
however, be connected to initiate helium bubbling. 

A successful operational helium bubbling control 
system will allow the conservation of helium, which is 
considered a critical item in short supply. Helium 
release would be initiated only as needed, rather than 
continuously as in some space systems. 


TABLE 4. FEASIBILITY TEST RESULTS 


Control Set 
Point (°C 
Below LOX 
Vaporiza- 
tion Point) 

Time From 
End of Thermal 
Pumping to 
Alarm Signal 
( seconds) 

LOX 

Line 

Pressure 

(N/cm 2 ) 

LOX 

Line 

Tempera- 
ture (°C) 

LOX 
Vapori- 
zation 
Tempera- 
ture (°C) 

System 

Error 

(°C) 

5. 8°C 

220 

30.48 

-176 

-170 

-0.2 

4. 2° C 

430 




-0. 35 

3. 0°C 

550 

29.79 

-173.6 

-170. 3 

-0. 3 
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C. AN UP-CONVERTER /DOWN-CONVERTER METHOD 
OF DATA TRANSMISSION by Edmund H. Gleason 

, N65 -23792 
4 ^ % 

An up-converter/ down-converter wideband micro- 
wave relay link is presented. The system involves con- 
verting VHF telemetry signals to microwave frequencies 
for transmission and then converting them back to VHF 
at the receiving site for demodulation. 

Also discussed are the advantages and disadvantages 
of the system as compared to existing means of relaying 
telemetry data via microwave links. 


24 VHF telemetry links in the 215 to 260 MHz band, 
causing checkout requirements and problems to be- 
come quite complex. Allowing each onboard trans- 
mitter to radiate from the vehicle antennas would re- 
quire radio frequency silence in the vicinity of the 
launch area. If the checkout is prolonged, other ve- 
hicle checkouts would be delayed. 

A proposed system to eliminate these problems 
involves up-conversion of the VHF telemetry frequen- 
cies, which can be transmitted within a tight direc- 
tional beam. To prove the feasibility of such a scheme, 
an experimental system was devised to transmit quick- 
look SA-6 telemetry data from the Green Mountain 
receiving site to the Astrionics Laboratory ground 
station at the Marshall Space Flight Center, a dis- 
tance of about ten miles. 


INTRODUCTION 

As launch vehicles increase in size and complex- 
ity, additional telemetry links are required to handle 
transmission of the increasing number of measure- 
ments. For example, Saturn V will have as many as 


SYSTEM CONCEPT 

The only change required within the vehicle is 
the addition of a coaxial switch, a directional cou- 
pler, and a dummy load (Fig. 1). For checkout, the 



Fi. 

F 2 . 

f 3 


FIGURE 1. VHF— TO— MICROWAVE UP-CONVERTER/TRANSMITTER FOR LAUNCH-PAD INSTALLATION 
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coaxial switch is placed in the TEST position and a 
portion of the power is sampled by the coupler and 
routed to the umbilical tower. These signals are ap- 
plied to a mixer along with the output of a stable local 
oscillator (stalo) operating at, for example, 7. 7 GHz. 
The mixer output is filtered to allow only the upper 
sideband frequencies in the 7. 915 to 7. 96 GHz band to 
pass on to the traveling wave tube amplifier. The 
amplified signals are then transmitted via a micro- 
wave link to the receiving site. 

The receiving end is shown in Figure 2. The in- 
coming microwave frequencies are heterodyned with 
a 7.7 GHz stalo signal to recover the original VHF 
telemetry frequencies, which are then amplified and 
routed to individual receivers for demodulation and 
recording. 

The advantages of such a system are obvious. 
Only the up-converted microwave frequencies are al- 
lowed to radiate and these are kept within a tight beam 
directed to the receiving site. Simultaneous checkout 
of several vehicles can take place simply by assign- 
ing different microwave bands to each pad. Checkout 
accuracy is maintained because the up-conversion and 
down-conversion processes do not deteriorate the 
signals in any way. 


Solid state multiplier chains that provide stabil- 
ities of one part in 10? long term are available. Sta- 
bilities one magnitude better are within the state of 
the art. The resulting long term system drift, as- 
suming the worst case, would not exceed 1600 Hz, 
which is well within the tolerances of telemetry re- 
ceivers having a bandwidth of 500 kHz. 

SYSTEM DESIGN 

Output Sampling 

The first requirement of the transmitting system 
is to sample the outputs of each VHF telemetry link 
within the vehicle. This feature can be accomplished 
by installing a remotely controlled, single pole, dou- 
ble throw coaxial switch in the telemetry antenna line 
and routing the sampled outputs through the umbilical 
disconnect to a level set attenuator. The attenuator 
can then be adjusted to provide a nominal 0. i milli- 
watt (-10 dbm) per channel signal for application to 
the up-converter. 

Up-Converter 

The up-converter consists of an Orthomode bal- 
anced mixer and a local oscillator. Figure 3 shows a 
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FIGURE 3. EFFECT OF LOCAL OSCILLATOR POWER ON MIXER 
CHARACTERISTICS 


family of curves of input versus output characteris- 
tics for different local oscillator power input levels. 
The output of the mixer, produced by heterodyning 
action, consists of the first order upper and lower 
sidebands, the suppressed carrier, and the higher 
order sidebands. 


The mixer characteristic becomes more linear as 
local oscillator power is increased (Fig. 3). It is 
important that the most linear characteristic be used 
to minimize cross modulation effects between chan- 
nels. Therefore, a local oscillator drive level in the 
range of 64 to 150 milliwatts should be used to oper- 
ate the mixer in the linear region. Diode burnout 
level is approximately 300 milliwatts; the maximum 
drive level is limited to 150 milliwatts by imposing a 


safety factor of two. Assuming a maximum of 24 
channels, there will be a total VHF signal power of 
2. 4 milliwatts (+3. 8 dbm) applied to the mixer. Fig- 
ure 3 shows that for this level the mixer will still be 
operating within the linear portion of the curve for a 
64-milliwatt local oscillator drive level and the re- 
sulting upper sideband output level will be 0. 045 
milliwatts (-13. 5 dbm). 

To achieve the desired stability of one part in 
10? a local oscillator source consisting of a crystal- 
controlled varactor multiplier chain is preferred. A 
klystron, phase locked to a crystal oscillator refer- 
ence, could also be used but requires much higher 
voltages and would be much larger. 

Since only the first order upper sidebands are to 
be transmitted, a bandpass filter must be provided at 
the mixer output to adequately attenuate the local os- 
cillator signals, the lower sidebands, and the higher 
order sidebands. The filter must have a 45 MHz flat 
bandwidth. All signals outside this passband must be 
attenuated at least 60 db, since the filter is followed 
by a wideband power amplifier. 


Power Amplifier 

To amplify the 0. 045 milliwatt (-13. 5 dbm) per 
channel signals to 100 milliwatts (+20 dbm) , a trav- 
eling wave tube amplifier with 33. 5 db gain is re- 
quired. This amplifier must be operating below sat- 
uration to minimize any cross modulation effects be- 
tween channels. Amplifiers with these capabilities 
are presently available with mean-time-before-failure 
rates higher than 10,000 hours and with completely 
solid state power supplies. 


Receiving System 

The receiver must be capable of receiving the 
transmitted signals, down converting to VHF, and pro- 
viding outputs for each of the VHF telemetry receiv- 
ers required for checkout of the vehicle. 

To prevent interference from other microwave 
sources, a preselector with characteristics identical 
to those of the filter in the transmitter is used. A 
low noise mixer preamplifier and a local oscillator 
source is required to convert the microwave signals 
back to VHF frequencies. The local oscillator must 
have the same frequency output and the same stability 
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as the one used in the transmitter. However, the re- 
quired power output is only 1 to 3 milliwatts, thereby 
reducing the size of this device as compared to the 
transmitter local oscillator. 

To provide the outputs for the telemetry receiv- 
ers, a passive multicoupler is used at the output of 
the mixer preamplifier. The receivers are then tuned 
separately to each of the 24 channels, enabling check- 
out of the systems. 

EXPERIMENTAL SYSTEM 

To substantiate the feasibility of the proposed 
system, an experimental link was constructed to re- 
lay telemetry signals received during the SA-6 launch. 
The VHF-to-microwave up-converter/ transmitter 
shown in Figure 4 was located at the Green Mountain 
receiving site some ten miles from Astrionics Labo- 
ratory at MSFC. The principle was the same as that 
of the launch pad scheme except that the signals from 
the vehicle were received by a tri-helix antenna, am- 
plified, and then converted to microwave frequencies 


for transmission to Astrionics Laboratory. The 
receiver /down-converter is shown in Figure 4 and is 
identical to the one proposed for launch pad checkout. 

Since low level received signals were being dealt 
with in this case, maintaining a desirable signal-to- 
noise ratio for the overall system required use of a 
VHF preamplifier with a noise figure equal to or 
better than that of the receivers used. Also, the 
overall gain in the transmitting system was greater 
than the loss caused by the frequency conversion 
process, the insertion loss of the components used, 
and the microwave path loss. 

Figure 5 shows the measured system signal-to- 
noise characteristics that were determined by apply- 
ing an FM modulated signal to the VHF preamplifier 
at the transmitting end and measuring the video signal- 
to-noise ratio at the receiving end. The resultant 
characteristics are nearly identical to those of a typ- 
ical telemetry receiver, thus demonstrating the supe- 
riority of this system to a typical microwave link 
utilizing subcarrier channels with no better than a * 
30 db signal- to -noise ratio per link. 



f, = 240 MHz 
f 2 = 250 MHz 
f 3 = 260 MHz 


F, » f, + fLO = 7.94 GHz 
F a - f 2 + fLO = 7.95 GHz 
F 3 = f 3 + fLO = 7.96 GHz 


note: frequencies 

SHOWN ARE 
EXAMPLES 


FIGURE 4. VHF— TO— MICROWAVE UP-CONVERTER/TRANSMITTER (GREEN MOUNTAIN) 
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FIGURE 5. FM-MODULATED VHF SIGNAL APPLIED AT GREEN 
* MOUNTAIN AND RECEIVED AT ASTRIONICS 

LABORATORY 


The only disadvantage of this system is high level 
interference in the 240 to 260 MHz band. Because 
this is picked up by the tri-helix antenna, saturation 
of the amplifier occurs, blocking out the lower level 
telemetry signals. However, RF silence was main- 
tained during the SA-6 flight and the system provided 
excellent quick-look data to the Astrionics Laboratory 
ground station. 

It should be noted that no interference problem 
would result in the launch pad scheme because the 
VHF signals are carried via n hard line' 1 and the mi- 
crowave frequencies are held within a directional 
beam. 


CONCLUSIONS 

1. The proposed microwave relay system indi- 
cates capabilities for eliminating problems encount- 
ered during telemetry checkouts for large boosters. 

2. The proposed system is simple to devise, 
should be relatively trouble free, and should not in- 
troduce any inaccuracies in the checkout results. 
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PUBLICATIONS 
AND PRESENTATIONS 


k. PUBLICATIONS 


SATURN INSTRUMENTATION SYSTEMS, J.T. Powell, 
ISA Journal, November 1964 (also presented at the 10th 
National Aerospace Instrumentation Symposium, New 
York, N. Y., May 4-6, 1964) 

A brief description of the Saturn vehicles is given, 
delineating the makeup of and differences between the 
Saturn I, Saturn IB, and Saturn V. Certain ground rules 
and constraints for instrumentation system design are dis- 
cussed. The measuring system (transducers and signal 
conditioning) is outlined, showing the signal paths into 
the telemetry systems. Various telemetry modulation and 
multiplexing methods are defined, and the combining of 
these techniques into a typical Saturn telemetry system 
is illustrated. Special emphasis is given to digital tech- 
niques and their interrelationship to vehicle checkout 
procedures. Radio frequency tracking systems used are 
briefly discussed, giving the characteristics of each and 
the coverage obtained. Command systems for range safety 
and guidance are explained in some detail. The use of 
optical instrumentation (television and film cameras) is 
covered, showing the arrangement planned for Saturn V. 
The "operational " Saturn V vehicles to be used for moon 
launches have somewhat reduced instrumentation and are 
briefly discussed. 

STABILITY OF A THIRD ORDER NONLINEAR DIFFER- 
ENTIAL EQUATION, Bernard Asner, Jr., IEEE Trans- 
actions on Automatic Control, October 1964, Vol . AC-9 

Utilizing the second method of Liapunov, the sta- 
bility properties of a third order nonlinear differential 
equation are investigated. The Liapunov function used 
establishes asymptotic stability in the large if the non- 
linear function is restricted to a sector for which the 
lower bound satisfies an inequality analogous to the 
Routh-Hurwitz conditions. 


RADIO FREQUENCY EVALUATION OF SA-4 VEHICLE, 
Olen Ely, MTP-ASTR- l-64-3,~J^^7'75, 1964 

Al I radi o frequency systems on board the SA-4 vehi- 
cle functioned properly throughout powered flight. The 
systems were still operating at the time of reentry; but 
some data were lost on the tracking systems because of 
the rolling and tumbling of the vehicle. This was ex- 
pected since the tracking system used only one antenna. 

Good telemetry data were received at the Cape 
Telemetry 2station throughout the entire powered flight, 
including retrorocket ignition. Some of the smaller sta- 
tions experienced interruptions during this time. 

Main engine flame attenuation was similar to pre-^ 
vious flights. Retrorocket attenuation was more severe 
in some instances than on SA-3. The most significant! 
result of retrorocket firing was that it caused the signal 
to be attenuated and noisy for approximately five to six 
seconds; yet the retrorockets burned for only two seconds . 

This was the first flight carrying the radar altimeter 
and the Mistram tracking system. The roll of the vehi- 
cle causedsome problems with the radar altimeter; how- 
ever, the equipment performed satisfactorily. A com- 
parison of the data from the Azusa system, the UDOP 
system, and the Mistram system shows that the data 
agreedwithin three to five meters during powered flight . 


DISCRETE ACTIVITY INDICATOR SYSTEM, George A. 
Bailey, NASA TM X-53045, May^,“l96? 

A discrete activity indicator system processes 768 
discrete signal inputs and presents information to a com- 
puter only when a change in an input occurs. Resolution 
time of the system is approximately 1 .2 milliseconds for 
1000 discrete inputs. 

A basic clock in the system drives a functional clock 
which, in turn, drives a binary counter and six 8-phase 
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drivers in parallel. Each 8-phase driver operates a 128- 
way current-steering switch, which sequentially interro- 
gates the cores of a 128-bit memory bank. If a change 
occurs in a discrete signal input (from 0 to 28 V dc or 
28 to 0 V dc), an output is obtained from the associated 
memory core (when the core is interrogated) and applied 
to a sense amplifier. The outputs from the sense ampli- 
fier and the binary counter are then applied to a system 
register for transfer to the computer. These signals allow 
the computer to determine which discrete signal changed 
and the direction of the change. When a change occurs, 
a priority-interrupt circuit stops the interrogation until 
the computer processes the information. After the infor- 
mation is processed, a clock-start multivibrator allows 
the interrogation to resume. 

The power requirement for the system is reduced 
75 percent by using a ground commutator to turn on the 
discrete drivers of the memory in groups of 16 when each 
particular group is being interrogated. 

Positive synchronization between the binary counter 
^ind the current-steering switches is provided by a binary 
synchronizer. 

I 

A NONDESTRUCTIVE MEMORY BUFFER USING 
SINGLE APERTURED CORES, - George A. Bailey, NASA 
TM X-53057, June 4, 1*6* 

A nondestructive memory buffer using single aper- 
tured cores is described. The principle of operation is 
that of using read currents and a decoding device to re- 
write into previously read cores. This memory operates 
faster and consumes less power than those using multi - 
apertured cores. The supporting logic system features 
unconventional use of magnetic devices to minimize the 
need for semiconductors. 


INFLIGHT GRAVITATIONAL COMPUTATION, Joseph 
T. "Dabbs, NASA TM X-53061 ," June 9, 1964 

A generalized scheme of gravitational computation 
was considered for solution by airborne digital vehicle 
computers. The only assumption was that the coordinate 
systems used were rectangular, orthogonal, and right- 
hand. A computational scheme was derived that should 
be valid for any Saturn class vehicle computer and re- 
quires only the change of certain constants for any change 
in mission. 


An added advantage of this scheme is that by in- 
serting different constants the gravitation function of the 
moon would be readily computable with the same set of 
equations . 

The scheme presented is readily programable and 
has inherent scheme accuracy sufficient for any foresee- 
able missions. 

RADIO FREQUENCY EVALUATION OF SA-5 VEHICLE , 
Olen Ely, R^ W. Hockenberger, Parley Howell, and 
Geddes Boone, NASA TM-X53073, June 10, 1964 

All radio frequency systems on the SA-5 vehicle 
performed satisfactorily except the radar altimeter, which 
was unable to remain in a locked condition. The only 
other discrepancy was an instantaneous drop of between 
8 and 10 db in the C-band radar signal level at 472.8 
seconds . 

Signal strength degradation caused by main engine 
flame was less severe than on any previous Saturn flight. 
The level of attenuation varied from station to station, 
reaching a maximum of 26 db on VHF telemetry at the 
Cape Telemetry 2 station. UDOP signal strength data 
showed 'a maximum attenuation of 15 db. Higher fre- 
quency equipment, such as Mistram, Azusa/Glotrac, and 
C-band radar, was not significantly affected by the main 
engine exhaust plume. 

Retrorocket attenuation was very severe, causing 
all signals to drop to threshold. However, C-band radar 
recovered before retrorocket burning was completed. 
Ullage rockets alone did not affect the signal level of 
any RF system. 

PRECISION OPTICAL TRACKING SYSTEM FOR AD- 
VANCED LAUNCH VEHICLES , "Charles L. Wyman, 
Klaus Juergensen, Robert L. Kurtz, Larry Hayes, and 
John M. Gould, NASA TM X-53076, June 25, 1964 

The development efforts on an experimental optical 
radar system for precision tracking for advanced launch 
vehicles during the early launch phase are described. 
Basic concepts for a prototype system are discussed along 
with some advanced techniques that are being investi- 
gated. Analyses and investigations concerning the proj- 
ect are discussed. An error analysis is made assuming a 
tracker that measures elevation angle, azimuth angle, 
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and distance in digital form. The analysis is concerned 
with errors in distance, velocity, and acceleration com- 
puted from distance measurements. The distance meas- 
urements are assumed to be inaccurate by a systematic 
error and a distributed error. Also, the error caused by 
linear interpolation between displayed points is consid- 
ered. No data smoothing is taken into account. Equa- 
tions are developed that can be used to obtain optimal 
data rate and highest frequency with which the true 
range and angular parameters may change without causing 
errors exceeding the required accuracies. 

A trajectory analysis is performed utilizing a typi- 
cal Saturn V trajectory to determine the dynamic range 
of the instrument parameters. Vehicle and instrument 
parameters are plotted against time, and the curves are 
used for design purposes. 

A range analysis develops optical radar equations. 
The system uses retrodirective corner reflectors on the 
booster to increase the return signal. The increase is 
calculated and included in the range equations along 
with calculations of atmospheric attenuation. A pre- 
liminary calculation of background radiation is made. 

A preliminary synthesis of the tracking mount servo 
assumes a type one servo with tachometer stabilization. 
The analysis indicates that the system has the desired 
static accuracy in the presence of high breakaway fric- 
tion, can operate at very small dynamic errors, and ex- 
hibits good stability. 

An investigation concerning refractive and turbu- 
lence effects of the atmosphere on laser beams is also 
discussed. 


DEVELOPMENT OF SATURN FIRE DETECTION SYS- 
TEMS , H. D. Burke, R. Stockard, and J. E. Zimmerman, 
NASA TM X— 53 141 , September 24, 1964 

The research and development leading to the ad- 
vancement of fire detection techniques are discussed, 
with emphasis on the flight fire detection system. This 
system, which has been designed for flight operation, is 
an extension of the prelaunch fire detection system pres- 
ently in use on the Saturn vehicle. 

The prelaunch fire detection system is given con- 
siderable attention, as well as the fire detection sensing 
element. Both systems use a thermocouple gauge as the 


sensing element and operate on the basis of the temper- 
ature rate of rise at the gauge, rather than the absolute 
temperature . 


EVALUATION OF ADHESIVES FOR INSTALLING FLAT 
CONDUCTOR CABLES, Bobby W. Kennedy7 NASA 
TM X-53080, October 8, 1964 

Fourteen adhesives for mylar-insulated flat conductor 
cable were tested for creep and peel strength, corrosion 
effects, heat distortion, vacuum effects, and cryogenic 
temperature effects. 

Minnesota Mining and Manufacturing Company's 
1099 adhesive was found to be the best liquid adhesive. 
Fasson Products Corporation's S-277 adhesive was found 
to be the best pressure sensitive adhesive. Dow Coming 
Company's Silastic 140 is suitable in areas where initial 
tack is not a problem. 

4 

SPACE VEHICLE SA-4, TELEMETRY SYSTEM, E. H. 
Reeves, Jr., J. R. Stovall, and W. B. Threlkeld, Jr. 
NASA TM X-53153, October 26, 1964 1 

The performance evaluation of the complete telem- 
etry system (ten links) used for flight testing Saturn 
vehicle SA-4 is presented. The eight operable telemetry 
links, the two experimental telemetry links, and the 
telemetry auxiliary equipment assembly have been tech- 
nically analyzed on an individual basis. 

Statistical analyses were performed on much of the 
telemetry data, and the results of these analyses are 
presented . 

No RF signal fade or dropout difficulties were en- 
countered on any of the links. 

A malfunction of the PCM system occurred during 
the flight. An analysis of the cause and effect of this 
malfunction is included. 

It is concluded that the overall performance of the 
telemetry system used for flight testing SA-4 was as 
anticipated. It is inferred that the derived test data 
will facilitate the development of more efficient telem- 
etry systems . 
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SECURE RANGE SAFETY COMMAND SYSTEM FOR 
SATURN , H. R. Lower/, NASA TM“1<-53162; 

. November 9, 1964 

A high alphabet range safety command concept 
using digital methods compatible with the present modu- 
lation technique used on the Eastern Test Range is dis- 
cussed. 

The new system provides a high degree of protection 
against intentional interrogation by unfriendly intruders 
and against unintentional interrogation (false alarms) by 
noise. The security against intentional interrogation is 
measured by the fraction of the total code combinations 
that an intelligent unfriendly interrogator, having un- 
limited technical resources, would be able to transmit 
during equipment access time. It is assumed that all sys- 
tem parameters are known to him except the code-of- 
the-mission. The security against unintentional or ran- 
domly generated false alarms is, as it should be, much 
greater than against intentional, intelligent interroga- 
tion. Operational difficulties of handling the security 
aspects of the system are minimized by using a remov- 
Ible , inexpensive, easily stored code plug to change the 
code-of-the-mis$ion . All hardware is unclassified ex- 
cept the actual flight code plugs. 

In contrast to a binary system, which utilizes a 
language consisting of only two distinct characters, the 
new secure range safety command system uses a language 
consisting of 21 characters. Each character of the 21- 
charaoter alphabet is formed from a subalphabet of seven 
symbols taken two at a time. Each symbol is a burst of 
an audio-frequency tone. A pair of tones transmitted 
simultaneously comprise a character. This high-alphabet 
technique provides the necessary high security without 
sacrificing essential performance characteristics of the 
present system. Further, the new high-alphabet system 
is, by design, tolerant of relatively wide radar or pulse 
interference . 

The security is realized mainly by the message for- 
mat, which consists of a nine-character address word, 
followed by a two character function (command) word. 
The limited time of access to the equipment is the final 
factor in the security. 

Each character of the nine-character address word 
is unique within a given address, and therefore the new 
system can operate without the necessity of transmitted 


clock or reference synchronizing information. The sym- 
bols are, however, synchronous with the character inter- 
val period (and with each other) and are spaced and 
phased to minimize intersymbol interference, which could 
create unwanted characters. 

The main RF carrier is frequency modulated by a 
subcarrier system which, in turn, employs a multiple fre- 
quency shift (MFS) technique. The term frequency shift 
keying (FSK) is also used to describe the technique. The 
portion of the baseband occupied by the subcarrier sys- 
tem is approximately the same as that now reserved and 
used for range safety purposes. The nominal tone fre- 
quencies, however, are not the same as the Inter-Range 
Instrumentation Group (IRIG) tone channel frequencies. 
The tone-frequency spacing is, by design, an integral 
multiple of the character repetition rate. This makes 
possible a simple phase-coherent tone keying scheme. 

The two main vehicular components of the new 
secure range safety command system (receiver and de- 
coder) are presently being manufactured. Provided all 
present schedules can be met, prototype receivers and 
decoders will fly as IU passengers on vehicle SA-8 and 
several later vehicles. Upon qualification by the Range 
Safety Activity at PAFB, the new system will become an 
active system on all Saturn vehicles, replacing the pres- 
ent non-secure equipment. 

A system design has been evolved which has a high 
degree of security, will work with reliability in its ex- 
pected communication environment, and lends itself to 
implementation of dependable hardware. 


RADIO FREQUENCY EVALUATION OF SA-6 VEHICLE , 
Olen Ely, R. W. Hockenberger, Parley HoweTl, and 
Geddes Boone, NASA TM X-53170, November 24, 1964 

All RF systems performed satisfactorily with the ex- 
ception of C-band radar, which was approximately 10 to 
20 db below the theoretical values of signal strength. 
This system suffered an additional 10 db drop in signal at 
480 seconds; however, it recovered at 560 seconds. This 
was almost identical to the phenomenon that occurred on 
SA-5 at 472. 8 seconds. This malfunction has been attri- 
buted to arcing either in the cabling connecting the ra- 
dar beacon to the antenna or in one of the connectors. 
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Signal strength degradation caused by main engine 
flame was very similar to that of SA-5. However, un- 
like SA-5, some of the signal strength data contained 
noise caused by main engine exhaust until inboard engine 
cutoff at 143 seconds. 

Retrorocket ignition caused the signal level of all 
RF systems to drop below threshold. The C-band radar 
system recovered within 2.4 seconds, prior to any other 
system . 

The ullage rocket exhaust did not noticeably affect 
the signal propagation. It is possible that the severe 
attenuation of the RF signals is caused by combining the 
ullage rocket and retrorocket exhausts. 


EVALUATION OF POLYMERIC MATERIALS FOR ELEC- 
TRONIC^ Bobby W. Kennedy, NASA 

TMX-53181, December 16, 1964 

Certain potting and conformal coating materials 
were contributing to the failure of critical glass-sealed 
components and solder joints in cordwood modules and 
on printed circuit boards. The cause of these failures 
was determined to be differences in coefficients of ther- 
mal expansion for materials. Analysis indicated the ne- 
cessity ofselecting materials with a lower coefficient of 
thermal expansion for potting cordwood modules and 
changing from a hard, rigid epoxy to a soft, flexible 
polyurethane for printed circuit board conformal coating . 

Of the four potting materials tested, Emerson and 
Cuming's Stycast 1090 epoxy was selected for potting 
cordwood modules. 

Three of the coating materials were too thin to sup- 
port components; therefore, no further testing of these 
materials is anticipated. However, these materials may 
be used when a thin protective coating is needed and 
when component support is not a problem. 

Products Research-1538 polyurethane material is re- 
commended and is being used for coating printed circuit 
boards; however, in high frequency circuits, an epoxy 
may be used after encapsulating the component with 
vinyl sleeving to prevent RF losses. 


B. PRESENTATIONS 


THE APPLICATION OF Z-TRANSFORM THEORY TO 
THE ANALYSIS OF SWITCHED-TYPE NONLINEAR 
SYSTEMS , E. S. McVey* and G. S. Nurre, presented 
at the IEEE Symposium, New York, N. Y., March 23, 
1964 

The stability of certain discontinuous feedback sys- 
tems is analyzed by the use of the z-transform theory. 

It is shown that there can be found in the z-plane a sys- 
tem which is equivalent to the nonlinear system under 
steady state conditions. The existence of any limit cycle 
can be found together with the frequency at which it 
occurs. An intuitive argument for the stability of a limit 
cycle is also presented. Several examples are worked 
and the results are compared with those obtained by other 
methods. The z-transform method is exact and is not re- 
stricted by the order of the system. 

A SELF-ADAPTIVE TRACKING FILTER , Michael T.| 
Bore 1 1 i and Hans H. Hosenthien, presented to the SAE 
Aerospace Vehicle Flight Controls System Committee, 
Missile and Pilotless Aircraft Subcommittee, New York,* 
N. Y., July 8-10, 1964 

This presentation describes an adaptive tracking 
notch (ATN) filter that has the capability of either gain- 
or phase-stabilizing the low frequency structural bending 
modes of large aerospace vehicle boosters, such as Saturn 

V. The basic principle of the ATN filter is to construct 
a re versed-phase bending signal, utilizing modulation 
techniques, and then to sum this constructed signal with 
the sensor output signal , suppressing the undesired bend- 
ing signal for gain-stabilization or shifting its phase for 
phase-stabilization. A brief discussion of the ATN fil- 
ter's characteristics is presented. Results from analog 
computer simulations of the Saturn V vehicle utilizing 
breadboard ATN filters are reviewed and evaluated. 

DEVELOPMENT OF SERVOVALVES WITH IMPROVED 
RELIABILITY FOR SPACE VEHICLES , M. A. Kalange, 

W. H. Pol lock**, and W. J. Thayer***, presented to the 
SAE Committee A-6, Aerospace Fluid Power Technolo- 
gies, Boston, Mass., September 14-18, 1964 

Considerations for improvement in the reliability of 
Saturn engine gimbal servosystems are briefly covered. 

*U. S. Army 
**Cadillac Gage Co. 

***Moog Servo Control, Inc. 
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Saturn I servovalves operate with increased electrical 
input power. Saturn V vehicle stages will use mechani- 
cal feedback actuators with increased electrical input 
power, larger orifices and nozzle sizes, larger torque 
motor wire size, and greater spool driving forces. 

Further considerations for improvement in* system 
reliability led to the development of servovalves with 
improved reliability. The primary objective was to de- 
velop servovalves that have a higher reliability factor 
by increasing electrical input power and minimizing 
sensitivity to fluid contaminatiori . Two approaches were 
taken; in the first approach,. a conventional two stage, 
double nozzle flapper electrohydraulic servovalve with 
mechanical feedback from the spool position was used. 
Methods for optimizing the characteristics effecting the 
reliability are shown, and the test results are given for 
the optimized functional model. In addition to being 
built for use on an electrical feedback actuator, the 
valve was modified to show that it was designed with 
sufficient power for its operation on a mechanical feed- 
back actuator. The second approach is also a two stage 
mechanical feedback servovalve; however, the hydrau- 
I lie amplifier in this case is of the jet pipe receiver de- 
sign. The study began with a theoretical analysis of the 
conventional state-of-the-art servovalves and selection 
* of the jet pipe as offering the greatest potential for high 
reliability. The study then progressed to a series of test 
fixtures and breadboard hardware to evaluate the opti- 
mized design, which through further development re- 
sulted in a functional model. The characteristics of the 
functional model are presented. 


STATUS OF GUIDANCE AND CONTROL METHODS,, 
INSTRUMENTATION, AND TECHNIQUES AS APPLIED 
IN THE APOLLO PROJECT , Walter Haeussermann and 
Robert C. Duncan*, presented at the Lecture Series on 
Orbit Optimization and Advanced Guidance Instrumen- 
tation, Advisory Group for Aeronautical Research and 
Development, North Atlantic Treaty Organization, 
Duesseldorf, Germany, October 21-22, 1964 

Future space missions will place numerous and ex- 
tremely variable requirements on launch vehicles and 
spacecraft. To meet these diverse requirements and pro- 
vide adequate safety for the space crew, guidance and 


control systems must be extremely flexible. As a means 
of presenting these requirements and the methods used to 
meet them, this paper describes the Apollo mission in 
detail. This description touches briefly on the Saturn- 
Apollo space vehicle and goes into deeper detail con- 
cerning the various phases of the lunar landing mission 
(including transit between the earth and moon). A de- 
tailed description of the guidance and control theory, 
schemes, and hardware to be used is presented. 

A short history and description of the development 
of the delta-minimum guidance scheme and its success- 
ful utilization in ballistic missiles are presented. The 
path-adaptive guidance theory used to guide the Saturn 
launch vehicle is discussed in some detail. This theory 
is based on utilization of the calculus of variations to 
solve for an optimum path to achieve a given mission. 
Several means of implementing this theory are being in- 
vestigated to determine which will be most suitable for 
achieving the wide spectrum of expected mission with 
the greatest simplicity and accuracy. The path-adaptive 
guidance schemes permit the guidance system to accept 
a wide range of nonstandard conditions and furnish the 
vehicle with commands to place it on an optimum path. 

The Saturn guidance and control hardware is de- 
scribed in detail, including the factors that determine 
the final selection of components. The guidance and 
control system is designed around a central system of 
basic components, principally a digital computer and a 
data adapter. The subsystems, such as the attitude con- 
trol and the guidance subsystems, all operate in con- 
junction with this central system. For unprecedented 
reliability, redundancy of components and circuits is 
applied almost throughout the entire guidance and con- 
trol system. 

Flight simulation for verification of component and 
system design is also discussed, in its early stages, this 
simulation is accomplished by simulating the components 
themselves. As development progresses, actual com- 
ponents are placed in the flight simulation. 

Prelaunch operations are discussed briefly, followed 
by a discussion of the operations that occur during earth 
orbit. One of these operations is the checkout of the 
onboard guidance and control equipment. Injection of 
the Saturn-Apollo space vehicle into a translunar conic 
will not be initiated until the checkout operations have 


*Manned Spacecraft Center 
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been successfully completed. A short description of the 
transposition and docking maneuver, which occurs after 
injection, is given. 

Following a brief description of the Apollo space- 
craft and the attached lunar excursion module (LEM), a 


general description of the guidance theory that governs 
their motion is presented. This theory embodies statisti- 
cally optimum navigation techniques for updating the 
knowledge of position and velocity in space. The guid- 
ance and control hardware used in the Apollo spacecraft 
and the LEM is described in detail. 



